CN114185095B - Method for suppressing multiple waves of three-dimensional plane wave domain seismic data - Google Patents
Method for suppressing multiple waves of three-dimensional plane wave domain seismic data Download PDFInfo
- Publication number
- CN114185095B CN114185095B CN202111461899.3A CN202111461899A CN114185095B CN 114185095 B CN114185095 B CN 114185095B CN 202111461899 A CN202111461899 A CN 202111461899A CN 114185095 B CN114185095 B CN 114185095B
- Authority
- CN
- China
- Prior art keywords
- domain
- dimensional
- tau
- seismic data
- data
- 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
- 238000000034 method Methods 0.000 title claims abstract description 64
- 230000001629 suppression Effects 0.000 claims abstract description 19
- 238000006243 chemical reaction Methods 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 14
- 230000009466 transformation Effects 0.000 claims description 71
- 239000011159 matrix material Substances 0.000 claims description 45
- 230000005855 radiation Effects 0.000 claims description 36
- 238000001514 detection method Methods 0.000 claims description 29
- 238000013507 mapping Methods 0.000 claims description 20
- 230000001131 transforming effect Effects 0.000 claims description 11
- 238000013506 data mapping Methods 0.000 claims description 6
- 238000003491 array Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 14
- 238000013016 damping Methods 0.000 description 6
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000000844 transformation 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
-
- 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/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for suppressing multiple waves of three-dimensional plane wave domain seismic data, which comprises the following steps: selecting a three-dimensional common shot point gather from the time-space domain seismic data, performing two-time three-dimensional Tau-p conversion, and performing multiple suppression processing to obtain multiple data of a time domain plane wave domain; extracting a common wave point ray parameter gather from the multiple wave data of the time domain plane wave domain, and respectively and alternately performing two-time inverse three-dimensional Tau-p conversion and inverse Fourier conversion to obtain multiple wave seismic data of a time-space domain; and performing self-adaptive subtraction operation on the multiple wave array of the time-space domain and the time-space domain seismic data to obtain a result after the multiple wave suppression of the three-dimensional plane wave domain seismic data. The method for suppressing the multiple waves of the three-dimensional plane wave domain seismic data greatly improves the accuracy of seismic exploration, and has important significance for improving the seismic processing level and the seismic imaging accuracy, enriching and developing the seismic data processing theory.
Description
Technical Field
The invention relates to the field of geophysical exploration digital signal processing, in particular to a method for suppressing multiple waves of three-dimensional plane wave domain seismic data.
Background
In marine seismic exploration, strong multiple interference exists in seismic data under the influence of strong reflection coefficients of the sea water surface. The existence of multiple waves complicates the wave field of the seismic wave, causes serious interference on the amplitude and energy of the reflected wave, limits the bandwidth of the seismic data, obscures the knowledge of the subsurface geologic structure, and is a serious obstacle to the processing and interpretation of the seismic data.
Processing three-dimensional data volumes acquired in three dimensions using conventional two-dimensional seismic data processing methods is error-prone, resulting in low accuracy in seismic exploration. With the increasing depth of seismic exploration, the existing two-dimensional exploration has not met the requirement of exploration precision.
Disclosure of Invention
The invention aims to provide a method for suppressing multiple waves of three-dimensional plane wave domain seismic data, which is used for solving the problem of low precision of the existing two-dimensional exploration data.
The invention provides a method for suppressing multiple waves of three-dimensional plane wave domain seismic data, which comprises the following steps:
from time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Selecting a three-dimensional common shot point gather, and obtaining a three-dimensional plane wave Tau-p domain earthquake through twice three-dimensional Tau-p transformationData
For the three-dimensional plane wave Tau-p domain seismic data 
Performing multiple suppression processing to obtain multiple data of time domain plane wave domain +.>
Multiple data from the time domain plane wave domain
Extracting the common wave-detecting point ray parameter gather, respectively and alternately performing two times of inverse three-dimensional Tau-p conversion and inverse Fourier conversion to obtain multiple wave seismic data M (t, x) of a time-space domain r ,y r ,x s ,y s )；
Multiple arrays M (t, x) of the time-space domain r ,y r ,x s ,y s ) And the time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Performing self-adaptive subtraction operation to obtain a result after multiple suppression of the three-dimensional plane wave domain seismic data;
wherein t is the longitudinal time, x r Seismic traces in the inline direction, y r Seismic traces in the crossline direction, x s A seismic source in the direction of the inline line, y s Is a seismic source in the cross-line direction,
for the inline direction detector radiation parameters, < +.>
The wave-detecting point radiation parameters are the wave-detecting point radiation parameters in the transverse measuring line direction; />
Is a longitudinal lineDirectional Source ray parameters, +.>
Is a cross-line direction source ray parameter.
Further, obtaining the three-dimensional plane wave Tau-p domain seismic data through twice three-dimensional Tau-p transformation
The method comprises the following steps:
step A1: first three-dimensional Tau-p transform
From time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Selecting three-dimensional common shot point gather d 1 (t,x r ,y r ) The method comprises the steps of carrying out a first treatment on the surface of the Gathering d of three-dimensional common shot points 1 (t,x r ,y r ) Performing Fourier transform, first three-dimensional Tau-p transform and inverse Fourier transform to obtain seismic data of the wave detection point Tau-p time domain
Seismic data of the wave detection point Tau-p time domain +.>
Time-space domain seismic data D (t, x) placed at corresponding locations r ,y r ,x s ,y s ) In (1) obtaining the seismic data of the wave detection point Tau-p domain +.>
Step A2: second order three-dimensional Tau-p transformation
From the detector points Tau-p domain seismic data
Selecting a three-dimensional common-detection-point ray parameter gather d 2 (t，x s ，y s ) And fourier transformingPerforming a second three-dimensional Tau-p transformation on the data corresponding to the main frequency in the obtained data, and then calculating three-dimensional Tau-p transformation data +.>
Formation of three-dimensional Tau-p domain seismic data in full frequency domain>
From a three-dimensional common-detector-point ray parameter gather d 2 (t，x s ，y s ) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data +.>
Firstly, after inverse Fourier transformation, the seismic data of the three-dimensional plane wave domain are put into the corresponding position +.>
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained>
Wherein t is longitudinal time, and f is frequency; x is x r Is an earthquake channel in the direction of the in-line,
for the radial parameter of the detecting point in the direction of the inline line, y r Is a seismic trace in the direction of a transverse survey line, +. >
Is the radiation parameter of the detector point in the transverse direction, x s Is a seismic source in the direction of the inline>
Is the source ray parameter in the direction of the inline line, y s Is a source in the transverse line direction +.>
Is a cross-line direction source ray parameter.
Further, the step A1 includes the steps of:
step A11: given a known time-space domain seismic data D (t, x r ,y r ,x s ,y s )；
Step A12: from time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Selecting a three-dimensional common shot point gather d 1 (t,x r ,y r )；
Step A13: gathering d of three-dimensional common shot points 1 (t,x r ,y r ) Performing Fourier transform, first three-dimensional Tau-p transform and inverse Fourier transform to obtain seismic data of the wave detection point Tau-p time domain
Wherein t is the longitudinal time, < >>
For the inline direction detector radiation parameters, < +.>
The wave-detecting point radiation parameters are the wave-detecting point radiation parameters in the transverse measuring line direction;
step A14: seismic data of a wave-detecting point Tau-p time domain
Time-space domain seismic data D (t, x) placed at corresponding locations r ,y r ,x s ,y s ) In (1) obtaining the seismic data of the wave detection point Tau-p domain +.>
Specifically, the step a13 includes the steps of:
step A131: three-dimensional common shot point gather d 1 (t,x r ,y r ) Obtaining three-dimensional common shot point gather data df of a frequency domain by utilizing Fourier transformation to the frequency domain 1 (f,x r ,y r ) Wherein f is frequency;
step A132: data df of three-dimensional common shot point gather in frequency domain 1 (f,x r ,y r ) Obtaining Tau-p domain data of the wave-detecting point frequency domain through primary three-dimensional Tau-p conversion
Step A133: tau-p domain data of wave detector frequency domain
Performing inverse Fourier transform to obtain seismic data of the detection point Tau-p time domain>
Specifically, the step A2 includes the steps of:
step A21: from the wave-detecting point Tau-p domain seismic data
Selecting a three-dimensional common-detection-point ray parameter gather d 2 (t，x s ，y s )；
Step A22: gather d of three-dimensional common-detection-point ray parameters 2 (t，x s ，y s ) Obtaining three-dimensional common shot point ray parameter gather data df of a frequency domain through Fourier transformation to the frequency domain 2 (f，x s ，y s )；
Step A23: according to frequencyThree-dimensional co-shot ray parameter gather data df for domain 2 (f，x s ，y s ) Firstly, carrying out a second three-dimensional Tau-p transformation according to the seismic data corresponding to the main frequency of the seismic data to obtain a shot Tau-p frequency domain gather of the seismic data corresponding to the main frequency
Then, three-dimensional Tau-p transformation data of the corresponding seismic data outside the main frequency is calculated according to the three-dimensional Tau-p transformation data
The shot Tau-p frequency domain gather of which the main frequency corresponds to the seismic data is +.>
Three-dimensional Tau-p transform data corresponding to seismic data outside the dominant frequency>
Combining to form three-dimensional Tau-p domain seismic data in full frequency domain>
Step A24: gather d of three-dimensional common-detection-point ray parameters 2 (t，x s ，y s ) Transforming to obtain three-dimensional plane wave domain seismic data 
Three-dimensional Tau-p domain seismic data in full frequency domain are->
Firstly, after inverse Fourier transformation, the seismic data of the three-dimensional plane wave domain are put into the corresponding position +.>
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained>
More specifically, the step a23 includes the steps of:
step A231: selecting three-dimensional common shot point ray parameter gather data df of frequency domain 2 (f，x s ，y s ) Data corresponding to the main frequency in the three-dimensional common shot gather data df0 of the seismic data corresponding to the main frequency 2 (f，x s ，y s ) And for three-dimensional common shot point gather data df0 of which main frequency corresponds to seismic data 2 (f，x s ，y s ) Performing a second three-dimensional Tau-p transformation to obtain a shot Tau-p frequency domain gather of the seismic data with main frequency
Step A232: shot Tau-p frequency domain gather corresponding to seismic data according to main frequency
Calculating a diagonal constraint matrix W of the inline and a diagonal constraint matrix V of the crossline;
step A233: according to the transverse line diagonal constraint matrix W and the transverse line diagonal constraint matrix V, calculating three-dimensional Tau-p conversion data of the corresponding seismic data except the main frequency
Step a234: shot Tau-p frequency domain gather of main frequency corresponding to seismic data
Three-dimensional Tau-p transform data corresponding to seismic data outside the dominant frequency >
Combining to form three-dimensional Tau-p domain seismic data in full frequency domain>
Further, performing multiple suppression processing to obtain multiple data of the time domain plane wave domain
The method comprises the following steps:
step B1: for the three-dimensional plane wave Tau-p domain seismic data
Performing Fourier transform along time direction to obtain seismic data +.>
Step B2: seismic data from the plane wave domain of the frequency domain
Selecting an initial frequency slice from the results of (a)>
For the selected initial frequency slice->
Performing linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slice ∈ ->
Step B3: slicing multiple data frequencies
Multiple data array for composing frequency domain plane wave domain>
Step B4: multiple data array for frequency domain plane wave domain
Performing inverse Fourier transform to obtain multiple data of time domain plane wave domain>
Wherein t is longitudinal time, and f is frequency; x is x r Is an earthquake channel in the direction of the in-line,
for the radial parameter of the detecting point in the direction of the inline line, y r Is a seismic trace in the direction of a transverse survey line, +.>
Is the radiation parameter of the detector point in the transverse direction, x s Is a seismic source in the direction of the inline>
Is the source ray parameter in the direction of the inline line, y s Is a source in the transverse line direction +.>
Is a cross-line direction source ray parameter. />
Specifically, the step B2 includes the steps of:
step B21: seismic data from the plane wave domain of the frequency domain
Selecting a frequency slice as an initial frequency slice +.>
Slicing the initial frequency>
Performing linear mapping to obtain mapped frequency slices +.>
Step B22: slicing the mapped frequency
Performing squaring operation to obtain frequency slice after multiple data mapping>
Step B23: frequency slicing after mapping the multiple data
Inverse linear mapping is carried out to obtain multiple wave data frequency slices +.>
Further, the inverse three-dimensional Tau-p transformation and the inverse Fourier transformation are alternately performed twice to obtain multiple seismic data M (t, x) of the time-space domain r ,y r ,x s ,y s ) The method comprises the following steps:
step C2: for common detection point ray parameter trace set
Obtaining a common-detection-point ray parameter gather m after the first inverse transformation of the frequency domain through the first inverse three-dimensional Tau-p transformation 2 (t,x s ,y s )；
Step C3: first inverse transformed common-detector-point ray parameter trace m for frequency domain 2 (t,x s ,y s ) Corresponding data
Performing the first inverse Fourier transform to obtain a common-detection-point ray parameter gather after the first inverse transform 
And the plurality of common-wave-point radiation parameter gathers after the first inverse transformation are +.>
Composition of three-dimensional Tau-p Domain dataset->
Step C4: from three-dimensional Tau-p domain data sets
Co-shot gather of three-dimensional Tau-p domain data is extracted>
Performing a second inverse three-dimensional Tau-p transformation to obtain a multiple three-dimensional common shot point gather m 4 (t,x r ,y r )；
Step C5: data mf corresponding to three-dimensional common shot point gather of frequency domain 4 (t,x r ,y r ) Performing a second inverse Fourier transform to obtain a three-dimensional common shot point gather m in a time domain 4 (t,x r ,y r ) And three-dimensional common shot point gather m of a plurality of time domains 4 (t,x r ,y r ) Multiple arrays M (t, x) forming a time-space domain r ,y r ,x s ,y s )。
Further, the method for obtaining the result of the multiple suppression of the three-dimensional plane wave domain seismic data comprises the following steps:
assuming that the overall energy of the seismic data after multiple suppression is minimum;
first input time-space domain seismic data D (t, x r ,y r ,x s ,y s )；
Then the time-space domain seismic data D (t, x r ,y r ,x s ,y s ) And multiple array M (t, x) of time-space domain r ,y r ,x s ,y s ) And performing self-adaptive subtraction operation.
Compared with the prior art, the invention has the beneficial effects that:
the invention discloses a method for pressing multiple waves of three-dimensional plane wave domain seismic data, which adopts a three-dimensional seismic method for collecting underground information on a certain area, can know the condition of an underground geological structure from a three-dimensional space in a three-dimensional way, and can reduce the requirement of a computer on the storage of the three-dimensional seismic data to a certain extent by utilizing the compression characteristic of the two-dimensional plane wave domain on the seismic data based on Tau-p transformation; under the condition that the seismic data transversely fluctuates, the seismic data can be effectively compressed by transforming the seismic data into a biplane wave domain, the amount of the compressed seismic data is reduced, the calculated amount can be reduced to a certain extent, the requirement on a computer storage medium is reduced, and the calculation efficiency is improved; the multiple pressing method does not need information of underground medium, and is a complete data driving method. The method for suppressing the multiple waves of the three-dimensional plane wave domain seismic data can provide cross-section, plane and three-dimensional underground geologic structure images, greatly improves the accuracy of seismic exploration, is particularly effective for areas with complicated and changeable underground geologic structures, and has important significance for improving the seismic processing level and the seismic imaging accuracy and enriching and developing the seismic data processing theory.
Drawings
FIG. 1 is a schematic view of an input three-dimensional common shot gather d according to embodiment 1 of the present invention 1 (t,x r ,y r ) Is a display result of (1);
FIG. 2 is a cross-sectional view of a three-dimensional common shot gather d according to embodiment 1 of the present invention 1 (t,x r ,y r ) Wave-detecting point Tau-p time domain seismic data obtained through first Tau-p conversion
Results of (2);
FIG. 3 shows the seismic data of the time domain of the detection points Tau-p provided in embodiment 1 of the present invention
Selecting three-dimensional common-detection-point ray parameter gathers according to keywordsd 2 (t，x s ，y s ) And obtaining the seismic data of the time domain three-dimensional Tau-p transformation after carrying out Fourier transformation, second three-dimensional Tau-p transformation and inverse Fourier transformation>
Results of (2);
FIG. 4 is a plot of seismic data from the dip Tau-p time domain of FIG. 3
After fourier transformation, a frequency slice is selected from the result as the initial frequency slice +.>
Slicing the initial frequency
Mapped frequency slice obtained by linear mapping>
Wherein FIG. 4 (a) is the initial frequency slice before linear mapping +.>
FIG. 4 (b) is a frequency slice after mapping
FIG. 5 is a multiple array of time-space domain provided in embodiment 1 of the present invention
Displaying results in a time-space domain;
FIG. 6 shows the method for providing the space domain seismic data D (t, x) according to embodiment 1 of the invention r ,y r ,x s ,y s ) And multiple array M (t, x) of time-space domain r ,y r ,x s ,y s ) A primary wave result obtained by self-adaptive subtraction operation, whereinThe position indicated by the arrow in fig. 6 is originally numerous.
Detailed Description
The invention provides a method for multiple suppression of three-dimensional plane wave domain seismic data, which comprises the steps of firstly transforming the seismic data into a three-dimensional plane wave domain, and according to the characteristics of dense detection points and sparse shots of the three-dimensional seismic data, transforming the three-dimensional plane wave domain seismic data by adopting a first Tau-p transformation algorithm and a second Tau-p transformation algorithm based on main frequency constraint in sequence, so as to improve the representation precision of the plane wave domain to time space domain seismic data; performing multiple prediction processing by mapping, multiplying and the like on the frequency slices in the plane wave domain; and then inversely converting the result of the multiple prediction back to a time-space domain, and carrying out self-adaptive subtraction to achieve the purposes of multiple suppression and data precision improvement.
Example 1
Embodiment 1 provides a method for multiple suppression of three-dimensional plane wave domain seismic data, comprising the following steps:
step A: two three-dimensional Tau-p transforms
From time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Selecting a three-dimensional common shot point gather, and obtaining seismic data of a three-dimensional plane wave domain through two three-dimensional Tau-p transformations 
The method specifically comprises the following steps:
step A1: first three-dimensional Tau-p transform
From time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Selecting three-dimensional common shot point gather d 1 (t,x r ,y r ) The method comprises the steps of carrying out a first treatment on the surface of the Gathering d of three-dimensional common shot points 1 (t,x r ,y r ) Performing Fourier transform, first three-dimensional Tau-p transform and inverse Fourier transform to obtain seismic data of the wave detection point Tau-p time domain
Seismic data of the wave detection point Tau-p time domain +.>
Time-space domain seismic data D (t, x) placed at corresponding locations r ,y r ,x s ,y s ) In (1) obtaining the seismic data of the wave detection point Tau-p domain +.>
Wherein t is the longitudinal time, x r Seismic traces in the inline direction, y r Seismic traces in the crossline direction, x s A seismic source in the direction of the inline line, y s For seismic source in transverse direction, p xr For the radial parameter of the detecting point in the direction of the longitudinal line, p yr Is in the transverse direction of the line the spot radiation parameters.
The detailed steps are as follows:
step A11: given a known time-space domain seismic data D (t, x r ,y r ,x s ,y s )；
Wherein t is the longitudinal time, x r Seismic traces in the inline direction, y r Seismic traces in the crossline direction, x s A seismic source in the direction of the inline line, y s A seismic source in the direction of a transverse survey line;
step A12: from time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Selecting a three-dimensional common shot point gather d 1 (t,x r ,y r )；
Step A13: gathering d of three-dimensional common shot points 1 (t,x r ,y r ) Performing Fourier transform, first three-dimensional Tau-p transform and inverse Fourier transform to obtain seismic data of the wave detection point Tau-p time domain
Wherein t is the longitudinal time, < >>
For the inline direction detector radiation parameters, < +.>
Is in the transverse direction of the line the spot radiation parameters.
Wherein, the first three-dimensional Tau-p transformation process is performed in the frequency domain, and more detailed steps are as follows:
step A131: three-dimensional common shot point gather d 1 (t,x r ,y r ) Obtaining three-dimensional common shot point gather data df of a frequency domain through Fourier transformation to the frequency domain 1 (f,x r ,y r ) Wherein f is frequency;
step A132: data df of three-dimensional common shot point gather in frequency domain 1 (f,x r ,y r ) Obtaining Tau-p domain data of the wave-detecting point frequency domain through primary three-dimensional Tau-p conversion
Wherein, tau-p domain data of the wave detection point frequency domain
The calculated expression of (2) is:
in the method, in the process of the invention,
tau-p domain data, df, for the frequency domain of the detector 1 Three-dimensional common shot gather df for frequency domain 1 (f,x r ,y r )；
i is an identity matrix; lambda is a damping operator; h represents a conjugate transpose;
Specifically, three-dimensional common shot gather df in frequency domain 1 Operators of inline directions
From a three-dimensional common shot gather d 1 (t，x r ，y r ) Is a seismic trace x in the inline direction r And the radial parameter in the direction of the inline>
The calculation formula is as follows:
wherein x is r Is x r1 ，x r2 ，…，x rn Wherein r is 1 ，r 2 ，…，r n Is that r xSubscript of (2)Serial numbers not only represent seismic traces in the inline direction, but also represent inline direction offset;
Three-dimensional common shot point gather df in frequency domain 1 Operator of crossline direction of (2)
From a three-dimensional common shot gather d 1 (t，x r ，y r ) Is a cross-line direction seismic trace y r And the radiation parameters in the crossline direction +.>
The calculation formula is as follows:
wherein y is r Representing y r1 ，y r2 ，…，y rn Wherein r is 1 ，r 2 ，…，r n Is y r The subscript sequence number in the cross-line direction not only represents the seismic trace in the cross-line direction, but also represents the offset distance in the cross-line direction;
Step A133: tau-p domain data of wave detector frequency domain
Performing inverse Fourier transform to obtain seismic data of the detection point Tau-p time domain>
As a specific embodiment, when three-dimensional common shot point gather d 1 (t，x r ，y r ) When the input record is shown in figure 1, obtaining the seismic data of the wave-detecting point Tau-p time domain through Fourier transformation, first three-dimensional Tau-p transformation and inverse Fourier transformation 
As shown in fig. 2.
Step A14: seismic data of a wave-detecting point Tau-p time domain
Time-space domain seismic data D (t, x) placed at corresponding locations r ,y r ,x s ,y s ) In (1) obtaining the seismic data of the wave detection point Tau-p domain +.>
Step A2: second order three-dimensional Tau-p transformation
From the wave-detecting point Tau-p domain seismic data
Selecting a three-dimensional common-detection-point ray parameter gather d 2 (t，x s ，y s ) After Fourier transformation, performing a second three-dimensional Tau-p transformation on the data corresponding to the main frequency in the obtained data, and then calculating the three-dimensional T of the seismic data corresponding to the main frequencyau-p transformation data->
Formation of three-dimensional Tau-p domain seismic data in full frequency domain>
From a three-dimensional common-detector-point ray parameter gather d 2 (t，x s ，y s ) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data +.>
Firstly, after inverse Fourier transformation, the seismic data of the three-dimensional plane wave domain are put into the corresponding position +.>
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained>
Wherein t is longitudinal time, and f is frequency; x is x r Is an earthquake channel in the direction of the in-line,
for the radial parameter of the detecting point in the direction of the inline line, y r Is a seismic trace in the direction of a transverse survey line, +.>
The wave-detecting point radiation parameters are the wave-detecting point radiation parameters in the transverse measuring line direction; x is x s Seismic source in inline direction，/>
Is the source ray parameter in the direction of the inline line, y s Is a source in the transverse line direction +.>
Is a cross-line direction source ray parameter.
The detailed steps are as follows:
step A21: from the wave-detecting point Tau-p domain seismic data
Selecting a three-dimensional common-detection-point ray parameter gather d 2 (t，x s ，y s )；
As a specific implementation mode, 1, 4 and 5 are selected as keywords, and the seismic data in the Tau-p domain is obtained from the wave detection points
Extracting corresponding dimension according to key words to select a three-dimensional common-detection-point ray parameter gather d 2 (t，x s ，y s ) Seismic data from the wave-detection point Tau-p domain +.>
Extracting the 1 st, 4 th and 5 th three-dimensional common-wave-point ray parameter trace d 2 (t，x s ，y s )。
Step A22: gather d of three-dimensional common-detection-point ray parameters 2 (t，x s ，y s ) Obtaining three-dimensional common shot point ray parameter gather data df of a frequency domain through Fourier transformation to the frequency domain 2 (f，x s ，y s )；
Step A23: three-dimensional common shot point ray parameter gather data df according to frequency domain 2 (f，x s ，y s ) Firstly, carrying out a second three-dimensional Tau-p transformation according to the seismic data corresponding to the main frequency of the seismic data to obtain a shot Tau-p frequency domain gather of the seismic data corresponding to the main frequency
Then, three-dimensional Tau-p transformation data of the corresponding seismic data outside the main frequency is calculated according to the three-dimensional Tau-p transformation data
The shot Tau-p frequency domain gather of which the main frequency corresponds to the seismic data is +.>
Three-dimensional Tau-p transform data corresponding to seismic data outside the dominant frequency >
Combining to form three-dimensional Tau-p domain seismic data in full frequency domain>
The method specifically comprises the following steps:
step A231: selecting three-dimensional common shot point ray parameter gather data df of frequency domain 2 (f，x s ，y s ) Data corresponding to the main frequency in the three-dimensional common shot gather data df0 of the seismic data corresponding to the main frequency 2 (f，x s ，y s ) And for three-dimensional common shot point gather data df0 of which main frequency corresponds to seismic data 2 (f，x s ，y s ) Performing a second three-dimensional Tau-p transformation to obtain a shot Tau-p frequency domain gather of the seismic data with main frequency
Wherein the main frequency corresponds to shot Tau-p frequency domain gather of the seismic data
The calculated expression of (2) is:
in the formula, df0 2taup Shot Tau-p frequency domain gather for main frequency corresponding to seismic data
df0 2 Three-dimensional common shot gather data df0 for primary frequency corresponding seismic data 2 (f，x s ，y s )；
i is an identity matrix; μ is a damping parameter; h represents a conjugate transpose;
Specifically, the dominant frequency corresponds toShot Tau-p frequency domain gather of seismic data
Operator of inline direction->
From a common-detector-point radial parameter gather d 2 (t，x s ，y s ) Seismic source x in the inline direction of (2) s And the radial parameter in the direction of the inline>
Calculating, wherein the calculation formula is as follows; />
Wherein x is s Represents x s1 ，x s2 ，…，x sn Wherein s is 1 ，s 2 ，…，s n Is x s Subscript number, x in (x) s Not only representing the seismic source in the direction of the inline line, but also representing the offset distance in the direction of the inline line;
Shot Tau-p frequency domain gather with main frequency corresponding to seismic data
Operator of crossline direction->
From a common-detector-point radial parameter gather d 2 (t，x s ，y s ) Source y in the crossline direction of (2) s And the radiation parameters in the crossline direction +.>
Calculating, wherein the calculation formula is as follows;
wherein y is s Representing y s1 ，y s2 ，…，y sn Wherein s is 1 ，s 2 ，…，s n Is y s Subscript number, y in (2) s Not only representing the seismic source in the transverse line direction, but also representing the transverse line direction offset;
Step A232: shot Tau-p frequency domain gather corresponding to seismic data according to main frequency
Calculating a diagonal constraint matrix W of the inline and a diagonal constraint matrix V of the crossline;
The detailed calculation process is as follows:
set shot Tau-p frequency domain gather df0 with main frequency corresponding to seismic data 2taup Operators of inline directions
The corresponding three-dimensional common shot point gather data is df0 XS The expression is:
in the method, in the process of the invention,
shot Tau-p frequency domain gather df0 for primary frequency corresponding seismic data 2taup Operators of the inline direction;
i is an identity matrix; μ is a damping parameter; h represents a conjugate transpose;
df0 2 three-dimensional common shot gather data df0 for primary frequency corresponding seismic data 2 (f，x s ，y s )；
Set shot Tau-p frequency domain gather df0 with main frequency corresponding to seismic data 2taup Operator L of crossline direction of (2) ys The corresponding three-dimensional common shot point gather data is df0 YS The expression is:
in the method, in the process of the invention,
shot Tau-p frequency domain gather df0 for primary frequency corresponding seismic data 2taup Operators of the crossline direction of (2);
i is an identity matrix; μ is a damping parameter; h represents a conjugate transpose;
Since the non-diagonal elements of the diagonal constraint matrix W of the inline are all 0, only the matrix W formed by the diagonal elements of the diagonal constraint matrix W of the inline is calculated during calculation ii The expression is:
wherein ε is the stability factor of the inline; i is the number of rows or columns of diagonal positions of the diagonal constraint matrix W of the inline; w (W) ii Is a matrix composed of diagonal elements of a diagonal constraint matrix W of the inline;
Since the non-diagonal elements of the diagonal constraint matrix V of the transverse line are all 0, only the matrix V formed by the diagonal elements of the diagonal constraint matrix V of the transverse line is needed to be calculated during calculation ii The expression is:
wherein ζ is a stabilizing factor of the crossline; i is the number of rows or columns of diagonal positions of the diagonal constraint matrix V of the crossline; v (V) ii Is a matrix composed of diagonal elements of a diagonal constraint matrix V of the transverse measuring lines;
Step A233: according to the transverse line diagonal constraint matrix W and the transverse line diagonal constraint matrix V, calculating three-dimensional Tau-p conversion data of the corresponding seismic data except the main frequency
Wherein the three-dimensional Tau-p transform data of the corresponding seismic data outside the primary frequency
The calculation formula of (2) is as follows:
wherein df1 2taup Three-dimensional Tau-p transform data for seismic data corresponding to frequencies outside of the dominant frequency
df0 2 Three-dimensional for seismic data corresponding to dominant frequenciesCo-shot gather data df0 2 (f，x s ，y s )；
W is the diagonal constraint matrix of the inline; v is the diagonal constraint matrix of the transverse line;
i is an identity matrix; h represents a conjugate transpose;
lambda' is a damping operator in the direction of the inline line; mu' is a damping operator in the transverse line direction;
step a234: shot Tau-p frequency domain gather of main frequency corresponding to seismic data
Three-dimensional Tau-p transform data corresponding to seismic data outside the dominant frequency>
Combining to form three-dimensional Tau-p domain seismic data in full frequency domain>
Step A24: gather d of three-dimensional common-detection-point ray parameters 2 (t，x s ，y s ) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data in full frequency domain are->
Firstly, after inverse Fourier transformation, the seismic data of the three-dimensional plane wave domain are put into the corresponding position +.>
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained>
The detailed steps are as follows:
step A241: gather d of three-dimensional common-detection-point ray parameters 2 (t，x s ，y s ) Transforming to obtain three-dimensional plane wave domain seismic data
Step a242: three-dimensional Tau-p domain seismic data in full frequency domain
Inverse Fourier transform to obtain time domain three-dimensional Tau-p transformed seismic data +.>
Step A243: seismic data transformed in time domain three-dimensional Tau-p
Three-dimensional plane wave domain seismic data +.>
Obtaining three-dimensional plane wave Tau-p domain seismic data
In particular, when three-dimensional common-detector-point ray parameter trace d 2 (t，x s ，y s ) Obtaining the seismic data of the time domain three-dimensional Tau-p transformation through Fourier transformation, secondary three-dimensional Tau-p transformation and inverse Fourier transformation 
The results of (2) are shown in FIG. 3. It can be seen that significant focal points occur in a limited region of the data, such as between 600-1200 of the source ray parameters, indicating that the three-dimensional Tau-p transform compresses the dominant energy of the data into these points, illustrating the compression characteristics of the method on the data.
And (B) step (B): for three-dimensional plane wave Tau-p domain seismic data
Performing multiple suppression processing to obtain multiple data of time domain plane wave domain +.>
The method specifically comprises the following steps:
firstly, for three-dimensional plane wave Tau-p domain seismic data
Fourier transforming in time directionThen selecting the initial frequency slice from the result>
Slicing the selected initial frequency
Performing linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slices
Then frequency slicing multiple wave data
Multiple data array for composing frequency domain plane wave domain>
Finally, multiple data array of frequency domain plane wave domain
Performing inverse Fourier transform to obtain multiple data of time domain plane wave domain>
Wherein t is longitudinal time, and f is frequency; x is x r Is an earthquake channel in the direction of the in-line,
for the radial parameter of the detecting point in the direction of the inline line, y r Is a seismic trace in the direction of a transverse survey line, +.>
Is the radiation parameter of the detector point in the transverse direction, x s Is a seismic source in the direction of the inline>
Is the source ray parameter in the direction of the inline line, y s Is a source in the transverse line direction +.>
Is a cross-line direction source ray parameter.
The more detailed steps are as follows:
step B1: seismic data for three-dimensional plane wave domain
Performing Fourier transform on the array along the time direction to obtain seismic data +.>
Step B2: seismic data from the plane wave domain of the frequency domain
Selecting initial frequency slices
For the selected initial frequency slice->
Performing linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slice ∈ ->
The method specifically comprises the following steps:
step B21: seismic data from the plane wave domain of the frequency domain
Selecting a frequency slice as an initial frequency slice +.>
To the beginningStart frequency slice->
Performing linear mapping to obtain mapped frequency slices +.>
Specifically, it will
And->
Respectively substituting the initial frequency slices +.>
Is->
And->
Obtaining mapped frequency slices +.>
Then the initial frequency slice +.>
And mapped frequency slice->
The calculated expression of (2) is:
in the method, in the process of the invention,
for the inline common uplink parameters, +.>
For the transverse line sharing the uplink ray parameters, +.>
For the inline co-downlink radiation parameters, +. >
The downlink ray parameters are shared for the transverse measuring line; />
For the inline direction detector radiation parameters, < +.>
The wave-detecting point radiation parameters are the wave-detecting point radiation parameters in the transverse measuring line direction; />
For the source ray parameter in the inline direction, +.>
Is a cross-line direction source ray parameter.
As can be seen from fig. 4, the initial frequency slice
And mapped frequency slices
Each of which comprises a plurality of small matrices, wherein each small matrix is a block matrix, and the mapped frequency slices are +.>
Each row of the small matrix of (a) is a common downlink ray parameter P d Trace set, each row of the small matrix is inline common downlink ray parameter +.>
Each small matrix occupies a common downstream ray parameter P d Is a crossline common downlink ray parameter +.>
A numerical value; mapped frequency slice->
Each column of the small matrix of (a) is a common uplink radiation parameter P o Trace set, each row of the small matrix is inline common uplink ray parameter +.>
Each small matrix is a block matrix and occupies the common uplink ray parameter P o Is a crossline common uplink radiation parameter +.>
Numerical values.
Step B22: slicing the mapped frequency
Performing squaring operation to obtain frequency slice after multiple data mapping>
Wherein for each mapped frequency slice
Multiplication operation is carried out, the multiplication operation of the small matrix represents convolution operation of an uplink ray parameter gather and a downlink ray parameter gather time-space domain, and frequency slices after multiplication to obtain multiple wave data mapping are obtained >
The calculation formula is as follows:
in the method, in the process of the invention,
representing each mapped frequency slice; />
Representing the frequency slice after the multiple wave data mapping; />
For the inline common uplink parameters, +.>
The uplink ray parameters are shared for the transverse measurement line,
for the inline co-downlink radiation parameters, +.>
The downlink ray parameters are shared for the transverse measuring line;
step B23: frequency slicing after multiple data mapping
Inverse linear mapping is carried out to obtain multiple wave data frequency slices +.>
Step B3: slicing multiple data frequencies
Multiple data array for composing frequency domain plane wave domain>
Wherein, the multiple wave data frequencyRate slice
The calculated expression of (2) is:
in the method, in the process of the invention,
representing multiple data frequency slices, ++>
For the inline direction detector radiation parameters, < +.>
The wave-detecting point radiation parameters are the wave-detecting point radiation parameters in the transverse measuring line direction; />
For the source ray parameter in the inline direction, +.>
The method is characterized in that the method is a transverse survey line direction seismic source ray parameter; />
For the inline common uplink parameters, +.>
For the transverse line sharing the uplink ray parameters, +.>
For the inline co-downlink radiation parameters, +.>
The downlink ray parameters are shared for the transverse measuring line;
step B4: multiple frequency domain plane wave domainSecondary wave data
Performing inverse Fourier transform to obtain multiple data of time domain plane wave domain >
Step C: multiple data from time domain plane wave domain
Extracting the common wave-detecting point ray parameter gather, respectively and alternately performing two times of inverse three-dimensional Tau-p conversion and inverse Fourier conversion to obtain multiple wave seismic data M (t, x) of a time-space domain r ,y r ,x s ,y s ) The method specifically comprises the following steps:
step C1: multiple data from time domain plane wave domain
Extracting a common detector point ray parameter trace set +.>
Step C2: first inverse three-dimensional Tau-p transform
For common detection point ray parameter trace set
Obtaining a common-detection-point ray parameter gather m after the first inverse transformation of the frequency domain through the first inverse three-dimensional Tau-p transformation 2 (t,x s ,y s )；
Wherein, the common detection point radial parameter gather
The corresponding frequency domain data is +.>
Co-detection after first inverse transformationWave point ray parameter trace set m 2 (t,x s ,y s ) The corresponding frequency domain data is +.>
The two can be obtained through calculation in a frequency domain, and the calculation formula is as follows:
wherein mf is 2 For the common-detector-point ray parameter gather m after the first inverse transformation of the frequency domain 2 (t,x s ,y s ) Corresponding data
mf 1 Common-detector-point ray parameter gather for frequency domain
Corresponding data +.>
Step C3: first inverse Fourier transform
First inverse transformed common-detector-point ray parameter trace m for frequency domain 2 (t,x s ,y s ) Corresponding data
Performing the first inverse Fourier transform to obtain a common-detection-point ray parameter gather after the first inverse transform
And the plurality of common-wave-point radiation parameter gathers after the first inverse transformation are +.>
Composition of three-dimensional Tau-p Domain dataset->
Step C4: second inverse three-dimensional Tau-p transformation
From three-dimensional Tau-p domain data sets
Co-shot gather of three-dimensional Tau-p domain data is extracted>
Performing a second inverse three-dimensional Tau-p transformation to obtain a multiple three-dimensional common shot point gather m 4 (t,x r ,y r )；
Wherein, the three-dimensional Tau-p domain data is common shot point gather
The corresponding frequency domain data is +.>
Three-dimensional common shot point gather m 4 (t,x r ,y r ) The corresponding frequency domain data is mf 4 (t,x r ,y r ) The two can be obtained through calculation in a frequency domain, and the calculation formula is as follows:
wherein mf is 4 For three-dimensional common shot point gather m 4 (t,x r ,y r ) The corresponding frequency domain data is mf 4 (t,x r ,y r )；
mf 3 co-shot gather for three-dimensional Tau-p domain data in the frequency domain 
Corresponding data
Step C5: second inverse Fourier transform
Data mf corresponding to three-dimensional common shot point gather of frequency domain 4 (t,x r ,y r ) Performing a second inverse Fourier transform to obtain a three-dimensional common shot point gather m in a time domain 4 (t,x r ,y r ) And three-dimensional common shot point gather m of a plurality of time domains 4 (t,x r ,y r ) Multiple arrays M (t, x) forming a time-space domain r ,y r ,x s ,y s )；
Wherein the multiple array M (t, x) r ,y r ,x s ,y s ) As shown in fig. 5.
Step D: time-space domain seismic data D (t, x r ,y r ,x s ,y s ) And multiple array M (t, x) of time-space domain r ,y r ,x s ,y s ) And performing self-adaptive subtraction operation to obtain a result after the multiple suppression of the three-dimensional plane wave domain seismic data. The method specifically comprises the following steps:
assuming that the overall energy of the seismic data after multiple suppression is minimum;
first input time-space domain seismic data D (t, x r ,y r ,x s ,y s )；
Then the time-space domain seismic data D (t, x r ,y r ,x s ,y s ) And multiple array M (t, x) of time-space domain r ,y r ,x s ,y s ) And performing self-adaptive subtraction operation.
Wherein, the multiple wave array M (t, x) r ,y r ,x s ,y s ) And time-space domain seismic data D (t, x r ,y r ,x s ,y s ) The subtracted part is adaptively subtracted from the input data of fig. 1 as shown in fig. 5, resulting in fig. 6. As is apparent from a comparison of fig. 6 and fig. 1, the multiples at the corresponding positions are suppressed by the multiple suppression.
While the invention has been described in detail in the foregoing general description and specific examples, it will be apparent to those skilled in the art that modifications and improvements can be made thereto. Accordingly, such modifications or improvements may be made without departing from the spirit of the invention and are intended to be within the scope of the invention as claimed.
Claims (9)
1. The method for suppressing the multiple waves of the three-dimensional plane wave domain seismic data is characterized by comprising the following steps of:
from time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Selecting a three-dimensional common shot point gather, and obtaining three-dimensional plane wave Tau-p domain seismic data through twice three-dimensional Tau-p transformation
For the three-dimensional plane wave Tau-p domain seismic numberAccording to
Performing multiple suppression processing to obtain multiple data of time domain plane wave domain +.>
Multiple data from the time domain plane wave domain
Extracting the common wave-detecting point ray parameter gather, respectively and alternately performing two times of inverse three-dimensional Tau-p conversion and inverse Fourier conversion to obtain multiple wave seismic data M (t, x) of a time-space domain r ,y r ,x s ,y s )；
Multiple arrays M (t, x) of the time-space domain r ,y r ,x s ,y s ) And the time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Performing self-adaptive subtraction operation to obtain a result after multiple suppression of the three-dimensional plane wave domain seismic data;
wherein t is the longitudinal time, x r Seismic traces in the inline direction, y r Seismic traces in the crossline direction, x s A seismic source in the direction of the inline line, y s Is a seismic source in the cross-line direction,
for the inline direction detector radiation parameters, < +.>
The wave-detecting point radiation parameters are the wave-detecting point radiation parameters in the transverse measuring line direction; />
For the source ray parameter in the inline direction, +. >
The method is characterized in that the method is a transverse survey line direction seismic source ray parameter;
obtaining three-dimensional plane wave Tau-p domain seismic data D through twice three-dimensional Tau-p transformation 2taup 
The method comprises the following specific steps:
step A1: first three-dimensional Tau-p transform
From time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Selecting three-dimensional common shot point gather d 1 (t,x r ,y r ) The method comprises the steps of carrying out a first treatment on the surface of the Gathering d of three-dimensional common shot points 1 (t,x r ,y r ) Performing Fourier transform, first three-dimensional Tau-p transform and inverse Fourier transform to obtain seismic data of the wave detection point Tau-p time domain
Seismic data of the wave detection point Tau-p time domain +.>
Time-space domain seismic data D (t, x) placed at corresponding locations r ,y r ,x s ,y s ) In (1) obtaining the seismic data of the wave detection point Tau-p domain +.>
Step A2: second order three-dimensional Tau-p transformation
From the detector points Tau-p domain seismic data
Selecting a three-dimensional common-detection-point ray parameter gather d 2 (t，x s ，y s ) After Fourier transformation, performing a second three-dimensional Tau-p transformation on the data corresponding to the main frequency in the obtained data, and calculating three-dimensional Tau-p transformation data +_ of the seismic data corresponding to the main frequency>
Formation of three-dimensional Tau-p domain seismic data in full frequency domain>
From a three-dimensional common-detector-point ray parameter gather d 2 (t，x s ，y s ) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data +.>
Firstly, after inverse Fourier transformation, the seismic data of the three-dimensional plane wave domain are put into the corresponding position +. >
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained>
Where f is the frequency.
2. The method according to claim 1, wherein the step A1 comprises the specific steps of:
step A11: given a known time-space domain seismic data D (t, x r ,y r ,x s ,y s )；
Step A12: from time-space domain seismic data D (t, x r ,y r ,x s ,y s ) Selecting a three-dimensional common shot point gather d 1 (t,x r ,y r )；
Step A13: gathering d of three-dimensional common shot points 1 (t,x r ,y r ) Performing Fourier transform, first three-dimensional Tau-p transform and inverse Fourier transform to obtain seismic data of the wave detection point Tau-p time domain
Wherein t is the longitudinal time, < >>
For the inline direction detector radiation parameters, < +.>
The wave-detecting point radiation parameters are the wave-detecting point radiation parameters in the transverse measuring line direction;
step A14: seismic data of a wave-detecting point Tau-p time domain
Time-space domain seismic data D (t, x) placed at corresponding locations r ,y r ,x s ,y s ) In (1) obtaining the seismic data of the wave detection point Tau-p domain +.>
/>
3. The method according to claim 2, wherein said step a13 comprises the specific steps of:
step A131: three-dimensional common shot point gather d 1 (t,x r ,y r ) Obtaining three-dimensional common shot point gather data df of a frequency domain by utilizing Fourier transformation to the frequency domain 1 (f,x r ,y r ) Wherein f is frequency;
step A132: data df of three-dimensional common shot point gather in frequency domain 1 (f,x r ,y r ) Obtaining the frequency of the wave detection point through the first three-dimensional Tau-p conversionTau-p domain data of domain
Step A133: tau-p domain data of wave detector frequency domain
Performing inverse Fourier transform to obtain seismic data of the detection point Tau-p time domain>
4. The method according to claim 1, wherein the step A2 comprises the specific steps of:
step A21: from the wave-detecting point Tau-p domain seismic data
Selecting a three-dimensional common-detection-point ray parameter gather d 2 (t，x s ，y s )；
Step A22: gather d of three-dimensional common-detection-point ray parameters 2 (t，x s ，y s ) Obtaining three-dimensional common shot point ray parameter gather data df of a frequency domain through Fourier transformation to the frequency domain 2 (f，x s ，y s )；
Step A23: three-dimensional common shot point ray parameter gather data df according to frequency domain 2 (f，x s ，y s ) Firstly, carrying out a second three-dimensional Tau-p transformation according to the seismic data corresponding to the main frequency of the seismic data to obtain a shot Tau-p frequency domain gather of the seismic data corresponding to the main frequency
Then, three-dimensional Tau-p transformation data of the corresponding seismic data outside the main frequency is calculated according to the three-dimensional Tau-p transformation data
The shot Tau-p frequency domain gather of which the main frequency corresponds to the seismic data is +.>
Three-dimensional Tau-p transform data corresponding to seismic data outside the dominant frequency>
Combining to form three-dimensional Tau-p domain seismic data in full frequency domain>
Step A24: gather d of three-dimensional common-detection-point ray parameters 2 (t，x s ，y s ) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data in full frequency domain are->
Firstly, after inverse Fourier transformation, three-dimensional plane wave domain seismic data of the corresponding position are put in +.>
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained>
5. The method of claim 4, wherein the step a23 comprises the specific steps of:
step A231: selecting three-dimensional common shot point ray parameter gather data df of frequency domain 2 (f，x s ，y s ) Data corresponding to the main frequency in the three-dimensional common shot gather data df0 of the seismic data corresponding to the main frequency 2 (f，x s ，y s ) And for three-dimensional common shot point gather data df0 of which main frequency corresponds to seismic data 2 (f，x s ，y s ) Performing a second three-dimensional Tau-p transformation to obtain a shot Tau-p frequency domain gather of the seismic data with main frequency
Step A232: shot Tau-p frequency domain gather corresponding to seismic data according to main frequency
Calculating a diagonal constraint matrix W of the inline and a diagonal constraint matrix V of the crossline;
step A233: according to the transverse line diagonal constraint matrix W and the transverse line diagonal constraint matrix V, calculating three-dimensional Tau-p conversion data of the corresponding seismic data except the main frequency
Step a234: shot Tau-p frequency domain gather of main frequency corresponding to seismic data 
Three-dimensional Tau-p transform data corresponding to seismic data outside the dominant frequency>
Combining to form three-dimensional Tau-p domain seismic data in full frequency domain>
6. The method of claim 1 wherein multiple suppression processing is performed to obtain multiple data of the time domain plane wave domain
The method comprises the following specific steps:
step B1: for the three-dimensional plane wave Tau-p domain seismic data
Performing Fourier transform along time direction to obtain seismic data +.>
Step B2: seismic data from the plane wave domain of the frequency domain
Selecting an initial frequency slice from the results of (a)>
For the selected initial frequency slice->
Performing linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slice ∈ ->
Step B3: slicing multiple data frequencies
Multiple data array for composing frequency domain plane wave domain>
Step B4: multiple data array for frequency domain plane wave domain
Performing inverse Fourier transform to obtainMultiple data of inter-domain plane wave domain +.>
7. The method of claim 6, wherein said step B2 comprises the specific steps of:
step B21: seismic data from the plane wave domain of the frequency domain 
Selecting a frequency slice as an initial frequency slice +.>
Slicing the initial frequency>
Performing linear mapping to obtain mapped frequency slices +.>
Step B22: slicing the mapped frequency
Performing squaring operation to obtain frequency slice after multiple data mapping>
Step B23: frequency slicing after mapping the multiple data
Inverse linear mapping to obtain multiple data frequencySlice->
8. The method of claim 1, wherein the time-space domain multiple seismic data M (t, x) is obtained by alternately performing two inverse three-dimensional Tau-p transforms and inverse fourier transforms, respectively r ,y r ,x s ,y s ) The method comprises the following steps:
step C2: for common detection point ray parameter trace set
Obtaining a common-detection-point ray parameter gather m after the first inverse transformation of the frequency domain through the first inverse three-dimensional Tau-p transformation 2 (t,x s ,y s )；
Step C3: first inverse transformed common-detector-point ray parameter trace m for frequency domain 2 (t,x s ,y s ) Corresponding data
Performing the first inverse Fourier transform to obtain a common-detection-point ray parameter gather after the first inverse transform
And the plurality of common-wave-point radiation parameter gathers after the first inverse transformation are +.>
Composition of three-dimensional Tau-p Domain dataset->
Step C4: from three-dimensional Tau-p domain data sets 
One three is extractedCo-shot gather of dimension Tau-p domain data>
Performing a second inverse three-dimensional Tau-p transformation to obtain a multiple three-dimensional common shot point gather m 4 (t,x r ,y r )；
Step C5: data mf corresponding to three-dimensional common shot point gather of frequency domain 4 (t,x r ,y r ) Performing a second inverse Fourier transform to obtain a three-dimensional common shot point gather m in a time domain 4 (t,x r ,y r ) And three-dimensional common shot point gather m of a plurality of time domains 4 (t,x r ,y r ) Multiple arrays M (t, x) forming a time-space domain r ,y r ,x s ,y s )。
9. The method of claim 1, wherein obtaining the result of multiple suppression of the three-dimensional plane wave domain seismic data comprises the specific steps of:
assuming that the overall energy of the seismic data after multiple suppression is minimum;
first input time-space domain seismic data D (t, x r ,y r ,x s ,y s )；
Then the time-space domain seismic data D (t, x r ,y r ,x s ,y s ) And multiple array M (t, x) of time-space domain r ,y r ,x s ,y s ) And performing self-adaptive subtraction operation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111461899.3A CN114185095B (en) | 2021-12-02 | 2021-12-02 | Method for suppressing multiple waves of three-dimensional plane wave domain seismic data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111461899.3A CN114185095B (en) | 2021-12-02 | 2021-12-02 | Method for suppressing multiple waves of three-dimensional plane wave domain seismic data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114185095A CN114185095A (en) | 2022-03-15 |
| CN114185095B true CN114185095B (en) | 2023-05-16 |
Family
ID=80542026
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111461899.3A Active CN114185095B (en) | 2021-12-02 | 2021-12-02 | Method for suppressing multiple waves of three-dimensional plane wave domain seismic data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114185095B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116755151B (en) * | 2023-06-16 | 2024-03-22 | 广东海洋大学 | Anti-diffraction wave field separation method based on iterative focusing optimization |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107884829A (en) * | 2017-12-06 | 2018-04-06 | 东华理工大学 | A kind of method for combining compacting shallow sea OBC Multiple Attenuation in Seismic Data |
| CN108983284A (en) * | 2018-06-22 | 2018-12-11 | 中国石油大学(华东) | A kind of domain f-p ghost reflection drawing method suitable for marine tiltedly cable data |
Family Cites Families (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NO20020984L (en) * | 2001-02-27 | 2002-08-28 | Petroleum Geo Services Inc | Angle dependent surface multiple attenuation for two-component sensor data from the seabed |
| US7477992B2 (en) * | 2005-02-18 | 2009-01-13 | Exxonmobil Upstream Research Company | Method for combining seismic data sets |
| US8082106B2 (en) * | 2007-08-16 | 2011-12-20 | Bp Corporation North America Inc. | 3D surface related multiple elimination for wide azimuth seismic data |
| CN102378922B (en) * | 2009-03-31 | 2015-12-16 | 国际壳牌研究有限公司 | For filtering in TAU-P territory and stablizing low-frequency method in deconvolution |
| EP3191872B1 (en) * | 2014-09-10 | 2023-07-26 | CGG Services SAS | Wave-field reconstruction using a reflection from a variable sea surface |
| CN104345343B (en) * | 2014-12-01 | 2017-02-22 | 中国海洋石油总公司 | Prediction method for complex seabed related interbed multiples |
| CN104914466B (en) * | 2015-06-26 | 2017-12-12 | 中国石油大学(华东) | A kind of method for improving seismic data resolution |
| WO2017218722A1 (en) * | 2016-06-15 | 2017-12-21 | Schlumberger Technology Corporation | Systems and methods for attenuating noise in seismic data and reconstructing wavefields based on the seismic data |
| US20180180755A1 (en) * | 2016-12-28 | 2018-06-28 | China Petroleum & Chemical Corporation | Method for angle-domain common image gather |
| CN109239772B (en) * | 2018-08-31 | 2020-02-14 | 中海石油(中国)有限公司湛江分公司 | Shallow water multi-wave model directivity prediction method for characteristic wave domain |
| WO2020072147A1 (en) * | 2018-10-03 | 2020-04-09 | Mayo Foundation For Medical Education And Research | Phase velocity imaging using an imaging system |
| CN110082823B (en) * | 2019-05-09 | 2020-08-14 | 中国石油大学(北京) | Seismic data interpolation method and device |
| CN110780348B (en) * | 2019-11-01 | 2021-07-09 | 中国石油大学(华东) | Primary wave and multiple combined imaging method and system based on stereo imaging conditions |
| CN111239827B (en) * | 2020-03-09 | 2021-07-30 | 吉林大学 | Three-dimensional seismic data multiple suppression method based on local similarity coefficient |
| CN111505719B (en) * | 2020-05-27 | 2022-05-27 | 中海石油(中国)有限公司湛江分公司 | Diffracted multiple suppression method based on wave field decomposition |
| CN112213775B (en) * | 2020-09-16 | 2023-01-24 | 中国石油天然气股份有限公司 | Fidelity frequency-boosting method for high-coverage-frequency pre-stack seismic data |
| CN112925023B (en) * | 2021-02-01 | 2022-03-22 | 中国石油大学(北京) | Full wave field inversion seismic data multiple suppression method |
| CN113391346B (en) * | 2021-06-03 | 2022-11-22 | 中国海洋石油集团有限公司 | Tau-p domain interlayer multiple prediction method based on inphase axis slope extrapolation edge expansion |
-
2021
- 2021-12-02 CN CN202111461899.3A patent/CN114185095B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107884829A (en) * | 2017-12-06 | 2018-04-06 | 东华理工大学 | A kind of method for combining compacting shallow sea OBC Multiple Attenuation in Seismic Data |
| CN108983284A (en) * | 2018-06-22 | 2018-12-11 | 中国石油大学(华东) | A kind of domain f-p ghost reflection drawing method suitable for marine tiltedly cable data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114185095A (en) | 2022-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4380059A (en) | F-K Filtering of multiple reflections from a seismic section | |
| Toomey et al. | Tomographic imaging of the shallow crustal structure of the East Pacific Rise at 9° 30′ N | |
| EP3339910B1 (en) | Device and method for model-based deblending | |
| CN104280777B (en) | Method for suppressing interference of seismic data multiples on land | |
| EA032186B1 (en) | Seismic adaptive focusing | |
| Van Dalen et al. | Retrieving surface waves from ambient seismic noise using seismic interferometry by multidimensional deconvolution | |
| CN112327362B (en) | Submarine multiple prediction and tracking attenuation method in velocity domain | |
| Abbad et al. | Automatic nonhyperbolic velocity analysis | |
| CN114185095B (en) | Method for suppressing multiple waves of three-dimensional plane wave domain seismic data | |
| An et al. | Auto-pick first breaks with complex raypaths for undulate surface conditions | |
| CN109738953B (en) | Complete multiple suppression method based on wavelet domain frequency division energy compensation | |
| Al-Heety et al. | Random and coherent noise attenuation for 2D land seismic reflection line acquired in Iraq | |
| Zhao et al. | Wavelet-crosscorrelation-based interferometric redatuming in 4D seismic | |
| CN110703332B (en) | Ghost wave compression method | |
| US5442591A (en) | Method for adaptively suppressing noise transients in summed co-sensor seismic recordings | |
| Egan et al. | The influence of spatial sampling on resolution | |
| CA2870892A1 (en) | A computerized method and a computer program product for determining a resulting data set representative of a geological region of interest | |
| Shang et al. | Broadband data with a new low frequency source—Acquisition and processing example from the Gulf-of-Mexico | |
| Santos et al. | High-resolution marine seismic data processing using a single-channel sparker system: Bahia shelf, Brazil | |
| CN112327361B (en) | Inclination interference elimination method based on linear same-phase axis iterative tracking attenuation | |
| Abbasi et al. | Attenuating long-period multiples in short-offset 2D streamer data: Gulf of California | |
| CN110398773B (en) | Recovery and reconstruction method for partial missing seismic data | |
| Liu et al. | Making true amplitude angle domain common image gathers using invertible radon transform | |
| WO2024016572A1 (en) | Noise identification method and apparatus, and device and storage medium | |
| Xu et al. | Reconstruction of the near‐offset gap in marine seismic data using seismic interferometric interpolation |
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 |