CN114185095A - 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
- CN114185095A CN114185095A CN202111461899.3A CN202111461899A CN114185095A CN 114185095 A CN114185095 A CN 114185095A CN 202111461899 A CN202111461899 A CN 202111461899A CN 114185095 A CN114185095 A CN 114185095A
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000009466 transformation Effects 0.000 claims abstract description 80
- 238000001514 detection method Methods 0.000 claims abstract description 75
- 230000001629 suppression Effects 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 44
- 238000013507 mapping Methods 0.000 claims description 22
- 230000001131 transforming effect Effects 0.000 claims description 9
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 238000013506 data mapping Methods 0.000 claims description 5
- 238000000844 transformation Methods 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 17
- 150000001875 compounds Chemical class 0.000 description 6
- 238000013016 damping Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 239000000523 sample Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000011144 upstream manufacturing 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 shot-point-sharing gather from the time-space domain seismic data, and performing two times of three-dimensional Tau-p transformation and multiple suppression to obtain multiple data of a time domain plane wave domain; extracting a co-detection wave point ray parameter gather from multi-time wave data of a time domain plane wave domain, and alternately performing two times of inverse three-dimensional Tau-p transformation and inverse Fourier transformation respectively to obtain multi-time wave seismic data of a time-space domain; and performing self-adaptive subtraction operation on the multiple array of the time-space domain and the time-space domain seismic data to obtain a result of the three-dimensional plane wave domain seismic data after multiple suppression. 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 and 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 seismic data of a three-dimensional plane wave domain.
Background
In marine seismic exploration, strong multiple interference exists in seismic data under the influence of a strong reflection coefficient of a seawater surface. The existence of multiple waves complicates the wave field of the seismic waves, causes serious interference to the amplitude and energy of reflected waves, limits the frequency bandwidth of seismic data, blurs the knowledge of underground geological structures, and is a serious obstacle to the processing and interpretation of the seismic data.
Three-dimensional data volume obtained by three-dimensional acquisition processed by a common two-dimensional seismic data processing method has errors, so that the seismic exploration precision is low. With the increasing depth of seismic exploration, the existing two-dimensional exploration cannot meet 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 of:
from time-space domain seismic data D (t, x)r,yr,xs,ys) Three-dimensional shot point-sharing gather is selected, and three-dimensional plane wave Tau-p domain seismic data are obtained through two times of three-dimensional Tau-p transformation
For the three-dimensional plane wave Tau-p domain seismic data
Multiple suppression processing is carried out to obtain multiple data of a time domain plane wave domain
Multiple data from the time domain plane wave domain
Extracting a co-detection point ray parameter gather, and respectively and alternately performing two times of three-dimensional Tau-p transformation and inverse Fourier transformation to obtain multiple seismic data M (t, x) of a time-space domainr,yr,xs,ys)；
The multiple array M (t, x) of the time-space domainr,yr,xs,ys) And said time-space domain seismic data D (t, x)r,yr,xs,ys) Performing self-adaptive subtraction operation to obtain a result of suppressing the multiple waves of the three-dimensional plane wave domain seismic data;
wherein t is the longitudinal time, xrSeismic traces in the inline direction, yrSeismic traces in the transverse direction, xsAs a source of inline, ysIs a seismic source in the transverse survey line direction,
is a ray parameter of a wave detection point in the longitudinal line direction,
detecting point ray parameters in the transverse measuring line direction;
is a seismic source ray parameter in the longitudinal survey line direction,
and the parameters are the seismic source ray parameters in the transverse survey line direction.
Further, the three-dimensional plane wave Tau-p domain seismic data are obtained through two times of 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,yr,xs,ys) Selecting three-dimensional common shot gather d1(t,xr,yr) (ii) a Collecting three-dimensional common shot gathers d1(t,xr,yr) Fourier transformation, first three-dimensional Tau-p transformation and inverse Fourier transformation are carried out to obtain a detection point Tau-p time domain seismic data
Seismic data of Tau-p time domain of detection point
Time-space domain seismic data D (t, x) put into corresponding positionsr,yr,xs,ys) In the method, seismic data of a wave detection point Tau-p domain are obtained
Step A2: second three-dimensional Tau-p transform
From the geophone point Tau-p domain seismic data
Selecting a three-dimensional co-detection point ray parameter gather d2(t，xs，ys) After Fourier transformation, performing secondary three-dimensional Tau-p transformation on data corresponding to main frequency in the obtained data, and calculating three-dimensional Tau-p transformation data corresponding to seismic data except the main frequency
Forming three-dimensional Tau-p domain seismic data under full frequency domain
From three-dimensional co-detector point ray parameter gather d2(t，xs，ys) Transforming to obtain three-dimensional plane wave domain seismic data
Full frequency domainSeismic data of lower three-dimensional Tau-p domain
After inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domain
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 the number ofrIs the seismic trace in the direction of the longitudinal measuring line,
for examining ray parameters, y, of point in the longitudinal direction of the linerIs the seismic trace in the transverse survey line direction,
for examining point ray parameters, x, in the transverse direction of the linesIs a seismic source in the direction of the longitudinal survey line,
for source ray parameters, y, in the inline directionsIs a seismic source in the transverse survey line direction,
and the parameters are the seismic source ray parameters in the transverse survey line direction.
Further, the step a1 includes the following steps:
step A11: given a known time-space domain seismic data D (t, x)r,yr,xs,ys)；
Step A12: from time-space domain seismic data D (t, x)r,yr,xs,ys) Selecting a three-dimensional common shot gather d1(t,xr,yr)；
Step A13: collecting three-dimensional common shot gathers d1(t,xr,yr) Fourier transformation, first three-dimensional Tau-p transformation and inverse Fourier transformation are carried out to obtain seismic data of a wave detection point Tau-p time domain
Where t is the time in the longitudinal direction,
is a ray parameter of a wave detection point in the longitudinal line direction,
detecting point ray parameters in the transverse measuring line direction;
step A14: seismic data of Tau-p time domain of detection point
Time-space domain seismic data D (t, x) put into corresponding positionsr,yr,xs,ys) In the method, seismic data of a wave detection point Tau-p domain are obtained
Specifically, the step a13 includes the following steps:
step A131: three-dimensional shot-sharing gather d1(t,xr,yr) Obtaining three-dimensional common shot gather data df of a frequency domain by utilizing Fourier transform to the frequency domain1(f,xr,yr) Wherein f is frequency;
step A132: three-dimensional common shot gather data df of frequency domain1(f,xr,yr) Obtaining Tau-p domain data of a frequency domain of a detection point through the first three-dimensional Tau-p transformation
Step A133: tau-p domain data of frequency domain of detection point
Performing inverse Fourier transform to obtain seismic data of a detection point Tau-p time domain
Specifically, the step a2 includes the following steps:
step A21: tau-p domain seismic data from geophone points
Selecting a three-dimensional co-detection point ray parameter gather d2(t，xs，ys)；
Step A22: collecting the ray parameter trace d of the three-dimensional common detection point2(t，xs，ys) Fourier transform is carried out to the frequency domain to obtain three-dimensional common shot point ray parameter gather data df of the frequency domain2(f，xs，ys)；
Step A23: three-dimensional common shot point ray parameter gather data df according to frequency domain2(f，xs，ys) Firstly, carrying out the 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 corresponding to the seismic data except the main frequency is calculated according to the three-dimensional Tau-p transformation data
Corresponding the main frequency to the shot Tau-p frequency domain gather of the seismic data
Three-dimensional Tau-p transform data corresponding to seismic data other than primary frequencies
Combining to form three-dimensional Tau-p domain seismic data under full frequency domain
Step A24: collecting the ray parameter trace d of the three-dimensional common detection point2(t，xs，ys) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data under full frequency domain
After inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domain
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained
More specifically, the step a23 includes the following steps:
step A231: selecting three-dimensional common shot point ray parameter gather data df of frequency domain2(f，xs，ys) Data corresponding to the medium main frequency form three-dimensional common shot gather data df0 of seismic data corresponding to the main frequency2(f，xs，ys) And three-dimensional common shot gather data df0 for seismic data corresponding to the dominant frequency2(f，xs，ys) Carrying out the second three-dimensional Tau-p transformation to obtain a shot Tau-p frequency domain gather of the seismic data corresponding to the main frequency
Step A232: shot Tau-p frequency domain gather corresponding to seismic data according to main frequency
Calculating a longitudinal measuring line diagonal constraint matrix W and a transverse measuring line diagonal constraint matrix V;
step A233: according to longitudinal directionCalculating three-dimensional Tau-p transformation data of corresponding seismic data except the main frequency by using a diagonal measuring constraint matrix W and a transverse diagonal measuring constraint matrix V
Step A234: corresponding the main frequency to the shot Tau-p frequency domain gather of the seismic data
Three-dimensional Tau-p transform data corresponding to seismic data other than primary frequencies
Combining to form three-dimensional Tau-p domain seismic data under full frequency domain
Further, multiple suppression processing is carried out 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
Fourier transform is carried out along the time direction to obtain the seismic data of the plane wave domain of the frequency domain
Step B2: seismic data from plane wave domain of frequency domain
Selecting an initial frequency slice from the results of (1)
Slicing the selected initial frequency
Performing linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slice
Step B3: frequency slicing multiple multiples of data
Multiple data array for forming 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 the number ofrIs the seismic trace in the direction of the longitudinal measuring line,
for examining ray parameters, y, of point in the longitudinal direction of the linerIs the seismic trace in the transverse survey line direction,
for examining point ray parameters, x, in the transverse direction of the linesIs a seismic source in the direction of the longitudinal survey line,
for source ray parameters, y, in the inline directionsIs a seismic source in the transverse survey line direction,
and the parameters are the seismic source ray parameters in the transverse survey line direction.
Specifically, the step B2 includes the following steps:
step B21: seismic data from the plane wave domain of the frequency domain
Selecting one frequency slice as initial frequency slice
Slicing the initial frequency
Linear mapping to obtain mapped frequency slice
Step B22: slicing the mapped frequencies
Performing squaring operation to obtain frequency slice after multiple wave data mapping
Step B23: mapping the multi-wave data to a frequency slice
Obtaining multi-wave data frequency slice by inverse linear mapping
Further, performing inverse three-dimensional Tau-p transformation and inverse Fourier transformation twice alternately to obtain multiple seismic data M (t, x) of the time-space domainr,yr,xs,ys) The method comprises the following steps:
step C2: for co-detection point ray parameter gather
Obtaining a common-detection-point ray parameter gather m after the first inverse transformation of a frequency domain through the first inverse three-dimensional Tau-p transformation2(t,xs,ys)；
Step C3: common detection point ray parameter gather m after first inverse transformation of frequency domain2(t,xs,ys) Corresponding data
Performing first inverse Fourier transform to obtain a co-detection point ray parameter gather after first inverse transform
And collecting a plurality of co-detection point ray parameter traces after first inverse transformation
Composing a three-dimensional Tau-p domain dataset
Step C4: from three-dimensional Tau-p domain data sets
Common shot gather for extracting three-dimensional Tau-p domain data
Performing a second inverse three-dimensional Tau-p transformation to obtain a multiple three-dimensional common shot gather m4(t,xr,yr)；
Step C5: data mf corresponding to three-dimensional common shot gather of frequency domain4(t,xr,yr) Performing second Fourier transform to obtain a time-domain three-dimensional common shot gather m4(t,xr,yr) And collecting m three-dimensional common shot point gathers of a plurality of time domains4(t,xr,yr) Multiple array M (t, x) forming time-space domainr,yr,xs,ys)。
Further, the step of obtaining the result of the three-dimensional plane wave domain seismic data after multiple suppression comprises the following steps:
assuming that the overall energy of the seismic data after the multiple suppression is the minimum;
firstly inputting time-space domain seismic data D (t, x)r,yr,xs,ys)；
Then the time-space domain seismic data D (t, x)r,yr,xs,ys) Multiple array M (t, x) of sum-time-space domainr,yr,xs,ys) And performing adaptive subtraction operation.
Compared with the prior art, the invention has the beneficial effects that:
the invention discloses a method for suppressing multiple waves of three-dimensional plane wave domain seismic data, which adopts a three-dimensional seismic method for acquiring underground information in a certain area, can three-dimensionally know the underground geological structure condition from a three-dimensional space, 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 a biplane wave domain on the seismic data based on Tau-p transformation; under the condition that the seismic data transversely have fluctuation, the seismic data can be effectively compressed by transforming the seismic data to a biplane wave domain, the amount of the compressed seismic data is reduced, the calculation amount can also be reduced to a certain extent, the requirement on a computer storage medium is reduced, and the calculation efficiency is improved; the multiple-suppression method does not require information from the subsurface medium and is a fully data-driven method. The method for suppressing the multiple waves of the three-dimensional plane wave domain seismic data can provide a sectional, planar and three-dimensional underground geological structure image, greatly improves the seismic exploration accuracy, is particularly effective for areas with complicated and changeable underground geological 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 diagram of an input three-dimensional common shot gather d provided in embodiment 1 of the present invention1(t,xr,yr) The display result of (1);
FIG. 2 is a diagram of a three-dimensional common shot gather d provided in embodiment 1 of the present invention1(t,xr,yr) Seismic data of demodulator probe Tau-p time domain obtained by first Tau-p conversion
The result of (1);
FIG. 3 is a diagram for providing Tau-p time domain seismic data for a survey point in accordance with embodiment 1 of the present invention
Selecting three-dimensional co-detection point ray parameter gather d according to keywords2(t，xs，ys) And after Fourier transformation, the second three-dimensional Tau-p transformation and inverse Fourier transformation are carried out to obtain the seismic data of the time domain three-dimensional Tau-p transformation
The result of (1);
FIG. 4 is a time domain seismic data plot from the demodulator probe Tau-p of FIG. 3
After Fourier transform, selecting a frequency slice from the result as an initial frequency slice
And slicing the initial frequency
Mapped frequency slice obtained by linear mapping
Where 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 in the time-space domain according to embodiment 1 of the present invention
Displaying the result in a time-space domain;
FIG. 6 is a schematic diagram of time-space domain seismic data D (t, x) provided in embodiment 1 of the present inventionr,yr,xs,ys) Multiple array M (t, x) of sum-time-space domainr,yr,xs,ys) The first-order wave result is obtained by performing the adaptive subtraction operation, wherein the position indicated by the arrow in fig. 6 originally has a plurality of times.
Detailed Description
The invention provides a method for suppressing multiple waves 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 sequentially and respectively adopting a first Tau-p transformation algorithm and a second Tau-p transformation algorithm based on main frequency constraint to transform the three-dimensional seismic data according to the characteristics of the three-dimensional seismic data that detection points are dense and shot points are sparse, so that the representation precision of the plane wave domain on time-space domain seismic data is improved; performing multiple prediction processing on the frequency slice in the plane wave domain by mapping, multiplying and the like; and then, the result of the multiple prediction is inversely transformed back to a time-space domain, and adaptive subtraction is carried out, so that the purposes of multiple suppression and data precision improvement are achieved.
Example 1
The embodiment 1 provides a method for suppressing multiple waves of three-dimensional plane wave domain seismic data, which comprises the following steps:
step A: two three-dimensional Tau-p transformations
From time-space domain seismic data D (t, x)r,yr,xs,ys) Selecting a three-dimensional common shot gather, and obtaining the seismic data of a three-dimensional plane wave domain through two times of three-dimensional Tau-p transformation
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,yr,xs,ys) Selecting three-dimensional common shot gather d1(t,xr,yr) (ii) a Collecting three-dimensional common shot gathers d1(t,xr,yr) Fourier transformation, first three-dimensional Tau-p transformation and inverse Fourier transformation are carried out to obtain seismic data of a wave detection point Tau-p time domain
Seismic data of Tau-p time domain of detection point
Time-space domain seismic data D (t, x) put into corresponding positionsr,yr,xs,ys) In the method, seismic data of a wave detection point Tau-p domain are obtained
Wherein t is the longitudinal time, xrSeismic traces in the inline direction, yrSeismic traces in the transverse direction, xsAs a source of inline, ysAs a seismic source in the transverse direction, pxrFor examining point ray parameters, p, for inline directionyrAnd the ray parameters are the ray parameters of the wave detection point in the transverse measuring line direction.
The detailed steps are as follows:
step A11: given a known time-space domain seismic data D (t, x)r,yr,xs,ys)；
Wherein t is the longitudinal time, xrSeismic traces in the inline direction, yrSeismic traces in the transverse direction, xsIn the longitudinal direction of the lineThe seismic source of (a), ysA seismic source in the transverse survey line direction;
step A12: from time-space domain seismic data D (t, x)r,yr,xs,ys) Selecting a three-dimensional common shot gather d1(t,xr,yr)；
Step A13: collecting three-dimensional common shot gathers d1(t,xr,yr) Fourier transformation, first three-dimensional Tau-p transformation and inverse Fourier transformation are carried out to obtain seismic data of a wave detection point Tau-p time domain
Where t is the time in the longitudinal direction,
is a ray parameter of a wave detection point in the longitudinal line direction,
and the ray parameters are the ray parameters of the wave detection point in the transverse measuring line direction.
Wherein, the first three-dimensional Tau-p transformation process is carried out in a frequency domain, and the more detailed steps are as follows:
step A131: three-dimensional shot-sharing gather d1(t,xr,yr) Obtaining three-dimensional common shot gather data df of the frequency domain through Fourier transformation to the frequency domain1(f,xr,yr) Wherein f is frequency;
step A132: three-dimensional common shot gather data df of frequency domain1(f,xr,yr) Obtaining Tau-p domain data of a frequency domain of a detection point through the first three-dimensional Tau-p transformation
Wherein, Tau-p domain data of frequency domain of wave detection point
The calculation expression of (a) is:
in the formula (I), the compound is shown in the specification,
for Tau-p domain data in the frequency domain of the detection point, df1As a three-dimensional common shot gather df of the frequency domain1(f,xr,yr)；
i is an identity matrix; lambda is a damping operator; h represents conjugate transpose;
Specifically, three-dimensional common shot gather df in frequency domain1Operator of the longitudinal direction of the line
Gather by three-dimensional common shot gather d1(t，xr，yr) Seismic trace x of the inline directionrAnd ray parameters of inline direction
And (4) calculating according to the following formula:
in the formula, xrIs xr1，xr2，…，xrnWherein r is1，r2，…，rnIs composed of rxThe subscript number in the middle not only represents the seismic trace in the longitudinal measuring line direction, but also represents the offset distance in the longitudinal measuring line direction;
Three-dimensional common shot gather df under frequency domain1Operator of transverse line direction
Gather by three-dimensional common shot gather d1(t，xr，yr) Transverse survey line direction of seismic trace yrAnd ray parameters in the transverse direction
And (4) calculating according to the following formula:
in the formula, yrDenotes yr1，yr2，…，yrnWherein r is1，r2，…，rnIs yrThe subscript number in (1) not only indicates the seismic trace in the crossline direction, but also indicates the offset in the crossline direction;
Step A133: tau-p domain data of frequency domain of detection point
Performing inverse Fourier transform to obtain seismic data of a detection point Tau-p time domain
As a specific embodiment, when three-dimensional common shot gather d1(t，xr，yr) When the input record is shown in figure 1, Fourier transform, first three-dimensional Tau-p transform and inverse Fourier transform are carried out to obtain seismic data of a demodulator probe Tau-p time domain
As shown in fig. 2.
Step A14: seismic data of Tau-p time domain of detection point
Time-space domain seismic data D (t, x) put into corresponding positionsr,yr,xs,ys) In the method, seismic data of a wave detection point Tau-p domain are obtained
Step A2: second three-dimensional Tau-p transform
Tau-p domain seismic data from geophone points
Selecting a three-dimensional co-detection point ray parameter gather d2(t，xs，ys) After Fourier transformation, performing secondary three-dimensional Tau-p transformation on data corresponding to main frequency in the obtained data, and calculating three-dimensional Tau-p transformation data corresponding to seismic data except the main frequency
Forming three-dimensional Tau-p domain seismic data under full frequency domain
From three-dimensional co-detector point ray parameter gather d2(t，xs，ys) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data under full frequency domain
After inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domain
In the method, a three-dimensional plane wave Tau-p domain is obtainedSeismic data
Wherein t is longitudinal time and f is frequency; x is the number ofrIs the seismic trace in the direction of the longitudinal measuring line,
for examining ray parameters, y, of point in the longitudinal direction of the linerIs the seismic trace in the transverse survey line direction,
detecting point ray parameters in the transverse measuring line direction; x is the number ofsIs a seismic source in the direction of the longitudinal survey line,
for source ray parameters, y, in the inline directionsIs a seismic source in the transverse survey line direction,
and the parameters are the seismic source ray parameters in the transverse survey line direction.
The detailed steps are as follows:
step A21: tau-p domain seismic data from geophone points
Selecting a three-dimensional co-detection point ray parameter gather d2(t，xs，ys)；
As a specific implementation mode, 1, 4 and 5 are selected as keywords, and seismic data from a detection point Tau-p domain
Extracting corresponding dimension according to keywords in the process of selecting a three-dimensional co-detection point ray parameter gather d2(t，xs，ys) From geophone point Tau-p domain seismic data
Extracting the 1 st, 4 th and 5 th dimensionsThree-dimensional co-detector point ray parameter gather d2(t，xs，ys)。
Step A22: collecting the ray parameter trace d of the three-dimensional common detection point2(t，xs，ys) Fourier transform is carried out to the frequency domain to obtain three-dimensional common shot point ray parameter gather data df of the frequency domain2(f，xs，ys)；
Step A23: three-dimensional common shot point ray parameter gather data df according to frequency domain2(f，xs，ys) Firstly, carrying out the 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 corresponding to the seismic data except the main frequency is calculated according to the three-dimensional Tau-p transformation data
Corresponding the main frequency to the shot Tau-p frequency domain gather of the seismic data
Three-dimensional Tau-p transform data corresponding to seismic data other than primary frequencies
Combining to form three-dimensional Tau-p domain seismic data under 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 domain2(f，xs，ys) Data corresponding to the medium main frequency form three-dimensional common shot gather data df0 of seismic data corresponding to the main frequency2(f，xs，ys) And three-dimensional common shot gather data df0 for seismic data corresponding to the dominant frequency2(f，xs，ys) Carrying out the second three-dimensional Tau-p transformation to obtain a shot Tau-p frequency domain gather of the seismic data corresponding to the main frequency
Wherein the main frequency corresponds to the shot Tau-p frequency domain gather of the seismic data
The calculation expression of (a) is:
in the formula, df02taupShot Tau-p frequency domain gathers for seismic data corresponding to primary frequencies
df02Three-dimensional common shot gather data df0 for primary frequency corresponding seismic data2(f，xs，ys)；
i is an identity matrix; mu is a damping parameter; h represents conjugate transpose;
Specifically, the dominant frequency corresponds to the shot Tau-p frequency domain gather of the seismic data
Operator of the longitudinal direction of the line
From the co-detector ray parameter gather d2(t，xs，ys) Of the longitudinal direction of the seismic source xsAnd ray parameters of inline direction
Calculating according to a calculation formula;
in the formula, xsDenotes xs1，xs2，…，xsnWherein s is1，s2，…，snIs xsSubscript number of (1), xsNot only represents the seismic source of the longitudinal line direction, but also represents the offset distance of the longitudinal line direction;
Shot Tau-p frequency domain gather of seismic data corresponding to main frequency
Operator of transverse line direction
From the co-detector ray parameter gather d2(t，xs，ys) Transverse survey line direction of seismic source ysAnd ray parameters in the transverse direction
Calculating according to a calculation formula;
in the formula, ysDenotes ys1，ys2，…，ysnWherein s is1，s2，…，snIs ysSubscript number of (1), ysNot only represents the seismic source in the transverse direction, but also represents the offset in the transverse direction;
Step A232: shot Tau-p frequency domain gather corresponding to seismic data according to main frequency
Calculating a longitudinal measuring line diagonal constraint matrix W and a transverse measuring line diagonal constraint matrix V;
the detailed calculation process is as follows:
setting the main frequency corresponding to the shot Tau-p frequency domain trace gather df0 of the seismic data2taupOperator of the longitudinal direction of the line
The corresponding three-dimensional common shot gather data is df0XSThe expression is as follows:
in the formula (I), the compound is shown in the specification,
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupAn operator in the longitudinal direction of (1);
i is an identity matrix; mu is a damping parameter; h represents conjugate transpose;
df02three-dimensional common shot gather data df0 for primary frequency corresponding seismic data2(f，xs，ys)；
Setting the main frequency corresponding to the shot Tau-p frequency domain trace gather df0 of the seismic data2taupOperator L of transverse direction of the lineysCorresponding three-dimensionalShot gather data is df0YSThe expression is as follows:
in the formula (I), the compound is shown in the specification,
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupThe operator of the transverse measuring line direction;
i is an identity matrix; mu is a damping parameter; h represents conjugate transpose;
Since the off-diagonal elements of the diagonal constraint matrix W of the longitudinal lines are all 0, the calculation only needs to calculate the matrix W formed by the diagonal elements of the diagonal constraint matrix W of the longitudinal linesiiThe expression is as follows:
wherein epsilon 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; wiiIs a matrix formed by diagonal elements of a diagonal constraint matrix W of a longitudinal line;
Since the off-diagonal elements of the diagonal constraint matrix V of the crosslines are all 0, the calculation only needs to calculate the matrix V formed by the diagonal elements of the diagonal constraint matrix V of the crosslinesiiThe expression is as follows:
where ζ is the stability factor across the line; i is the number of rows or columns of diagonal positions of the diagonal-diagonal constraint matrix V of the crossline; viiIs a matrix formed by diagonal elements of a diagonal constraint matrix V of a horizontal line;
Step A233: according to the longitudinal measuring line diagonal constraint matrix W and the transverse measuring line diagonal constraint matrix V, three-dimensional Tau-p transformation data corresponding to the seismic data except the main frequency are calculated
Wherein the three-dimensional Tau-p transform data corresponds to seismic data other than the primary frequencies
The calculation formula of (A) is as follows:
in the formula, df12taupThree-dimensional Tau-p transform data for seismic data corresponding to frequencies other than the dominant frequency
df02Three-dimensional common shot gather data df0 for primary frequency corresponding seismic data2(f，xs，ys)；
W is a diagonal constraint matrix of the longitudinal lines; v is a diagonal constraint matrix of the transverse measuring line;
i is an identity matrix; h represents conjugate transpose;
lambda' is a damping operator in the longitudinal line direction; mu' is a damping operator in the transverse measuring direction;
step A234: corresponding the main frequency to the shot Tau-p frequency domain gather of the seismic data
Three-dimensional Tau-p transform data corresponding to seismic data other than primary frequencies
Combining to form three-dimensional Tau-p domain seismic data under full frequency domain
Step A24: collecting the ray parameter trace d of the three-dimensional common detection point2(t，xs，ys) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data under full frequency domain
After inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domain
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained
The detailed steps are as follows:
step A241: collecting the ray parameter trace d of the three-dimensional common detection point2(t，xs，ys) Transforming to obtain three-dimensional plane wave domain seismic data
Step A242: three-dimensional Tau-p domain seismic data under full frequency domain
Performing inverse Fourier transform to obtain seismic data of time domain three-dimensional Tau-p transform
Step A243: seismic data transformed with time domain three-dimensional Tau-p
Three-dimensional plane wave domain seismic data embedded into corresponding positions
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained
Specifically, when the three-dimensional co-detection point ray parameter gather d2(t，xs，ys) Obtaining seismic data of time domain three-dimensional Tau-p transformation through Fourier transformation, second three-dimensional Tau-p transformation and inverse Fourier transformation
The results are shown in FIG. 3. It can be seen that significant focal points occur within a limited region of the data, such as between the source ray parameters to 600-.
And B: for three-dimensional plane wave Tau-p domain seismic data
Multiple suppression processing is carried out to obtain multiple data of a time domain plane wave domain
The method specifically comprises the following steps:
firstly aligning three-dimensional plane wave Tau-p domain seismic data
Fourier transform is performed along the time direction, and then an initial frequency slice is selected from the result
Slicing the selected initial frequency
Performing linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slice
Then frequency slicing the multiple multiples data
Multiple data array for forming frequency domain plane wave domain
Finally, the multiple data array of the 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 the number ofrIs the seismic trace in the direction of the longitudinal measuring line,
for examining ray parameters, y, of point in the longitudinal direction of the linerIs the seismic trace in the transverse survey line direction,
for examining point ray parameters, x, in the transverse direction of the linesIs a seismic source in the direction of the longitudinal survey line,
for source ray parameters, y, in the inline directionsIs a seismic source in the transverse survey line direction,
and the parameters are the seismic source ray parameters in the transverse survey line direction.
The more detailed procedure is as follows:
step B1: seismic data for three-dimensional plane wave domain
Fourier transform is carried out on the array along the time direction to obtain the seismic data of the plane wave domain of the frequency domain
Step B2: seismic data from plane wave domain of frequency domain
In selecting an initial frequency slice
Slicing the selected initial frequency
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 plane wave domain of frequency domain
Selecting one frequency slice as initial frequency slice
Slicing the initial frequency
Linear mapping to obtain mapped frequency slice
Specifically, will
And
respectively substituting into the initial frequency slice
In (1)
And
to obtain the mapped frequency slice
Then the initial frequency slice
And mapped frequency slice
The calculation expression of (a) is:
in the formula (I), the compound is shown in the specification,
the parameters of the vertical survey line and the uplink rays are shared,
the parameters of the horizontal line and the upward ray are determined,
is a parameter of a downlink ray shared by the longitudinal survey line,
the parameters are the horizontal survey line and the downlink ray parameters;
is a ray parameter of a wave detection point in the longitudinal line direction,
detecting point ray parameters in the transverse measuring line direction;
is a seismic source ray parameter in the longitudinal survey line direction,
and the parameters are the seismic source ray parameters in the transverse survey line direction.
As can be seen from FIG. 4, the initial frequency slice
And mapped frequency slice
Each of which comprises a plurality of small matrices,wherein each small matrix is a block matrix, and the frequency slice after mapping
Each row of the small matrix is a common downlink ray parameter PdGather, each row of the small matrix is a longitudinal survey line common downlink ray parameter
Each small matrix occupies a common downlink ray parameter PdOne crossline of the parallel ray parameters
A numerical value; mapped frequency slice
Each column of the small matrix is a common upstream ray parameter PoTrace gather, each row of small matrix is longitudinal measuring line common up ray parameter
Each small matrix is a block matrix and occupies the common uplink ray parameter PoOne transverse line of the same uplink ray parameter
Numerical values.
Step B22: slicing the mapped frequencies
Performing squaring operation to obtain frequency slice after multiple wave data mapping
Wherein each mapped frequency slice is sliced
Performing multiplication operation, and performing multiplication operation representation of small matrixPerforming convolution operation on the time-space domain of the uplink ray parameter gather and the downlink ray parameter gather to obtain frequency slices after multiple data mapping
The calculation formula is as follows:
in the formula (I), the compound is shown in the specification,
representing each mapped frequency slice;
representing the frequency slice after the multiple data mapping;
the parameters of the vertical survey line and the uplink rays are shared,
the parameters of the horizontal line and the upward ray are determined,
is a parameter of a downlink ray shared by the longitudinal survey line,
the parameters are the horizontal survey line and the downlink ray parameters;
step B23: mapping a post-frequency slice to multiple data
Obtaining multi-wave data frequency slice by inverse linear mapping
Step B3: frequency slicing multiple multiples of data
Multiple data array for forming frequency domain plane wave domain
Wherein, the multiple data frequency slice
The calculation expression of (a) is:
in the formula (I), the compound is shown in the specification,
representing a frequency slice of the multiple data,
is a ray parameter of a wave detection point in the longitudinal line direction,
detecting point ray parameters in the transverse measuring line direction;
is a seismic source ray parameter in the longitudinal survey line direction,
the seismic source ray parameters are in the transverse survey line direction;
the parameters of the vertical survey line and the uplink rays are shared,
the parameters of the horizontal line and the upward ray are determined,
is a parameter of a downlink ray shared by the longitudinal survey line,
the parameters are the horizontal survey line and the downlink ray parameters;
step B4: multiple data for frequency domain plane wave domain
Performing inverse Fourier transform to obtain multiple data of time domain plane wave domain
And C: multiple data from time domain plane wave domain
Extracting a co-detection point ray parameter gather, and respectively and alternately performing two times of three-dimensional Tau-p transformation and inverse Fourier transformation to obtain multiple seismic data M (t, x) of a time-space domainr,yr,xs,ys) The method specifically comprises the following steps:
step C1: multiple data from time domain plane wave domain
Extracting a common detection point ray parameter gather
Step C2: first inverse three-dimensional Tau-p transform
For co-detection point ray parameter gather
Obtaining a common-detection-point ray parameter gather m after the first inverse transformation of a frequency domain through the first inverse three-dimensional Tau-p transformation2(t,xs,ys)；
Wherein, the common detection point ray parameter gather
Corresponding frequency domain data is
Common-detection-point ray parameter gather m after first inverse transformation2(t,xs,ys) Corresponding frequency domain data is
The two can be obtained by calculation in the frequency domain, and the calculation formula is as follows:
in the formula, mf2Is a common detection point ray parameter gather m after the first inverse transformation of the frequency domain2(t,xs,ys) Corresponding data
mf1co-detector point ray parameter gathers for frequency domain
Corresponding data
Step C3: first inverse Fourier transform
Common detection point ray parameter gather m after first inverse transformation of frequency domain2(t,xs,ys) Corresponding data
Performing first inverse Fourier transform to obtain a co-detection point ray parameter gather after first inverse transform
And collecting a plurality of co-detection point ray parameter traces after first inverse transformation
Composing a three-dimensional Tau-p domain dataset
Step C4: second inverse three-dimensional Tau-p transform
From three-dimensional Tau-p domain data sets
Common shot gather for extracting three-dimensional Tau-p domain data
Performing a second inverse three-dimensional Tau-p transformation to obtain a multiple three-dimensional common shot gather m4(t,xr,yr)；
Wherein, the common shot gather of the three-dimensional Tau-p domain data
Corresponding frequency domain data is
Three-dimensional common shot gather m4(t,xr,yr) Corresponding frequency domain data is mf4(t,xr,yr) The two can be obtained by calculation in the frequency domain, and the calculation formula is as follows:
in the formula, mf4For three-dimensional common shot gather m4(t,xr,yr) Corresponding frequency domain data is mf4(t,xr,yr)；
mf3co-shot gathers for three-dimensional Tau-p domain data in frequency domain
Corresponding data
Step C5: second inverse Fourier transform
Data mf corresponding to three-dimensional common shot gather of frequency domain4(t,xr,yr) Performing second Fourier transform to obtain a time-domain three-dimensional common shot gather m4(t,xr,yr) And collecting m three-dimensional common shot point gathers of a plurality of time domains4(t,xr,yr) Multiple array M (t, x) forming time-space domainr,yr,xs,ys)；
Wherein whenMultiple array M (t, x) in space domainr,yr,xs,ys) As shown in fig. 5.
Step D: time-space domain seismic data D (t, x)r,yr,xs,ys) Multiple array M (t, x) of sum-time-space domainr,yr,xs,ys) And performing self-adaptive subtraction operation to obtain a result of suppressing the multiple waves 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 the multiple suppression is the minimum;
firstly inputting time-space domain seismic data D (t, x)r,yr,xs,ys)；
Then the time-space domain seismic data D (t, x)r,yr,xs,ys) Multiple array M (t, x) of sum-time-space domainr,yr,xs,ys) And performing adaptive subtraction operation.
Wherein, the multiple array M (t, x) of the time-space domainr,yr,xs,ys) And time-space domain seismic data D (t, x)r,yr,xs,ys) The subtracted portion is adaptively subtracted from the input data in fig. 1 as shown in fig. 5, resulting in fig. 6. As is apparent from comparison of fig. 6 with fig. 1, the multiples at the corresponding positions are suppressed by the multiple suppression.
Although the invention has been described in detail above with reference to a general description and specific examples, it will be apparent to one skilled in the art that modifications or improvements may be made thereto based on the invention. Accordingly, such modifications and improvements are intended to be within the scope of the invention as claimed.
Claims (10)
1. A method for suppressing multiple waves of three-dimensional plane wave domain seismic data is characterized by comprising the following steps:
from time-space domain seismic data D (t, x)r,yr,xs,ys) Selecting three-dimensional common shot gather, and repeating the steps of twice and three timesObtaining three-dimensional plane wave Tau-p domain seismic data by dimensional Tau-p transformation
For the three-dimensional plane wave Tau-p domain seismic data
Multiple suppression processing is carried out to obtain multiple data of a time domain plane wave domain
Multiple data from the time domain plane wave domain
Extracting a co-detection point ray parameter gather, and respectively and alternately performing two times of three-dimensional Tau-p transformation and inverse Fourier transformation to obtain multiple seismic data M (t, x) of a time-space domainr,yr,xs,ys)；
The multiple array M (t, x) of the time-space domainr,yr,xs,ys) And said time-space domain seismic data D (t, x)r,yr,xs,ys) Performing self-adaptive subtraction operation to obtain a result of suppressing the multiple waves of the three-dimensional plane wave domain seismic data;
wherein t is the longitudinal time, xrSeismic traces in the inline direction, yrSeismic traces in the transverse direction, xsAs a source of inline, ysIs a seismic source in the transverse survey line direction,
is a ray parameter of a wave detection point in the longitudinal line direction,
detecting point ray parameters in the transverse measuring line direction;
is a seismic source ray parameter in the longitudinal survey line direction,
and the parameters are the seismic source ray parameters in the transverse survey line direction.
2. The method of claim 1, wherein the three-dimensional plane wave Tau-p domain seismic data is acquired through two three-dimensional Tau-p transformations
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,yr,xs,ys) Selecting three-dimensional common shot gather d1(t,xr,yr) (ii) a Collecting three-dimensional common shot gathers d1(t,xr,yr) Fourier transformation, first three-dimensional Tau-p transformation and inverse Fourier transformation are carried out to obtain seismic data of a wave detection point Tau-p time domain
Seismic data of Tau-p time domain of detection point
Time-space domain seismic data D (t, x) put into corresponding positionsr,yr,xs,ys) In the method, seismic data of a wave detection point Tau-p domain are obtained
Step A2: second three-dimensional Tau-p transform
From the geophone point Tau-p domain seismic data
Selecting a three-dimensional co-detection point ray parameter gather d2(t，xs，ys) After Fourier transformation, performing secondary three-dimensional Tau-p transformation on data corresponding to main frequency in the obtained data, and calculating three-dimensional Tau-p transformation data corresponding to seismic data except the main frequency
Forming three-dimensional Tau-p domain seismic data under full frequency domain
From three-dimensional co-detector point ray parameter gather d2(t，xs，ys) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data under full frequency domain
After inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domain
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 the number ofrIs the seismic trace in the direction of the longitudinal measuring line,
for examining ray parameters, y, of point in the longitudinal direction of the linerIs the seismic trace in the transverse survey line direction,
for examining point ray parameters, x, in the transverse direction of the linesIs a seismic source in the direction of the longitudinal survey line,
for source ray parameters, y, in the inline directionsIs a seismic source in the transverse survey line direction,
and the parameters are the seismic source ray parameters in the transverse survey line direction.
3. The method as claimed in claim 2, wherein said step a1 includes the following specific steps:
step A11: given a known time-space domain seismic data D (t, x)r,yr,xs,ys)；
Step A12: from time-space domain seismic data D (t, x)r,yr,xs,ys) Selecting a three-dimensional common shot gather d1(t,xr,yr)；
Step A13: collecting three-dimensional common shot gathers d1(t,xr,yr) Fourier transformation, first three-dimensional Tau-p transformation and inverse Fourier transformation are carried out to obtain seismic data of a wave detection point Tau-p time domain
Where t is the time in the longitudinal direction,
is a ray parameter of a wave detection point in the longitudinal line direction,
detecting point ray parameters in the transverse measuring line direction;
step A14: seismic data of Tau-p time domain of detection point
Time-space domain seismic data D (t, x) put into corresponding positionsr,yr,xs,ys) In the method, seismic data of a wave detection point Tau-p domain are obtained
4. The method as claimed in claim 3, wherein the step A13 includes the following specific steps:
step A131: three-dimensional shot-sharing gather d1(t,xr,yr) Obtaining three-dimensional common shot gather data df of a frequency domain by utilizing Fourier transform to the frequency domain1(f,xr,yr) Wherein f is frequency;
step A132: three-dimensional common shot gather data df of frequency domain1(f,xr,yr) Obtaining Tau-p domain data of a frequency domain of a detection point through the first three-dimensional Tau-p transformation
Step A133: tau-p domain data of frequency domain of detection point
Performing inverse Fourier transform to obtain seismic data of a detection point Tau-p time domain
5. The method as claimed in claim 2, wherein said step a2 includes the following specific steps:
step A21: tau-p domain seismic data from geophone points
Selecting a three-dimensional co-detection point ray parameter gather d2(t，xs，ys)；
Step A22: collecting the ray parameter trace d of the three-dimensional common detection point2(t，xs，ys) Fourier transform is carried out to the frequency domain to obtain three-dimensional common shot point ray parameter gather data df of the frequency domain2(f，xs，ys)；
Step A23: three-dimensional common shot point ray parameter gather data df according to frequency domain2(f，xs，ys) Firstly, carrying out the 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 corresponding to the seismic data except the main frequency is calculated according to the three-dimensional Tau-p transformation data
Corresponding the main frequency to the shot Tau-p frequency domain gather of the seismic data
Three-dimensional Tau-p transform data corresponding to seismic data other than primary frequencies
Combining to form three-dimensional Tau-p domain seismic data under full frequency domain
Step A24: collecting the ray parameter trace d of the three-dimensional common detection point2(t，xs，ys) Transforming to obtain three-dimensional plane wave domain seismic data
Three-dimensional Tau-p domain seismic data under full frequency domain
After inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domain
In the method, three-dimensional plane wave Tau-p domain seismic data are obtained
6. The method as claimed in claim 5, wherein the step A23 includes the following specific steps:
step A231: selecting three-dimensional common shot point ray parameter gather data df of frequency domain2(f，xs，ys) Data corresponding to the medium main frequency form three-dimensional common shot gather data df0 of seismic data corresponding to the main frequency2(f，xs，ys) And three-dimensional common shot gather data df0 for seismic data corresponding to the dominant frequency2(f，xs，ys) Carrying out the second three-dimensional Tau-p transformation to obtain a shot Tau-p frequency domain gather of the seismic data corresponding to the main frequency
Step A232: shot Tau-p frequency domain gather corresponding to seismic data according to main frequency
Calculating a longitudinal measuring line diagonal constraint matrix W and a transverse measuring line diagonal constraint matrix V;
step A233: according to the longitudinal measuring line diagonal constraint matrix W and the transverse measuring line diagonal constraint matrix V, the meterThree-dimensional Tau-p transform data for seismic data other than computed dominant frequencies
Step A234: corresponding the main frequency to the shot Tau-p frequency domain gather of the seismic data
Three-dimensional Tau-p transform data corresponding to seismic data other than primary frequencies
Combining to form three-dimensional Tau-p domain seismic data under full frequency domain
7. The method according to claim 1, wherein the multiple-suppressing process is performed to acquire 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
Fourier transform is carried out along the time direction to obtain the seismic data of the plane wave domain of the frequency domain
Step B2: seismic data from plane wave domain of frequency domain
Selecting an initial frequency slice from the results of (1)
Slicing the selected initial frequency
Performing linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slice
Step B3: frequency slicing multiple multiples of data
Multiple data array for forming 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 the number ofrIs the seismic trace in the direction of the longitudinal measuring line,
for examining ray parameters, y, of point in the longitudinal direction of the linerIs the seismic trace in the transverse survey line direction,
for examining point ray parameters, x, in the transverse direction of the linesIs a seismic source in the direction of the longitudinal survey line,
for source ray parameters, y, in the inline directionsIs a seismic source in the transverse survey line direction,
and the parameters are the seismic source ray parameters in the transverse survey line direction.
8. The method as claimed in claim 3, wherein the step B2 includes the following specific steps:
step B21: seismic data from the plane wave domain of the frequency domain
Selecting one frequency slice as initial frequency slice
Slicing the initial frequency
Linear mapping to obtain mapped frequency slice
Step B22: slicing the mapped frequencies
Performing squaring operation to obtain frequency slice after multiple wave data mapping
Step B23: mapping the multi-wave data to a frequency slice
Obtaining multi-wave data frequency slice by inverse linear mapping
9. The method of claim 1, wherein the inverse three-dimensional Tau-p transform and the inverse fourier transform are alternately performed twice to obtain the multiple seismic data M (t, x) of the time-space domainr,yr,xs,ys) The method comprises the following steps:
step C2: for co-detection point ray parameter gather
Obtaining a common-detection-point ray parameter gather m after the first inverse transformation of a frequency domain through the first inverse three-dimensional Tau-p transformation2(t,xs,ys)；
Step C3: common detection point ray parameter gather m after first inverse transformation of frequency domain2(t,xs,ys) Corresponding data
Performing first inverse Fourier transform to obtain a co-detection point ray parameter gather after first inverse transform
And collecting a plurality of co-detection point ray parameter traces after first inverse transformation
Composing a three-dimensional Tau-p domain dataset
Step C4: from three-dimensional Tau-p domain data sets
Common shot gather for extracting three-dimensional Tau-p domain data
Performing a second inverse three-dimensional Tau-p transformation to obtain a multiple three-dimensional common shot gather m4(t,xr,yr)；
Step C5: data mf corresponding to three-dimensional common shot gather of frequency domain4(t,xr,yr) Performing second Fourier transform to obtain a time-domain three-dimensional common shot gather m4(t,xr,yr) And collecting m three-dimensional common shot point gathers of a plurality of time domains4(t,xr,yr) Multiple array M (t, x) forming time-space domainr,yr,xs,ys)。
10. The method of claim 1, wherein obtaining the multi-suppressed results of the three-dimensional plane wave domain seismic data comprises the specific steps of:
assuming that the overall energy of the seismic data after the multiple suppression is the minimum;
firstly inputting time-space domain seismic data D (t, x)r,yr,xs,ys)；
Then the time-space domain seismic data D (t, x)r,yr,xs,ys) Multiple array M (t, x) of sum-time-space domainr,yr,xs,ys) And performing 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 true CN114185095A (en) | 2022-03-15 |
| CN114185095B 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) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116755151A (en) * | 2023-06-16 | 2023-09-15 | 广东海洋大学 | Anti-diffraction wave field separation method based on iterative focusing optimization |
Citations (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0204625D0 (en) * | 2001-02-27 | 2002-04-10 | Petroleum Geo Services Us Inc | Angle dependent surface multiple attenuation for two-component marine bottom sensor data |
| US20060190181A1 (en) * | 2005-02-18 | 2006-08-24 | Exxonmobil Upstream Research Company | Method for combining seismic data sets |
| GB0814824D0 (en) * | 2007-08-16 | 2008-09-17 | Bp Corp North America Inc | 3D surface related multiple elimination for wide azimuth seismic data |
| WO2010120505A2 (en) * | 2009-03-31 | 2010-10-21 | Shell Oil Company | Method for stabilization of low frequencies in tau-p domain filtering and deconvolution |
| CN104345343A (en) * | 2014-12-01 | 2015-02-11 | 中国海洋石油总公司 | Prediction method for complex seabed related interbed multiples |
| CN104914466A (en) * | 2015-06-26 | 2015-09-16 | 中国石油大学(华东) | Method for improving seismic data resolution |
| WO2016038458A2 (en) * | 2014-09-10 | 2016-03-17 | Cgg Services Sa | Wave-field reconstruction using a reflection from a variable sea surface |
| 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 |
| CN107884829A (en) * | 2017-12-06 | 2018-04-06 | 东华理工大学 | A kind of method for combining compacting shallow sea OBC Multiple Attenuation in Seismic Data |
| US20180180755A1 (en) * | 2016-12-28 | 2018-06-28 | China Petroleum & Chemical Corporation | Method for angle-domain common image gather |
| CN108983284A (en) * | 2018-06-22 | 2018-12-11 | 中国石油大学(华东) | A kind of domain f-p ghost reflection drawing method suitable for marine tiltedly cable data |
| CN109239772A (en) * | 2018-08-31 | 2019-01-18 | 中海石油(中国)有限公司湛江分公司 | A kind of multiple wave pattern directive property prediction technique of the shallow water of feature wave zone |
| CN110082823A (en) * | 2019-05-09 | 2019-08-02 | 中国石油大学(北京) | A kind of seismic data interpolation method and device |
| CN110780348A (en) * | 2019-11-01 | 2020-02-11 | 中国石油大学(华东) | Primary wave and multiple combined imaging method and system based on stereo imaging conditions |
| WO2020072147A1 (en) * | 2018-10-03 | 2020-04-09 | Mayo Foundation For Medical Education And Research | Phase velocity imaging using an imaging system |
| CN111239827A (en) * | 2020-03-09 | 2020-06-05 | 吉林大学 | Three-dimensional seismic data multiple suppression method based on local similarity coefficient |
| CN111505719A (en) * | 2020-05-27 | 2020-08-07 | 中海石油(中国)有限公司湛江分公司 | Diffracted multiple suppression method based on wave field decomposition |
| CN112213775A (en) * | 2020-09-16 | 2021-01-12 | 中国石油天然气股份有限公司 | Fidelity frequency-boosting method for high-coverage-frequency pre-stack seismic data |
| CN112925023A (en) * | 2021-02-01 | 2021-06-08 | 中国石油大学(北京) | Full wave field inversion seismic data multiple suppression method |
| CN113391346A (en) * | 2021-06-03 | 2021-09-14 | 中国海洋石油集团有限公司 | 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 (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0204625D0 (en) * | 2001-02-27 | 2002-04-10 | Petroleum Geo Services Us Inc | Angle dependent surface multiple attenuation for two-component marine bottom sensor data |
| US20060190181A1 (en) * | 2005-02-18 | 2006-08-24 | Exxonmobil Upstream Research Company | Method for combining seismic data sets |
| GB0814824D0 (en) * | 2007-08-16 | 2008-09-17 | Bp Corp North America Inc | 3D surface related multiple elimination for wide azimuth seismic data |
| WO2010120505A2 (en) * | 2009-03-31 | 2010-10-21 | Shell Oil Company | Method for stabilization of low frequencies in tau-p domain filtering and deconvolution |
| WO2016038458A2 (en) * | 2014-09-10 | 2016-03-17 | Cgg Services Sa | Wave-field reconstruction using a reflection from a variable sea surface |
| CN104345343A (en) * | 2014-12-01 | 2015-02-11 | 中国海洋石油总公司 | Prediction method for complex seabed related interbed multiples |
| CN104914466A (en) * | 2015-06-26 | 2015-09-16 | 中国石油大学(华东) | 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 |
| 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 |
| CN109239772A (en) * | 2018-08-31 | 2019-01-18 | 中海石油(中国)有限公司湛江分公司 | A kind of multiple wave pattern directive property prediction technique of the shallow water of feature wave zone |
| WO2020072147A1 (en) * | 2018-10-03 | 2020-04-09 | Mayo Foundation For Medical Education And Research | Phase velocity imaging using an imaging system |
| CN110082823A (en) * | 2019-05-09 | 2019-08-02 | 中国石油大学(北京) | A kind of seismic data interpolation method and device |
| CN110780348A (en) * | 2019-11-01 | 2020-02-11 | 中国石油大学(华东) | Primary wave and multiple combined imaging method and system based on stereo imaging conditions |
| CN111239827A (en) * | 2020-03-09 | 2020-06-05 | 吉林大学 | Three-dimensional seismic data multiple suppression method based on local similarity coefficient |
| CN111505719A (en) * | 2020-05-27 | 2020-08-07 | 中海石油(中国)有限公司湛江分公司 | Diffracted multiple suppression method based on wave field decomposition |
| CN112213775A (en) * | 2020-09-16 | 2021-01-12 | 中国石油天然气股份有限公司 | Fidelity frequency-boosting method for high-coverage-frequency pre-stack seismic data |
| CN112925023A (en) * | 2021-02-01 | 2021-06-08 | 中国石油大学(北京) | Full wave field inversion seismic data multiple suppression method |
| CN113391346A (en) * | 2021-06-03 | 2021-09-14 | 中国海洋石油集团有限公司 | Tau-p domain interlayer multiple prediction method based on inphase axis slope extrapolation edge expansion |
Non-Patent Citations (3)
| Title |
|---|
| 于海铖;徐辉;张学涛;张军华;郑旭刚;: "基于频率慢度域二阶算子的多次波去噪方法研究及应用" * |
| 张红梅;: "线性Radon域预测反褶积在近海地震去除多次波中的应用" * |
| 马继涛;陈小宏;薛亚茹;: "三维表面多次波压制方法" * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116755151A (en) * | 2023-06-16 | 2023-09-15 | 广东海洋大学 | Anti-diffraction wave field separation method based on iterative focusing optimization |
| CN116755151B (en) * | 2023-06-16 | 2024-03-22 | 广东海洋大学 | Anti-diffraction wave field separation method based on iterative focusing optimization |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114185095B (en) | 2023-05-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2356492B1 (en) | Method for separating independent simultaneous sources | |
| US8619498B2 (en) | Device and method for calculating 3D angle gathers from reverse time migration | |
| US8687463B2 (en) | Method for processing at least two sets of seismic data | |
| WO2017106127A1 (en) | Automated near surface analysis by surface-consistent refraction methods | |
| EP3029491B1 (en) | Increasing similarity between seismic datasets | |
| US9250341B2 (en) | Device and method for extrapolating specular energy of reverse time migration three dimensional angle gathers | |
| US7987054B2 (en) | Efficient multiple prediction in two and three dimensions | |
| Schneider Jr et al. | Wave-equation velocity replacement of the low-velocity layer for overthrust-belt data | |
| US6278950B1 (en) | Turning-wave amplitude inversion | |
| GB2460175A (en) | Interpolating seismic data by projection onto convex sets | |
| EP3259619A1 (en) | Method and system of processing seismic data by providing surface aperture common image gathers | |
| Meyerholtz et al. | Convolutional quelling in seismic tomography | |
| Van Dalen et al. | Retrieving surface waves from ambient seismic noise using seismic interferometry by multidimensional deconvolution | |
| EA032186B1 (en) | Seismic adaptive focusing | |
| TSELENTIS et al. | A spectral approach to Moho depths estimation from gravity measurements in Epirus (NW Greece) | |
| KR101931092B1 (en) | Hybrid method for full waveform inversion using simultaneous and sequential source method | |
| CN113805237A (en) | Method and system for offset land cross-spread seismic using compressed sensing models | |
| Xu et al. | Radon-domain interferometric interpolation for reconstruction of the near-offset gap in marine seismic data | |
| Zelt et al. | Prestack depth migration of dense wide-angle seismic data | |
| CN114185095B (en) | Method for suppressing multiple waves of three-dimensional plane wave domain seismic data | |
| Al-Heety et al. | Random and coherent noise attenuation for 2D land seismic reflection line acquired in Iraq | |
| CN113687417B (en) | Three-dimensional prestack seismic data interbed multiple prediction and suppression method | |
| Egan et al. | The influence of spatial sampling on resolution | |
| Alshuhail et al. | Robust estimation of vertical symmetry axis models via joint migration inversion: including multiples in anisotropic parameter estimation | |
| EP2841965A1 (en) | A computerized method and a computer program product for determining a resulting data set representative of a geological region of interest |
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 |