CN109696705B - Deep sea OBS up-and-down wave field imaging data correction method and device - Google Patents
Deep sea OBS up-and-down wave field imaging data correction method and device Download PDFInfo
- Publication number
- CN109696705B CN109696705B CN201811547866.9A CN201811547866A CN109696705B CN 109696705 B CN109696705 B CN 109696705B CN 201811547866 A CN201811547866 A CN 201811547866A CN 109696705 B CN109696705 B CN 109696705B
- Authority
- CN
- China
- Prior art keywords
- imaging data
- wave field
- going
- time
- field imaging
- 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
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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Oceanography (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method and a device for correcting imaging data of an up-going wave field and a down-going wave field of deep sea OBS (on-board diagnostics), wherein the method comprises the following steps: obtaining imaging data of an up-going wave field and imaging data of a down-going wave field of the deep sea OBS, and assuming that multiples of more than two orders of the up-going wave field and the down-going wave field are basically eliminated; dividing the up-going wave field imaging data and the down-going wave field imaging data into at least two time windows respectively; performing cross-correlation operation on the uplink wave field imaging data and the downlink wave field imaging data in the same time window to obtain cross-correlation coefficients in each time window; determining the time shift amount corresponding to the maximum value of the cross correlation coefficient in each time window; and determining an imaging data correction value according to the time shift amount, and correcting the imaging data of the deep sea OBS up-and-down wave field. The invention obtains the up-going and down-going wave field imaging data which can be accurately matched, and solves the problems that the up-going and down-going wave field imaging data can not be completely matched and the final imaging result is inaccurate due to different offset speeds and noise of the data.
Description
Technical Field
The invention relates to the technical field of seismic data processing, in particular to a method and a device for correcting imaging data of an OBS (on-board diagnostics) up-and-down wave field in deep sea.
Background
The seismic exploration method comprises the technical methods of single-channel earthquake, multi-channel earthquake, high-resolution two-dimensional and three-dimensional earthquake, Ocean Bottom Seismograph (OBS) and the like. In recent years, geophysicists at home and abroad have conducted a series of studies on the above-mentioned technology from different sides and have achieved some results. Among them, the OBS seismic prospecting technology has become the mainstream seismic prospecting method with the advantage that it can provide rich converted shear wave information, plays an important role in marine prospecting, and is regarded by the industry. The OBS seismic exploration technology is a method for receiving and recording seismic waves generated by natural earthquakes and artificial earthquakes by utilizing a plurality of seismographs arranged on the sea bottom and obtaining a seabed geological structure through tomography or migration imaging, overcomes the troubles of conventional offshore reflection seismic exploration on adverse factors such as environment, cables and the like, adopts a unique arrangement mode, can obtain wave field information of a water component, a vertical component and two horizontal components of a hydrophone, and has important significance for improving stratum resolution by obtaining seismic wave information of various effective waves such as abundant longitudinal waves, converted transverse waves and the like.
In OBS seismic data acquisition, a shot point and a demodulator probe are positioned on different datum planes, the shot point is positioned several meters below the sea surface, the demodulator probe is positioned on the seabed and the demodulator probe is sparsely arranged, in the prior art, seismic wave upgoing wave field data are generally utilized, and then tomography imaging or migration imaging is carried out to obtain imaging data for analyzing the seabed geological structure. However, for deep sea OBS seismic data, the illumination of the up-going wave field is poor, which results in poor shallow imaging effect and serious acquisition footprint. The mirror migration technique based on the down-going wavefield data solves these problems well because the down-going wavefield has a broader illumination. However, the imaging effect of the mid-depth layer of the down-going wave field is not as good as that of the up-going wave field, and in order to solve this problem, the up-going wave field data and the down-going wave field data are imaged separately, and then the up-going wave field imaging data and the down-going wave field imaging data are added to form the final imaging result.
However, when the up-and-down wave field imaging data are added, the up-and-down wave field imaging data cannot be completely matched due to the possible difference of offset imaging speeds and the possible existence of noise in the data, which is mainly reflected in phase mismatch, and the final imaging result obtained after the addition is inaccurate. There is no method available in the prior art that can effectively correct deep sea OBS up-and-down wavefield imaging data and the present invention seeks to propose a solution to this problem.
Disclosure of Invention
The embodiment of the invention provides a method for correcting imaging data of an up-going and down-going wave field of a deep sea OBS, which is used for correcting the imaging data of the up-going and down-going wave field of the deep sea OBS to obtain the imaging data of the up-going and down-going wave field which can be accurately matched, so that an accurate imaging result of a submarine geological structure is obtained, and the problems that the imaging data of the up-going and down-going wave fields cannot be completely matched and the final imaging result is inaccurate due to different offset speeds and noise in the data are solved, and the method comprises the following:
obtaining imaging data of an up-going wave field and imaging data of a down-going wave field of the deep sea OBS, wherein the imaging data of the up-going wave field and the imaging data of the down-going wave field are imaging data generated by different receiving modes of the same seismic response;
dividing the up-going wave field imaging data and the down-going wave field imaging data into at least two time windows respectively;
performing cross-correlation operation on the uplink wave field imaging data and the downlink wave field imaging data in the same time window to obtain cross-correlation coefficients in each time window;
determining the time shift amount corresponding to the maximum value of the cross correlation coefficient in each time window;
and determining an imaging data correction value according to the time shift amount, and correcting the imaging data of the deep sea OBS up-and-down wave field.
The embodiment of the invention provides a deep sea OBS up-and-down wave field imaging data correction device, which is used for correcting deep sea OBS up-and-down wave field imaging data to obtain up-and-down wave field imaging data capable of being accurately matched, so that an accurate seabed geological structure imaging result is obtained, and the problems that the up-and-down wave field imaging data cannot be completely matched and the final imaging result is inaccurate due to different offset speeds and data noise are solved, the device comprises:
the data acquisition module is used for acquiring imaging data of an up-going wave field and imaging data of a down-going wave field of the deep sea OBS, wherein the imaging data of the up-going wave field and the imaging data of the down-going wave field are imaging data generated by different receiving modes of the same seismic response;
the time window division module is used for respectively dividing the uplink wave field imaging data and the downlink wave field imaging data into at least two time windows;
the data operation module is used for performing cross-correlation operation on the uplink wave field imaging data and the downlink wave field imaging data in the same time window to obtain cross-correlation coefficients in each time window;
the time shift amount determining module is used for determining the time shift amount corresponding to the maximum value of the cross correlation coefficient in each time window;
and the data correction module is used for determining an imaging data correction value according to the time shift amount and correcting the imaging data of the deep sea OBS up-and-down wave field.
Compared with the scheme that the up-going wave field imaging data and the down-going wave field imaging data are directly added to serve as the final imaging result in the prior art, the up-going wave field imaging data and the down-going wave field imaging data are respectively divided into at least two time windows by obtaining the deep sea OBS up-going wave field imaging data and the down-going wave field imaging data, cross-correlation operation is carried out on the up-going wave field imaging data and the down-going wave field imaging data in the same time window to obtain the cross-correlation coefficient in each time window, and then the time shift amount corresponding to the maximum value of the cross-correlation coefficient is determined in each time window; and determining an imaging data correction value according to the time shift amount, and correcting the deep sea OBS up-going and down-going wave field imaging data to obtain up-going and down-going wave field imaging data which can be accurately matched, and further adding the corrected up-going and down-going wave field imaging data to obtain an accurate final imaging result, so that the problems that the up-going and down-going wave field imaging data cannot be completely matched and the final imaging result is inaccurate due to different offset speeds and noise in the data are solved.
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 or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
FIG. 1 is a schematic diagram of a deep sea OBS up-down wave field imaging data correction method in an embodiment of the invention;
FIG. 2 is a data matching result before the imaging data of the deep sea OBS up-and-down wave field is corrected in the embodiment of the invention;
FIG. 3 is a data matching result after the imaging data of the deep sea OBS up-and-down wave field is corrected according to the embodiment of the invention;
FIG. 4 is a structural diagram of a deep sea OBS up-and-down wave field imaging data correction device in the embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
In order to correct imaging data of an up-and-down wave field of a deep sea OBS, obtain imaging data of the up-and-down wave field which can be accurately matched, obtain an accurate imaging result of a seabed geological structure, and overcome the problems that the imaging data of the up-and-down wave field cannot be completely matched and a final imaging result is inaccurate due to different migration speeds and noise in the data, an embodiment of the invention provides a method for correcting the imaging data of the up-and-down wave field of the deep sea OBS, and as shown in fig. 1, the method can comprise the following steps:
101, obtaining imaging data of an up-going wave field and imaging data of a down-going wave field of deep sea OBS, wherein the imaging data of the up-going wave field and the imaging data of the down-going wave field are imaging data generated by different receiving modes of the same seismic response;
102, dividing the uplink wave field imaging data and the downlink wave field imaging data into at least two time windows respectively;
103, performing cross-correlation operation on the uplink wave field imaging data and the downlink wave field imaging data in the same time window to obtain cross-correlation coefficients in each time window;
104, determining a time shift amount corresponding to the maximum value of the cross correlation coefficient in each time window;
and 105, determining an imaging data correction value according to the time shift amount, and correcting the imaging data of the up-and-down wave field of the deep sea OBS.
As shown in fig. 1, in the embodiment of the present invention, deep sea OBS upgoing wave field imaging data and downgoing wave field imaging data are obtained, the upgoing wave field imaging data and the downgoing wave field imaging data are respectively divided into at least two time windows, cross-correlation operation is performed on the upgoing wave field imaging data and the downgoing wave field imaging data in the same time window to obtain cross-correlation coefficients in the time windows, and then time shift amounts corresponding to maximum values of the cross-correlation coefficients are determined in the time windows; and determining an imaging data correction value according to the time shift amount, and correcting the deep sea OBS up-going and down-going wave field imaging data to obtain up-going and down-going wave field imaging data which can be accurately matched, and further adding the corrected up-going and down-going wave field imaging data to obtain an accurate final imaging result, so that the problems that the up-going and down-going wave field imaging data cannot be completely matched and the final imaging result is inaccurate due to different offset speeds and noise in the data are solved.
In specific implementation, deep sea OBS (on-board diagnostics) upgoing wave field imaging data and downgoing wave field imaging data are obtained, wherein the upgoing wave field imaging data and the downgoing wave field imaging data are imaging data generated by different receiving modes of the same seismic response;
in the embodiment, the acquired imaging data of the deep sea OBS upgoing wave field and the imaging data of the downgoing wave field are preprocessed, and upgoing wave fields and downgoing wave fields are eliminated by multiples of more than two orders. There are many methods for data preprocessing, such as: gain filtering processing, signal-to-noise ratio processing or signal enhancement processing. In this embodiment, the gain filtering processing, the signal-to-noise ratio improving processing, and the signal enhancing processing are all the prior art, and those skilled in the art can understand the gain filtering processing method by looking up the data, and the present invention is not described in detail. It will be understood by those skilled in the art that the foregoing data preprocessing is an exemplary illustration, and various data preprocessing methods can be determined according to the requirements, and all the related modifications are within the scope of the present invention.
In the embodiment, the absolute value is obtained from the acquired imaging data of the deep sea OBS up-going wave field and the imaging data of the down-going wave field, specifically, the imaging data X of the deep sea OBS up-going wave field is obtainedtAnd deep sea OBS down going wave field imaging data YtRespectively calculating absolute values to obtain up-going wave field imaging data | Xt| and down-going wavefield imaging data | YtAnd the data length is N.
In specific implementation, the up-going wave field imaging data and the down-going wave field imaging data are divided into at least two time windows respectively.
In an embodiment, the up-going wave field imaging data and the down-going wave field imaging data are respectively divided into at least two time windows, the time windows are respectively divided into N time windows, and the data length in each time window is N1To obtain | Xt1|，|Xt2|，...，|Xti|，...，|XtnI and Yt1|，|Yt2|，...，|Yti|，...，|YtnL, wherein adjacent time windows can be overlapped by a first set time length N2The first set time length is less than the time window length, i.e. N2<N1。
In specific implementation, performing cross-correlation operation on the uplink wave field imaging data and the downlink wave field imaging data in the same time window to obtain cross-correlation coefficients in each time window; the inventor finds that when the up-and-down wave field imaging data are added, the up-and-down wave field imaging data cannot be completely matched due to different offset imaging speeds and noise of the data, which is mainly reflected in phase mismatch, and the final imaging result obtained after the addition is inaccurate. Therefore, the embodiment of the invention provides a cross-correlation dynamic matching technology based on time window sliding on the assumption that seismic traces of an up-going wave field and a down-going wave field at the same position have stronger similarity.
Performing cross-correlation operation on the uplink wave field imaging data and the downlink wave field imaging data in the same time window according to the following formula:
wherein R isi(τi) Cross correlation coefficient, τiFor the amount of time shift, | XtiL is the up-going wavefield imaging data,
for time-shifted down-going wave field imaging data, N1Is the length of data within the time window.
Thus, for n time windows, n sets of cross-correlation coefficients, R, can be obtained1(τ1)，R2(τ2)，...，Ri(τi)，...，Rn(τn)。
In specific implementation, the time shift amount corresponding to the maximum value of the cross-correlation coefficient is determined in each time window
In the embodiment, n groups of cross correlation coefficients, R, are obtained1(τ1)，R2(τ2)，...，Ri(τi)，...，Rn(τn) Thereafter, the maximum value of each set of cross-correlation coefficients and the time shift amount corresponding to the maximum value are respectively obtained, for example, for the i-th set of cross-correlation coefficients Ri(τi) Calculating Ri(τi) Maximum value of Rimax(τimax) And the amount of time shift τ corresponding to the maximumimaxObtaining R by calculating the maximum value of n groups of cross correlation coefficients and the time shift corresponding to the maximum value1max(τ1max)，R2max(τ2max)，...，Rimax(τimax)，...，Rnmax(τnmax) And τ1max，τ2max，...，τimax，...，τnmax。
In an embodiment, if the determined time shift amount corresponding to the maximum value of the cross-correlation coefficient is greater than a first set threshold, the time shift amount corresponding to the maximum value of the cross-correlation coefficient is adjusted to the first set threshold.
In specific implementation, according to the time shift amount, an imaging data correction amount is determined, and imaging data of the deep sea OBS up-and-down wave field are corrected.
In the embodiment, the time shift τ corresponding to the maximum value of each obtained set of cross-correlation coefficients is first obtained1max，τ2max，...，τimax，...，τnmaxInterpolation is carried out to obtain the uplinkAnd correcting the imaging data of the same data length of the wave field imaging data or the downlink wave field imaging data. For example, for the up-going wave field imaging data or the down-going wave field imaging data having a data length of N, the data length of the imaging data correction amount obtained by interpolation is also N.
In the embodiment, if the adjacent time windows overlap for the first set time length, for the region where the time windows overlap, the midpoint of the time window overlap in the imaging data correction amount is set to 0, and then a slope is set, and transition processing is performed between the value of the imaging data correction amount at the non-overlapping portion of the time windows and the midpoint value of the time window overlap, which is 0.
In the embodiment, the smoothing process may be performed on the imaging data correction amount after the transition process is performed. In the embodiment, the smoothing process is a prior art, and a person skilled in the art can understand the transition processing method by looking up the data, and the invention is not described in detail.
In an embodiment, the up-going wave field imaging data is time-shifted according to an imaging data correction amount to obtain corrected up-going wave field imaging data matched with the down-going wave field imaging data, or the down-going wave field imaging data is time-shifted according to an imaging data correction amount to obtain corrected down-going wave field imaging data matched with the up-going wave field imaging data. Therefore, accurate final imaging results can be obtained after the corrected deep sea OBS up-and-down wave field imaging data are added.
The following provides a specific embodiment, which illustrates a specific application of the deep sea OBS up-and-down wave field imaging data correction method in the embodiment of the present invention. And taking the red sea OBS up-and-down wave field imaging data as an example, correcting the red sea OBS up-and-down wave field imaging data to obtain the up-and-down wave field imaging data which can be accurately matched. First, imaging data of an up-going wave field and imaging data of a down-going wave field of the deep sea OBS are acquired, as shown in fig. 2, the left side in the figure is the imaging data of the up-going wave field, and the right side is the imaging data of the down-going wave field. And respectively carrying out gain filtering processing on the uplink wave field data and the downlink wave field data, improving the signal-to-noise ratio and enhancing the signals, and then dividing the uplink wave field imaging data and the downlink wave field imaging data into 20 time windows, wherein the size of each time window is 250ms, and the time windows are overlapped for 50 ms. The cross-correlation process is performed in each time window, and the first set threshold is set to 20 ms. And determining the time shift amount of each time window according to the maximum cross correlation coefficient and the corresponding maximum time shift amount, carrying out interpolation processing on the time shift amount to obtain an imaging data correction amount, and carrying out transition processing on the imaging data correction amount by using a slope of 25ms at the overlapping part of the time windows. And smoothing the imaging data correction amount after the transition processing, thereby obtaining a final imaging data correction amount. Time shifting is carried out on the up-going wave field imaging data according to the imaging data correction value to obtain corrected up-going wave field imaging data matched with the down-going wave field imaging data, as shown in fig. 3, the up-going wave field imaging data on the left side in fig. 3 can be accurately matched with the down-going wave field imaging data on the right side, therefore, the data are corrected by using the deep sea OBS up-going wave field imaging data correction method, the accurately matched up-going wave field imaging data and down-going wave field imaging data can be obtained, the corrected up-going wave field imaging data and down-going wave field imaging data can be added to obtain an accurate final imaging result, and the problems that the up-going wave field imaging data cannot be completely matched and the final imaging result is inaccurate due to the.
Based on the same inventive concept, the embodiment of the invention also provides a deep sea OBS up-and-down wave field imaging data correction device, which is described in the following embodiment. Because the principles of solving the problems are similar to the method for correcting the imaging data of the up-going and down-going wave fields of the deep sea OBS, the implementation of the device can be referred to the implementation of the method, and repeated details are not repeated.
Fig. 4 is a block diagram of a deep sea OBS up-and-down wave field imaging data correction apparatus according to an embodiment of the present invention, as shown in fig. 4, the apparatus includes:
the data acquisition module 401 is configured to acquire imaging data of an up-going wave field and imaging data of a down-going wave field of the deep sea OBS, where the imaging data of the up-going wave field and the imaging data of the down-going wave field are imaging data generated by different receiving modes of the same seismic response;
a time window dividing module 402, configured to divide the uplink wave field imaging data and the downlink wave field imaging data into at least two time windows respectively;
the data operation module 403 is configured to perform cross-correlation operation on the uplink wave field imaging data and the downlink wave field imaging data in the same time window to obtain cross-correlation coefficients in each time window;
a time shift amount determining module 404, configured to determine, in each time window, a time shift amount corresponding to a maximum value of the cross correlation coefficient;
and the data correction module 405 is configured to determine an imaging data correction amount according to the time shift amount, and correct the imaging data of the up-and-down wave field of the deep sea OBS.
In one embodiment, adjacent time windows overlap for a first set length of time.
In one embodiment, the time shift amount determination module 404 is further configured to: determining the time shift amount corresponding to the maximum value of the cross-correlation coefficient, including: and if the determined time shift amount corresponding to the maximum value of the cross-correlation coefficient is larger than a first set threshold, adjusting the time shift amount corresponding to the maximum value of the cross-correlation coefficient to be the first set threshold.
In one embodiment, the data correction module 405 is further configured to:
interpolating the time shift amount to obtain an imaging data correction amount with the same data length as the up-going wave field imaging data or the down-going wave field imaging data;
time shifting the up-going wave field imaging data according to the imaging data correction amount to obtain corrected up-going wave field imaging data matched with the down-going wave field imaging data, or
And time shifting the downlink wave field imaging data according to the imaging data correction value to obtain corrected downlink wave field imaging data matched with the uplink wave field imaging data.
In summary, in the embodiment of the present invention, the imaging data of the up-going wave field and the imaging data of the down-going wave field of the deep sea OBS are obtained, the imaging data of the up-going wave field and the imaging data of the down-going wave field are respectively divided into at least two time windows, the imaging data of the up-going wave field and the imaging data of the down-going wave field in the same time window are subjected to cross-correlation operation to obtain the cross-correlation coefficient in each time window, and then the time shift corresponding to the maximum value of the cross-correlation coefficient is determined in each time; and determining an imaging data correction value according to the time shift amount, and correcting the deep sea OBS up-going and down-going wave field imaging data to obtain up-going and down-going wave field imaging data which can be accurately matched, and further adding the corrected up-going and down-going wave field imaging data to obtain an accurate final imaging result, so that the problems that the up-going and down-going wave field imaging data cannot be completely matched and the final imaging result is inaccurate due to different offset speeds and noise in the data are solved.
As will be appreciated by one skilled in the art, 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 flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams 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 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 deep sea OBS up-and-down wave field imaging data correction method is characterized by comprising the following steps:
obtaining imaging data of an up-going wave field and imaging data of a down-going wave field of the deep sea OBS, wherein the imaging data of the up-going wave field and the imaging data of the down-going wave field are imaging data generated by different receiving modes of the same seismic response;
dividing the up-going wave field imaging data and the down-going wave field imaging data into at least two time windows respectively;
performing cross-correlation operation on the uplink wave field imaging data and the downlink wave field imaging data in the same time window according to the following formula to obtain cross-correlation coefficients in the time windows:
wherein R isi(τi) Cross correlation coefficient, τiFor the amount of time shift, | XtiL is the up-going wavefield imaging data,
for time-shifted down-going wave field imaging data, N1Is the length of data within the time window;
determining the time shift amount corresponding to the maximum value of the cross correlation coefficient in each time window;
and determining an imaging data correction value according to the time shift amount, and correcting the imaging data of the deep sea OBS up-and-down wave field.
2. The method of claim 1, wherein adjacent time windows overlap for a first set length of time.
3. The method of claim 1, wherein determining the amount of time shift corresponding to the maximum value of the cross-correlation coefficient within each time window comprises: and if the determined time shift amount corresponding to the maximum value of the cross-correlation coefficient is larger than a first set threshold, adjusting the time shift amount corresponding to the maximum value of the cross-correlation coefficient to be the first set threshold.
4. The method of claim 1, wherein determining a correction amount for imaging data based on the amount of time shift to correct for deep sea OBS up and down going wavefield imaging data comprises:
interpolating the time shift amount to obtain an imaging data correction amount with the same data length as the up-going wave field imaging data or the down-going wave field imaging data;
time shifting the up-going wave field imaging data according to the imaging data correction value to obtain corrected up-going wave field imaging data matched with the down-going wave field imaging data, or
And time shifting the downlink wave field imaging data according to the imaging data correction value to obtain corrected downlink wave field imaging data matched with the uplink wave field imaging data.
5. A device that deep sea OBS up-and-down wave field imaging data matches, characterized in that includes:
the data acquisition module is used for acquiring imaging data of an up-going wave field and imaging data of a down-going wave field of the deep sea OBS, wherein the imaging data of the up-going wave field and the imaging data of the down-going wave field are imaging data generated by different receiving modes of the same seismic response;
the time window division module is used for respectively dividing the uplink wave field imaging data and the downlink wave field imaging data into at least two time windows;
the data operation module is used for performing cross-correlation operation on the uplink wave field imaging data and the downlink wave field imaging data in the same time window according to the following formula to obtain cross-correlation coefficients in the time windows:
wherein R isi(τi) Cross correlation coefficient, τiFor the amount of time shift, | XtiL is the up-going wavefield imaging data,
for time-shifted down-going wave field imaging data, N1Is the length of data within the time window;
the time shift amount determining module is used for determining the time shift amount corresponding to the maximum value of the cross correlation coefficient in each time window;
and the data correction module is used for determining an imaging data correction value according to the time shift amount and correcting the imaging data of the deep sea OBS up-and-down wave field.
6. The apparatus of claim 5, wherein adjacent time windows overlap for a first set length of time.
7. The apparatus of claim 5, wherein the time shift amount determination module is further to: determining the time shift amount corresponding to the maximum value of the cross-correlation coefficient, including: and if the determined time shift amount corresponding to the maximum value of the cross-correlation coefficient is larger than a first set threshold, adjusting the time shift amount corresponding to the maximum value of the cross-correlation coefficient to be the first set threshold.
8. The apparatus of claim 5, wherein the data correction module is further to:
interpolating the time shift amount to obtain an imaging data correction amount with the same data length as the up-going wave field imaging data or the down-going wave field imaging data;
time shifting the up-going wave field imaging data according to the imaging data correction value to obtain corrected up-going wave field imaging data matched with the down-going wave field imaging data, or
And time shifting the downlink wave field imaging data according to the imaging data correction value to obtain corrected downlink wave field imaging data matched with the uplink wave field imaging 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 method of any one of claims 1 to 4 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 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811547866.9A CN109696705B (en) | 2018-12-18 | 2018-12-18 | Deep sea OBS up-and-down wave field imaging data correction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811547866.9A CN109696705B (en) | 2018-12-18 | 2018-12-18 | Deep sea OBS up-and-down wave field imaging data correction method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109696705A CN109696705A (en) | 2019-04-30 |
| CN109696705B true CN109696705B (en) | 2020-11-06 |
Family
ID=66231840
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811547866.9A Active CN109696705B (en) | 2018-12-18 | 2018-12-18 | Deep sea OBS up-and-down wave field imaging data correction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109696705B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110824551B (en) * | 2019-11-14 | 2020-12-11 | 自然资源部第二海洋研究所 | OBS accurate positioning method based on acoustic ranging and multi-beam terrain |
| CN111239827B (en) * | 2020-03-09 | 2021-07-30 | 吉林大学 | Three-dimensional seismic data multiple suppression method based on local similarity coefficient |
| CN111856580B (en) * | 2020-07-27 | 2022-02-15 | 广州海洋地质调查局 | Ultra-long offset OBS data first-arrival wave energy enhancement method and processing terminal |
| CN112255685B (en) * | 2020-09-28 | 2022-03-01 | 广州海洋地质调查局 | OBS and sea surface streamer seismic data combined imaging method and processing terminal |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4935903A (en) * | 1989-05-30 | 1990-06-19 | Halliburton Geophysical Services, Inc. | Reinforcement of surface seismic wavefields |
| US4937793A (en) * | 1989-05-30 | 1990-06-26 | Halliburton Geophysical Services, Inc. | Processing method for marine seismic surveying utilizing dual streamers |
| US8345510B2 (en) * | 2008-06-02 | 2013-01-01 | Pgs Geophysical As | Method for aquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source(s) |
| CN102890290B (en) * | 2012-09-25 | 2016-01-06 | 中国石油天然气股份有限公司 | Prestack depth migration method under a kind of relief surface condition |
| CN105572735B (en) * | 2014-10-17 | 2018-02-02 | 中国石油天然气股份有限公司 | A kind of method and device for improving pre-stack depth migration imaging precision |
| US11092709B2 (en) * | 2016-11-17 | 2021-08-17 | Saudi Arabian Oil Company | Use of wavelet cross-correlation for virtual source denoising |
| CN106772616A (en) * | 2016-12-20 | 2017-05-31 | 中国石油天然气股份有限公司 | A kind of processing method and processing device of seismic imaging |
-
2018
- 2018-12-18 CN CN201811547866.9A patent/CN109696705B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109696705A (en) | 2019-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109696705B (en) | Deep sea OBS up-and-down wave field imaging data correction method and device | |
| Shen et al. | Full-waveform inversion: The next leap forward in subsalt imaging | |
| US8792299B2 (en) | Method and device for processing seismic data | |
| US8451682B2 (en) | Method and apparatus for deghosting seismic data | |
| EP3259619B1 (en) | Method and system of processing seismic data by providing surface aperture common image gathers | |
| CN105093301B (en) | The generation method and device of common imaging point angle of reflection angle gathers | |
| CN111025386B (en) | Vertical and horizontal wave separation method without separation false image | |
| CN114779340B (en) | Deep-towed high-resolution multi-channel seismic data fine processing method | |
| Guo et al. | Becoming effective velocity-model builders and depth imagers, Part 2—The basics of velocity-model building, examples and discussions | |
| CN111781635B (en) | Seabed four-component elastic wave Gaussian beam depth migration method and device | |
| US8380440B2 (en) | 3D residual binning and flatness error correction | |
| US10338248B2 (en) | Method for processing seismic images | |
| Zhang et al. | WARpick: A 3D warping tool for seismic gather flattening | |
| US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
| Zdraveva et al. | Improving interpretability of seismic images with directional image partitions and model-based techniques | |
| Perez et al. | Warping prestack imaged data to improve stack quality and resolution | |
| Alwon et al. | Broadband acquisition and processing for high-resolution imaging in the eastern Gulf of Mexico | |
| CN116027403A (en) | Imaging method and device based on wave field decomposition | |
| CN115598696A (en) | Water body velocity model quality control method and device, electronic equipment and readable storage medium | |
| CN115542387A (en) | Method and system for pre-stack migration to produce flare angle migration gathers | |
| Fried et al. | Removing distortions caused by water velocity variations: Water velocity determination | |
| Thiel et al. | 2d Acoustic Full Waveform Inversion of Submarine Salt Layer Using Dual Sensor Streamer Data | |
| CN115598698A (en) | Pre-stack time migration method, device, equipment and medium of downlink wave seismic data | |
| CN111505711A (en) | 3D refracted wave field migration imaging method | |
| CN112946744A (en) | Least square offset imaging method and system based on dynamic time difference warping |
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 |