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 dataMultiple suppression processing is carried out to obtain multiple data of a time domain plane wave domain
Multiple data from the time domain plane wave domainExtracting 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 dataSeismic data of Tau-p time domain of detection pointTime-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 dataSelecting 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 domainAfter inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domainIn 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 domainWhere 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 pointTime-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 pointPerforming 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 pointsSelecting 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 frequencyThen 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 dataCorresponding the main frequency to the shot Tau-p frequency domain gather of the seismic dataThree-dimensional Tau-p transform data corresponding to seismic data other than primary frequenciesCombining 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 dataThree-dimensional Tau-p domain seismic data under full frequency domainAfter inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domainIn 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 frequencyCalculating 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 dataThree-dimensional Tau-p transform data corresponding to seismic data other than primary frequenciesCombining 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 domainThe method comprises the following steps:
step B1: for the three-dimensional plane wave Tau-p domain seismic dataFourier 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 domainSelecting an initial frequency slice from the results of (1)Slicing the selected initial frequencyPerforming linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slice
Step B3: frequency slicing multiple multiples of dataMultiple data array for forming frequency domain plane wave domain
Step B4: multiple data array for frequency domain plane wave domainPerforming 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 domainSelecting one frequency slice as initial frequency sliceSlicing the initial frequency Linear mapping to obtain mapped frequency slice
Step B22: slicing the mapped frequenciesPerforming squaring operation to obtain frequency slice after multiple wave data mapping
Step B23: mapping the multi-wave data to a frequency sliceObtaining 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 gatherObtaining 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 dataPerforming first inverse Fourier transform to obtain a co-detection point ray parameter gather after first inverse transformAnd collecting a plurality of co-detection point ray parameter traces after first inverse transformationComposing a three-dimensional Tau-p domain dataset
Step C4: from three-dimensional Tau-p domain data setsCommon shot gather for extracting three-dimensional Tau-p domain dataPerforming 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 conversionThe 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 inventionSelecting 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 transformationThe result of (1);
FIG. 4 is a time domain seismic data plot from the demodulator probe Tau-p of FIG. 3After Fourier transform, selecting a frequency slice from the result as an initial frequency sliceAnd slicing the initial frequencyMapped frequency slice obtained by linear mappingWhere FIG. 4(a) is the initial frequency slice before linear mappingFIG. 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 inventionDisplaying 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 transformationThe 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 domainSeismic 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 domainWhere 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 pointThe 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);
As a three-dimensional common shot gather df of the frequency domain1An operator in the longitudinal direction of (1);
as a three-dimensional common shot gather df of the frequency domain1The operator of the transverse measuring line direction;
i is an identity matrix; lambda is a damping operator; h represents conjugate transpose;
as a three-dimensional common shot gather df of the frequency domain1Of the longitudinal line direction operatorThe conjugate transpose of (1);
as a three-dimensional common shot gather df of the frequency domain1Transverse survey direction operatorThe conjugate transpose of (c).
Specifically, three-dimensional common shot gather df in frequency domain1Operator of the longitudinal direction of the lineGather by three-dimensional common shot gather d1(t,xr,yr) Seismic trace x of the inline directionrAnd ray parameters of inline directionAnd (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 directionGather by three-dimensional common shot gather d1(t,xr,yr) Transverse survey line direction of seismic trace yrAnd ray parameters in the transverse directionAnd (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 pointPerforming 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 domainAs shown in fig. 2.
Step A14: seismic data of Tau-p time domain of detection pointTime-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 pointsSelecting 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 domainAfter inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domainIn 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 pointsSelecting 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 domainExtracting 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 dataExtracting 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 frequencyThen 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 dataCorresponding the main frequency to the shot Tau-p frequency domain gather of the seismic dataThree-dimensional Tau-p transform data corresponding to seismic data other than primary frequenciesCombining 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 dataThe calculation expression of (a) is:
in the formula, df02taupShot Tau-p frequency domain gathers for seismic data corresponding to primary frequenciesdf02Three-dimensional common shot gather data df0 for primary frequency corresponding seismic data2(f,xs,ys);
Shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupAn operator in the longitudinal direction of (1);
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;
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupLongitudinal direction ofOperator ofThe conjugate transpose of (1);
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupOperator of transverse line directionThe conjugate transpose of (c).
Specifically, the dominant frequency corresponds to the shot Tau-p frequency domain gather of the seismic dataOperator of the longitudinal direction of the lineFrom the co-detector ray parameter gather d2(t,xs,ys) Of the longitudinal direction of the seismic source xsAnd ray parameters of inline directionCalculating 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 frequencyOperator of transverse line directionFrom the co-detector ray parameter gather d2(t,xs,ys) Transverse survey line direction of seismic source ysAnd ray parameters in the transverse directionCalculating 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 frequencyCalculating 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 lineThe 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;
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupOperator of the longitudinal direction of the lineThe conjugate transpose of (1);
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;
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupOperator of transverse line directionThe conjugate transpose of (c).
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;
matrix W composed of diagonal elements of diagonal constraint matrix W representing vertical linesiiRow i and column i of (a), the number being equal in value to the primary frequency corresponding seismic eventShot Tau-p frequency domain gather df0 of data2taupOperator of the longitudinal direction of the lineCorresponding three-dimensional common shot gather data df0XS;
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;
matrix V composed of diagonal elements of diagonal constraint matrix V representing horizontal linesiiRow i and column i of (a), which is equal in value to the shot Tau-p frequency domain gather df0 of the primary frequency corresponding seismic data2taupOperator L of transverse direction of the lineysCorresponding three-dimensional common shot gather data df0YS;
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 frequenciesThe 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;
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupAn operator in the longitudinal direction of (1);
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; h represents conjugate transpose;
lambda' is a damping operator in the longitudinal line direction; mu' is a damping operator in the transverse measuring direction;
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupOperator of the longitudinal direction of the lineThe conjugate transpose of (1);
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupOperator of transverse line directionThe conjugate transpose of (1);
step A234: corresponding the main frequency to the shot Tau-p frequency domain gather of the seismic dataThree-dimensional Tau-p transform data corresponding to seismic data other than primary frequenciesCombining 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 dataThree-dimensional Tau-p domain seismic data under full frequency domainAfter inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domainIn 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 domainPerforming 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-pThree-dimensional plane wave domain seismic data embedded into corresponding positionsIn 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 transformationThe 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 dataMultiple 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 dataFourier transform is performed along the time direction, and then an initial frequency slice is selected from the resultSlicing the selected initial frequencyPerforming linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slice
Then frequency slicing the multiple multiples dataMultiple data array for forming frequency domain plane wave domain
Finally, the multiple data array of the frequency domain plane wave domainPerforming 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 domainFourier 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 domainIn selecting an initial frequency sliceSlicing the selected initial frequencyPerforming 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 domainSelecting one frequency slice as initial frequency sliceSlicing the initial frequencyLinear mapping to obtain mapped frequency slice
Specifically, willAndrespectively substituting into the initial frequency sliceIn (1)Andto obtain the mapped frequency sliceThen the initial frequency sliceAnd mapped frequency sliceThe 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 sliceAnd mapped frequency sliceEach of which comprises a plurality of small matrices,wherein each small matrix is a block matrix, and the frequency slice after mappingEach 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 parameterEach small matrix occupies a common downlink ray parameter PdOne crossline of the parallel ray parametersA numerical value; mapped frequency sliceEach 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 parameterEach small matrix is a block matrix and occupies the common uplink ray parameter PoOne transverse line of the same uplink ray parameterNumerical values.
Step B22: slicing the mapped frequenciesPerforming squaring operation to obtain frequency slice after multiple wave data mapping
Wherein each mapped frequency slice is slicedPerforming 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 mappingThe 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 dataObtaining multi-wave data frequency slice by inverse linear mapping
Step B3: frequency slicing multiple multiples of dataMultiple data array for forming frequency domain plane wave domain
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 domainPerforming inverse Fourier transform to obtain multiple data of time domain plane wave domain
And C: multiple data from time domain plane wave domainExtracting 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 domainExtracting a common detection point ray parameter gather
Step C2: first inverse three-dimensional Tau-p transform
For co-detection point ray parameter gatherObtaining 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 gatherCorresponding frequency domain data isCommon-detection-point ray parameter gather m after first inverse transformation2(t,xs,ys) Corresponding frequency domain data isThe 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
Shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupAn operator in the longitudinal direction of (1);
shot Tau-p frequency domain gather df0 for seismic data corresponding to primary frequencies2taupThe operator of the transverse measuring line direction;
Step C3: first inverse Fourier transform
Common detection point ray parameter gather m after first inverse transformation of frequency domain2(t,xs,ys) Corresponding dataPerforming first inverse Fourier transform to obtain a co-detection point ray parameter gather after first inverse transformAnd collecting a plurality of co-detection point ray parameter traces after first inverse transformationComposing a three-dimensional Tau-p domain dataset
Step C4: second inverse three-dimensional Tau-p transform
From three-dimensional Tau-p domain data setsCommon shot gather for extracting three-dimensional Tau-p domain dataPerforming 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 dataCorresponding 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);
As a three-dimensional common shot gather df of the frequency domain1An operator in the longitudinal direction of (1);
as a three-dimensional common shot gather df of the frequency domain1The operator of the transverse measuring line direction;
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 dataMultiple suppression processing is carried out to obtain multiple data of a time domain plane wave domain
Multiple data from the time domain plane wave domainExtracting 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 transformationsThe 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 domainSeismic data of Tau-p time domain of detection pointTime-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 dataSelecting 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 domainAfter inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domainIn 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 domainWhere 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;
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
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 pointsSelecting 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 frequencyThen 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 dataCorresponding the main frequency to the shot Tau-p frequency domain gather of the seismic dataThree-dimensional Tau-p transform data corresponding to seismic data other than primary frequenciesCombining 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 dataThree-dimensional Tau-p domain seismic data under full frequency domainAfter inverse Fourier transform, the seismic data are put into corresponding positions to three-dimensional plane wave domainIn 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 frequencyCalculating 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
7. The method according to claim 1, wherein the multiple-suppressing process is performed to acquire multiple data of the time-domain plane wave domainThe method comprises the following specific steps:
step B1: for the three-dimensional plane wave Tau-p domain seismic dataFourier 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 domainSelecting an initial frequency slice from the results of (1)Slicing the selected initial frequencyPerforming linear mapping, squaring operation and inverse linear mapping to obtain multiple data frequency slice
Step B3: frequency slicing multiple multiples of dataMultiple data array for forming frequency domain plane wave domain
Step B4: multiple data array for frequency domain plane wave domainPerforming 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 domainSelecting one frequency slice as initial frequency sliceSlicing the initial frequency Linear mapping to obtain mapped frequency slice
Step B22: slicing the mapped frequenciesPerforming squaring operation to obtain frequency slice after multiple wave data 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 gatherObtaining 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 dataPerforming first inverse Fourier transform to obtain a co-detection point ray parameter gather after first inverse transformAnd collecting a plurality of co-detection point ray parameter traces after first inverse transformationComposing a three-dimensional Tau-p domain dataset
Step C4: from three-dimensional Tau-p domain data setsCommon shot gather for extracting three-dimensional Tau-p domain dataPerforming 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 | |
US8103453B2 (en) | Method of seismic data interpolation by projection on convex sets | |
WO2017106127A1 (en) | Automated near surface analysis by surface-consistent refraction methods | |
US20120075954A1 (en) | Device and Method for Calculating 3D Angle Gathers from Reverse Time Migration | |
EP3029491B1 (en) | Increasing similarity between seismic datasets | |
US9250341B2 (en) | Device and method for extrapolating specular energy of reverse time migration three dimensional angle gathers | |
CN104280777B (en) | Method for suppressing interference of seismic data multiples on land | |
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 | |
WO2006054181A1 (en) | Method for processing at least two sets of seismic data | |
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) | |
CN113805237A (en) | Method and system for offset land cross-spread seismic using compressed sensing models | |
KR101931092B1 (en) | Hybrid method for full waveform inversion using simultaneous and sequential source method | |
Xu et al. | Radon-domain interferometric interpolation for reconstruction of the near-offset gap in marine seismic data | |
CN113687417B (en) | Three-dimensional prestack seismic data interbed multiple prediction and suppression method | |
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 | |
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 |