CN106842321B - Reconstruction of seismic data method and apparatus - Google Patents
Reconstruction of seismic data method and apparatus Download PDFInfo
- Publication number
- CN106842321B CN106842321B CN201510885713.5A CN201510885713A CN106842321B CN 106842321 B CN106842321 B CN 106842321B CN 201510885713 A CN201510885713 A CN 201510885713A CN 106842321 B CN106842321 B CN 106842321B
- Authority
- CN
- China
- Prior art keywords
- data
- seismic data
- reconstruct
- threshold
- frequency component
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000004364 calculation method Methods 0.000 claims abstract description 59
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 25
- 230000009466 transformation Effects 0.000 claims description 36
- 230000008569 process Effects 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 13
- 230000000694 effects Effects 0.000 abstract description 9
- 230000001788 irregular Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000002513 implantation Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000013501 data transformation Methods 0.000 description 1
- 238000010230 functional analysis Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000001755 vocal effect Effects 0.000 description 1
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/32—Transforming one recording into another or one representation into another
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 kind of Reconstruction of seismic data method and apparatus, belong to field of seismic exploration.The described method includes: acquisition seismic data;Obtain default seismic data threshold value;The seismic data is reconstructed by convex set projection algorithm in wave atom domain according to the default seismic data threshold value, the seismic data after obtaining reconstruct.The present invention is by being reconstructed the seismic data irregularly lacked with convex set projection algorithm in wave atom domain, it solves when original earthquake data being reconstructed by warp wavelet in the related technology, since the anisotropy of warp wavelet keeps its calculation amount larger, the lower problem of the reconstruction efficiency of seismic data.Under the premise of having reached the precision for guaranteeing the seismic data of reconstruct, the effect of higher computational efficiency is possessed when reconstructing seismic data.
Description
Technical field
The present invention relates to field of seismic exploration, in particular to a kind of Reconstruction of seismic data method and apparatus.
Background technique
Seismic prospecting is the important technology of exploration oil and gas resource before probing, and the technology is in coalfield and engineering geology
Reconnoitre, Study on Regional Geology and the earth's crust research etc., be also used widely.And due to the complexity of environment, the mistake of instrument
The limitation such as difference and economic cost, collected seismic data is often lack sampling, and the seismic data of lack sampling will lead to letter
Breath is irregular and lacks.Irregular and missing original earthquake data, which can generate data processing, to be seriously affected.Therefore, it is necessary to right
Original earthquake data is rebuild.
In the related technology when being rebuild to original earthquake data, usually using the seismic data weight based on sparse transformation
Construction method in the method rebuilds original earthquake data by Qu Bo (curvelet) transformation in sparse transformation,
Warp wavelet can obtain preferable reconstruction effect due to its good directionality, locality and anisotropy.
In the implementation of the present invention, the inventor finds that the existing technology has at least the following problems: the above method is logical
When crossing warp wavelet original earthquake data being reconstructed, since the anisotropy of warp wavelet keeps its calculation amount larger, earthquake
The reconstruction efficiency of data is lower.
Summary of the invention
When in order to solve that original earthquake data is reconstructed by warp wavelet in the prior art, due to warp wavelet
Anisotropy keeps its calculation amount larger, and the lower problem of the reconstruction efficiency of seismic data, the present invention provides a kind of seismic datas
Method for reconstructing and device.The technical solution is as follows:
According to the first aspect of the invention, a kind of Reconstruction of seismic data method is provided, which comprises
Acquire seismic data;
Obtain default seismic data threshold value;
The seismic data is carried out by convex set projection algorithm in wave atom domain according to the default seismic data threshold value
Reconstruct, the seismic data after obtaining reconstruct.
Optionally, the acquisition seismic data, comprising:
The original earthquake data in acquisition time domain;
One-dimensional Fourier transform is done along time orientation to the original earthquake data of the time-domain, obtains the described of frequency domain
Seismic data.
Optionally, the seismic data includes multiple frequency components, described to be existed according to the default seismic data threshold value
The seismic data is reconstructed by convex set projection algorithm for wave atom domain, the seismic data after obtaining reconstruct, comprising:
The first frequency component in the seismic data is preset by the convex set projection algorithm in wave atom domain
The iterative calculation of number, the first frequency component are any frequency component in the seismic data;
Detection iterates to calculate the signal-to-noise ratio of the reconstruct data obtained every time, until signal-to-noise ratio occur is greater than default signal-to-noise ratio threshold
When the reconstruct data of value, stop iterative calculation;
The signal-to-noise ratio is greater than the reconstruct data of default snr threshold as the reconstruct of the first frequency component
Shake data.
Optionally, described that the first frequency in the seismic data is divided by the convex set projection algorithm in wave atom domain
Amount carries out the iterative calculation of preset times, comprising:
The initial value for obtaining kth time iterative calculation, when the k is equal to 1, the initial value is in the seismic data
The first frequency component, when the k is greater than 1, the initial value is the reconstruct data that kth -1 time iterative calculation obtains;
Wave Atom Transformation is carried out to the initial value, obtains wave atom numeric field data;
Threshold process is carried out to the wave atom numeric field data by the default seismic data threshold value;
To treated, wave atom numeric field data carries out wave atom inverse transformation, obtains processing data;
The initial value of the first frequency component is implanted into the processing data by convex set projection formula and obtains the kth
The reconstruct data of secondary iterative calculation,
The convex set projection formula are as follows:
The Dk(f, x, y) is the reconstruct data of the kth time iterative calculation, the D (f, x, y) is the seismic data
In the first frequency component, the S (x, y) be identify given data and missing data sample operator, the WAx,yFor
Positive wave Atom Transformation operator, the Tk(f, x, y) is threshold operator, the Dk-1The reconstruct data iterated to calculate for kth -1 time,
In the k=1, D0=D (f, x, y).
Optionally, described that threshold process, packet are carried out to the wave atom numeric field data by the default seismic data threshold value
It includes:
The wave atom numeric field data is handled according to threshold formula,
The threshold formula are as follows:
The PkFor the threshold value of the kth time iterative calculation, the PkIncluded in the default seismic data threshold value P.
Optionally, described that the signal-to-noise ratio is greater than the reconstruct data of default snr threshold as the first frequency point
After the reconstruct seismic data of amount, the method also includes:
One-dimensional Fourier inversion is done along frequency direction to the reconstruct seismic data of the first frequency component.
According to the second aspect of the invention, a kind of Reconstruction of seismic data device is provided, described device includes:
Data acquisition module is configured as acquisition seismic data;
Threshold value obtains module, is configured as obtaining default seismic data threshold value;
Data reconstruction module is configured as according to the default seismic data threshold value in wave atom domain to the seismic data
It is reconstructed by convex set projection algorithm, the seismic data after obtaining reconstruct.
Optionally, the data acquisition module, is configured as:
The original earthquake data in acquisition time domain;
One-dimensional Fourier transform is done along time orientation to the original earthquake data of the time-domain, obtains the described of frequency domain
Seismic data.
Optionally, the seismic data includes multiple frequency components,
The data reconstruction module, comprising:
Submodule is iterated to calculate, is configured as in wave atom domain through the convex set projection algorithm in the seismic data
First frequency component carry out preset times iterative calculation, the first frequency component be the seismic data in any frequency
Rate component;
Signal-to-noise ratio computation submodule is configured as detecting the signal-to-noise ratio for iterating to calculate the reconstruct data obtained every time, until
When there are the reconstruct data that signal-to-noise ratio is greater than default snr threshold, stop iterative calculation;
Reconstruct data determine submodule, and the reconstruct data for being configured as the signal-to-noise ratio being greater than default snr threshold are made
For the reconstruct seismic data of the first frequency component.
Optionally, the iterative calculation submodule, comprising:
Initial value acquiring unit is configured as obtaining the initial value of kth time iterative calculation, described when the k is equal to 1
Initial value is the first frequency component in the seismic data, and when the k is greater than 1, the initial value changes for kth -1 time
The reconstruct data that generation is calculated;
Wave Atom Transformation unit is configured as carrying out wave Atom Transformation to the initial value, obtains wave atom numeric field data;
Threshold process unit is configured as carrying out threshold to the wave atom numeric field data by the default seismic data threshold value
Value processing;
Wave atom inverse transformation unit is configured as carrying out wave atom inverse transformation to treated wave atom numeric field data, obtain
Handle data;
Data are implanted into unit, are configured as that the initial value of the first frequency component is implanted into institute by convex set projection formula
It states processing data and obtains the reconstruct data of the kth time iterative calculation,
The convex set projection formula are as follows:
The Dk(f, x, y) is the reconstruct data of the kth time iterative calculation, the D (f, x, y) is the seismic data
In the first frequency component, the S (x, y) be identify given data and missing data sample operator, the WAx,yFor
Positive wave Atom Transformation operator, the Tk(f, x, y) is threshold operator, the Dk-1The reconstruct data iterated to calculate for kth -1 time,
In the k=1, D0=D (f, x, y).
Optionally, the threshold process unit, is configured as:
The wave atom numeric field data is handled according to threshold formula,
The threshold formula are as follows:
The PkFor the threshold value of the kth time iterative calculation, the PkIncluded in the default seismic data threshold value P.
Optionally, described device further include:
Fourier inversion module is configured as doing the reconstruct seismic data of the first frequency component along frequency direction
One-dimensional Fourier inversion.
Technical solution provided by the invention has the benefit that
By the way that the seismic data irregularly lacked is reconstructed with convex set projection algorithm in wave atom domain, solves phase
When original earthquake data being reconstructed by warp wavelet in the technology of pass, since the anisotropy of warp wavelet makes its calculation amount
It is larger, the lower problem of the reconstruction efficiency of seismic data.Under the premise of having reached the precision for guaranteeing the seismic data of reconstruct, in weight
Possess the effect of higher computational efficiency when structure seismic data.
Detailed description of the invention
To describe the technical solutions in the embodiments of the present invention more clearly, make required in being described below to embodiment
Attached drawing is briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for
For those of ordinary skill in the art, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is a kind of flow chart of Reconstruction of seismic data method provided in an embodiment of the present invention;
Fig. 2-1 is the flow chart of another Reconstruction of seismic data method provided in an embodiment of the present invention;
Fig. 2-2 is the flow chart iterated to calculate in Fig. 2-1 illustrated embodiment;
Fig. 3-1 is a kind of block diagram of Reconstruction of seismic data device provided in an embodiment of the present invention;
Fig. 3-2 is the block diagram of data reconstruction module in Fig. 3-1 illustrated embodiment;
Fig. 3-3 is the block diagram that submodule is iterated to calculate in Fig. 3-1 illustrated embodiment;
Fig. 3-4 is the block diagram of another Reconstruction of seismic data device provided in an embodiment of the present invention.
Through the above attached drawings, it has been shown that the specific embodiment of the disclosure will be hereinafter described in more detail.These attached drawings
It is not intended to limit the scope of this disclosure concept by any means with verbal description, but is by referring to specific embodiments
Those skilled in the art illustrate the concept of the disclosure.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with attached drawing to embodiment party of the present invention
Formula is described in further detail.
Fig. 1 is a kind of flow chart of Reconstruction of seismic data method provided in an embodiment of the present invention, under this method may include
The several steps in face:
In a step 101, seismic data is acquired.
In a step 102, default seismic data threshold value is obtained.
In step 103, convex set projection algorithm is passed through to seismic data in wave atom domain according to default seismic data threshold value
It is reconstructed, the seismic data after obtaining reconstruct.
In conclusion Reconstruction of seismic data method provided in an embodiment of the present invention, by wave atom domain with convex set projection
The seismic data irregularly lacked is reconstructed in algorithm, solves in the related technology by warp wavelet to original earthquake number
When according to being reconstructed, since the anisotropy of warp wavelet keeps its calculation amount larger, the reconstruction efficiency of seismic data is lower to be asked
Topic.Under the premise of having reached the precision for guaranteeing the seismic data of reconstruct, possess higher computational efficiency when reconstructing seismic data
Effect.
Fig. 2-1 is the flow chart of another Reconstruction of seismic data method provided in an embodiment of the present invention, and this method can wrap
Include following several steps:
In step 201, the original earthquake data in acquisition time domain.
It, can acquisition time domain be primitively first when using Reconstruction of seismic data method provided in an embodiment of the present invention
Shake data.The method of the original earthquake data in acquisition time domain can refer to the relevant technologies, illustratively, the time-domain of acquisition
Original earthquake data can be the three-dimensional irregular missing seismic data of prestack in a certain work area, and the record length of interception can be 1
Second, the sampling interval can be 1 millisecond, totally 1000 sampled points.
In step 202, one-dimensional Fourier transform is done along time orientation to the original earthquake data of time-domain, obtains frequency
The seismic data in domain.
After acquiring the original earthquake data of time-domain, the original earthquake data of time-domain can be done along time orientation
One-dimensional Fourier transform obtains the seismic data of frequency domain.After original earthquake data is transformed into frequency domain, due to only needing to calculate
The data of positive frequency, thus can reduce in subsequent processing by 50% calculation amount.
In step 203, default seismic data threshold value is obtained.
After obtaining the seismic data of frequency domain, available default seismic data threshold value, the default seismic data threshold value
It can be configured in advance, specific set-up mode can refer to the relevant technologies, and this embodiment is not repeated.
In step 204, the first frequency component in seismic data is carried out by convex set projection algorithm in wave atom domain
The iterative calculation of preset times.
After obtaining the seismic data of frequency domain, convex set projection algorithm can be passed through to seismic data in wave atom domain
In first frequency component carry out preset times iterative calculation.Wherein, first frequency component is any frequency in seismic data
Rate component.
As shown in Fig. 2-2, this step may include following several sub-steps:
In sub-step 2041, the initial value of kth time iterative calculation is obtained.
Changing for preset times is carried out to the first frequency component in seismic data by convex set projection algorithm in wave atom domain
When generation calculates, the initial value of kth available first time iterative calculation, when k is equal to 1, initial value is to record in seismic data
First frequency component, when k is greater than 1, initial value is to iterate to calculate obtained reconstruct data kth -1 time.
In sub-step 2042, wave Atom Transformation is carried out to initial value, obtains wave atom numeric field data.
When obtaining the initial value of kth time iterative calculation, wave Atom Transformation can be carried out to initial value, obtain wave atom
Numeric field data.Wave Atom Transformation is the Functional Analysis method with multi-resolution analysis, locality and directionality, compares warp wavelet, wave
Atom Transformation more can sparsely indicate seismic data, better Approaching earthquake lineups.Illustratively, wave original is carried out to initial value
The mathematic(al) representation of son transformation can be with are as follows:
Dk(f, x, y)=WAx,yDk-1。
Wherein, Dk(f, x, y) is the wave atom numeric field data obtained when kth time iterative calculation, WAx,yFor for space variable
X, the positive wave Atom Transformation operator of y, Dk-1For the initial value of kth subwave Atom Transformation, in k=1, D0=D (f, x, y), i.e.,
The first frequency component that record in shake data.
In sub-step 2043, threshold process is carried out to wave atom numeric field data by default seismic data threshold value.
After obtaining wave atom numeric field data, wave atom numeric field data can be carried out at threshold value by default seismic data threshold value
Reason,.Optionally, wave atom numeric field data can be handled according to threshold formula.
Threshold formula are as follows:
Wherein, Tk(f, x, y) is threshold operator, PkFor the threshold value of kth time iterative calculation, PkIncluded in default seismic data
In threshold value P.
In sub-step 2044, to treated, wave atom numeric field data carries out wave atom inverse transformation, obtains processing data.
After carrying out threshold process to wave atom numeric field data, it is anti-wave atom can be carried out to treated wave atom numeric field data
Transformation obtains processing data.Treated wave atom numeric field data and the wave atom contravariant for space variable x, y can be converted
Son is multiplied to obtain processing data, and illustratively, wave atom inverse transformation operator can be WAx,y -1。
In sub-step 2045, the initial value implantation processing data of first frequency component are obtained by convex set projection formula
The reconstruct data of kth time iterative calculation.
It, can will be at the initial value implantation of first frequency component by convex set projection formula after obtaining processing data
Reason data obtain the reconstruct data of kth time iterative calculation.And convex set projection formula can be with are as follows:
Wherein, Dk(f, x, y) is the reconstruct data of kth time iterative calculation, D (f, x, y) be in seismic data first frequently
Rate component, S (x, y) are sample operator, the WA for identifying given data and missing datax,yFor positive wave Atom Transformation operator, Tk(f,
X, y) it is threshold operator, Dk-1The reconstruct data iterated to calculate for kth -1 time, in k=1, D0=D (f, x, y).
In step 205, detection iterates to calculate the signal-to-noise ratio of the reconstruct data obtained every time, is greater than until there is signal-to-noise ratio
When the reconstruct data of default snr threshold, stop iterative calculation.
When being iterated calculating, the signal-to-noise ratio for iterating to calculate the reconstruct data obtained every time can detecte, until occurring
When signal-to-noise ratio is greater than the reconstruct data of default snr threshold, stop iterative calculation.
Illustratively, signal-to-noise ratio computation formula can be with are as follows:
Wherein, SNR is signal-to-noise ratio, Dobs(f, x, y) is to observe data, the i.e. irregular missing seismic data of frequency domain,
Drec(f, x, y) is the reconstruction data observed after data implantation, Dcal(f, x, y) is the reconstruction data observed before data implantation.
It should be noted that Reconstruction of seismic data method provided in an embodiment of the present invention, compared to weight in the related technology
The seismic data built, signal-to-noise ratio can promote 1 to 2 times.
In step 206, signal-to-noise ratio is greater than the reconstruct data of default snr threshold as the weight of first frequency component
Structure seismic data.
When signal-to-noise ratio occur and being greater than the reconstruct data for presetting snr threshold, signal-to-noise ratio can be greater than to default signal-to-noise ratio
Reconstruct seismic data of the reconstruct data of threshold value as first frequency component.So far, with having obtained the reconstruct of first frequency component
Shake data.
In step 207, one-dimensional Fourier inversion is done along frequency direction to the reconstruct seismic data of first frequency component.
It, can reconstruct seismic data edge to first frequency component after obtaining the reconstruct earthquake number of first frequency component
Frequency direction does one-dimensional Fourier inversion, converts it to time-domain, subsequent to seismic data progress correlation analysis to facilitate.
It should be noted that can be by step 204 to step 207, to the earthquake number of each frequency component of frequency domain
According to being rebuild, the reconstruct seismic data of time-domain is obtained.In addition it is also possible to first pass through step 204 to step 206 obtain it is each
The reconstruct seismic data of frequency domain, then the seismic data of the reconstruct of the frequency domains of all frequency components is done along frequency direction one-dimensional
Fourier inversion, the embodiment of the present invention do not restrict.
It should be added that Reconstruction of seismic data method provided in an embodiment of the present invention, by being lacked to irregular
Before the seismic data of mistake carries out wave Atom Transformation, frequency domain is first converted it to, has reached reduction calculation amount, promotes seismic data
Rebuild the effect of efficiency.
In conclusion Reconstruction of seismic data method provided in an embodiment of the present invention, by wave atom domain with convex set projection
The seismic data irregularly lacked is reconstructed in algorithm, solves in the related technology by warp wavelet to original earthquake number
When according to being reconstructed, since the anisotropy of warp wavelet keeps its calculation amount larger, the reconstruction efficiency of seismic data is lower to be asked
Topic.Under the premise of having reached the precision for guaranteeing the seismic data of reconstruct, possess higher computational efficiency when reconstructing seismic data
Effect.
Following is apparatus of the present invention embodiment, can be used for executing embodiment of the present invention method.For apparatus of the present invention reality
Undisclosed details in example is applied, embodiment of the present invention method is please referred to.
Fig. 3-1 is a kind of block diagram of Reconstruction of seismic data device provided in an embodiment of the present invention, Reconstruction of seismic data dress
It sets and includes:
Data acquisition module 310 is configured as acquisition seismic data.
Threshold value obtains module 320, is configured as obtaining default seismic data threshold value.
Data reconstruction module 330 is configured as passing through seismic data in wave atom domain according to default seismic data threshold value
Convex set projection algorithm is reconstructed, the seismic data after obtaining reconstruct.
In conclusion Reconstruction of seismic data device provided in an embodiment of the present invention, by wave atom domain with convex set projection
The seismic data irregularly lacked is reconstructed in algorithm, solves in the related technology by warp wavelet to original earthquake number
When according to being reconstructed, since the anisotropy of warp wavelet keeps its calculation amount larger, the reconstruction efficiency of seismic data is lower to be asked
Topic.Under the premise of having reached the precision for guaranteeing the seismic data of reconstruct, possess higher computational efficiency when reconstructing seismic data
Effect.
Optionally, data acquisition module 310 are configured as:
The original earthquake data in acquisition time domain;
One-dimensional Fourier transform is done along time orientation to the original earthquake data of time-domain, obtains the earthquake number of frequency domain
According to.
Optionally, seismic data includes multiple frequency components, as shown in figure 3-2, data reconstruction module 330, comprising:
Submodule 331 is iterated to calculate, is configured as in wave atom domain through convex set projection algorithm to the in seismic data
One frequency component carries out the iterative calculation of preset times, and first frequency component is any frequency component in seismic data;
Signal-to-noise ratio computation submodule 332 is configured as detecting the signal-to-noise ratio for iterating to calculate the reconstruct data obtained every time, directly
To when there are reconstruct data of the signal-to-noise ratio greater than default snr threshold, stop iterative calculation;
Reconstruct data determine submodule 333, and the reconstruct data for being configured as signal-to-noise ratio being greater than default snr threshold are made
For the reconstruct seismic data of first frequency component.
Optionally, as shown in Fig. 3-3, submodule 331 is iterated to calculate, comprising:
Initial value acquiring unit 331a is configured as obtaining the initial value of kth time iterative calculation, when k is equal to 1, initially
Value is the first frequency component in seismic data, and when k is greater than 1, initial value is the reconstruct number that kth -1 time iterative calculation obtains
According to;
Wave Atom Transformation unit 331b is configured as carrying out wave Atom Transformation to initial value, obtains wave atom numeric field data;
Threshold process unit 331c is configured as carrying out at threshold value wave atom numeric field data by default seismic data threshold value
Reason;
Wave atom inverse transformation unit 331d is configured as carrying out wave atom inverse transformation to treated wave atom numeric field data,
Obtain processing data;
Data are implanted into unit 331e, and being configured as will be at the initial value implantation of first frequency component by convex set projection formula
Reason data obtain the reconstruct data of kth time iterative calculation, convex set projection formula are as follows:
Dk(f, x, y) is the reconstruct data of kth time iterative calculation, D (f, x, y) be that first frequency in seismic data divides
Amount, S (x, y) are sample operator, the WA for identifying given data and missing datax,yFor positive wave Atom Transformation operator, Tk(f,x,y)
For threshold operator, Dk-1The reconstruct data iterated to calculate for kth -1 time, in k=1, D0=D (f, x, y).
Optionally, threshold process unit 331c, is configured as:
Wave atom numeric field data is handled according to threshold formula,
Threshold formula are as follows:
PkFor the threshold value of kth time iterative calculation, PkIncluded in default seismic data threshold value.
Optionally, as shown in Figure 3-4, the device further include:
Fourier inversion module 340 is configured as doing the reconstruct seismic data of first frequency component along frequency direction
One-dimensional Fourier inversion.
It should be added that Reconstruction of seismic data method provided in an embodiment of the present invention, by being lacked to irregular
Before the seismic data of mistake carries out wave Atom Transformation, frequency domain is first converted it to, has reached reduction calculation amount, promotes seismic data
Rebuild the effect of efficiency.
In conclusion Reconstruction of seismic data device provided in an embodiment of the present invention, by wave atom domain with convex set projection
The seismic data irregularly lacked is reconstructed in algorithm, solves in the related technology by warp wavelet to original earthquake number
When according to being reconstructed, since the anisotropy of warp wavelet keeps its calculation amount larger, the reconstruction efficiency of seismic data is lower to be asked
Topic.Under the premise of having reached the precision for guaranteeing the seismic data of reconstruct, possess higher computational efficiency when reconstructing seismic data
Effect.
Those of ordinary skill in the art will appreciate that realizing that all or part of the steps of above-described embodiment can pass through hardware
It completes, relevant hardware can also be instructed to complete by program, the program can store in a kind of computer-readable
In storage medium, storage medium mentioned above can be read-only memory, disk or CD etc..
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all in spirit of the invention and
Within principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (8)
1. a kind of Reconstruction of seismic data method, which is characterized in that the described method includes:
Acquire seismic data;
Obtain default seismic data threshold value;
The seismic data is reconstructed by convex set projection algorithm in wave atom domain according to the default seismic data threshold value,
Seismic data after obtaining reconstruct;
The seismic data includes multiple frequency components, it is described according to the default seismic data threshold value in wave atom domain to institute
It states seismic data to be reconstructed by convex set projection algorithm, the seismic data after obtaining reconstruct, comprising:
Preset times are carried out to the first frequency component in the seismic data by the convex set projection algorithm in wave atom domain
Iterative calculation, the first frequency component be the seismic data in any frequency component;
Detection iterates to calculate the signal-to-noise ratio of the reconstruct data obtained every time, until signal-to-noise ratio occur is greater than default snr threshold
When reconstructing data, stop iterative calculation;
The signal-to-noise ratio is greater than the reconstruct data of default snr threshold as the reconstruct earthquake number of the first frequency component
According to;
It is described that the first frequency component in the seismic data is preset by the convex set projection algorithm in wave atom domain
The iterative calculation of number, comprising:
The initial value for obtaining kth time iterative calculation, when the k is equal to 1, the initial value is described in the seismic data
First frequency component, when the k is greater than 1, the initial value is the reconstruct data that kth -1 time iterative calculation obtains;
Wave Atom Transformation is carried out to the initial value, obtains wave atom numeric field data;
Threshold process is carried out to the wave atom numeric field data by the default seismic data threshold value;
To treated, wave atom numeric field data carries out wave atom inverse transformation, obtains processing data;
The initial value of the first frequency component is implanted into the processing data by convex set projection formula and obtains the kth time repeatedly
The reconstruct data that generation calculates,
The convex set projection formula are as follows:
The Dk(f, x, y) is the reconstruct data of the kth time iterative calculation, the D (f, x, y) is in the seismic data
The first frequency component, the S (x, y) are sample operator, the WA for identifying given data and missing datax,yFor positive wave
Atom Transformation operator, the Tk(f, x, y) is threshold operator, the Dk-1The reconstruct data iterated to calculate for kth -1 time, in institute
When stating k=1, D0=D (f, x, y).
2. the method according to claim 1, wherein the acquisition seismic data, comprising:
The original earthquake data in acquisition time domain;
One-dimensional Fourier transform is done along time orientation to the original earthquake data of the time-domain, obtains the earthquake of frequency domain
Data.
3. the method according to claim 1, wherein it is described by the default seismic data threshold value to the wave
Atom numeric field data carries out threshold process, comprising:
The wave atom numeric field data is handled according to threshold formula,
The threshold formula are as follows:
The PkFor the threshold value of the kth time iterative calculation, the PkIncluded in the default seismic data threshold value P.
4. the method according to claim 1, wherein described be greater than default snr threshold for the signal-to-noise ratio
After reconstructing reconstruct seismic data of the data as the first frequency component, the method also includes:
One-dimensional Fourier inversion is done along frequency direction to the reconstruct seismic data of the first frequency component.
5. a kind of Reconstruction of seismic data device, which is characterized in that described device includes:
Data acquisition module is configured as acquisition seismic data;
Threshold value obtains module, is configured as obtaining default seismic data threshold value;
Data reconstruction module is configured as passing through the seismic data in wave atom domain according to the default seismic data threshold value
Convex set projection algorithm is reconstructed, the seismic data after obtaining reconstruct;
The seismic data includes multiple frequency components,
The data reconstruction module, comprising:
Submodule is iterated to calculate, is configured as in wave atom domain through the convex set projection algorithm to the in the seismic data
One frequency component carries out the iterative calculation of preset times, and the first frequency component is any frequency point in the seismic data
Amount;
Signal-to-noise ratio computation submodule is configured as detecting the signal-to-noise ratio for iterating to calculate the reconstruct data obtained every time, until occurring
When signal-to-noise ratio is greater than the reconstruct data of default snr threshold, stop iterative calculation;
Reconstruct data determine submodule, are configured as the signal-to-noise ratio being greater than the reconstruct data of default snr threshold as institute
State the reconstruct seismic data of first frequency component;
The iterative calculation submodule, comprising:
Initial value acquiring unit is configured as obtaining the initial value of kth time iterative calculation, described initial when the k is equal to 1
Value is the first frequency component in the seismic data, and when the k is greater than 1, the initial value is -1 iteration meter of kth
Obtained reconstruct data;
Wave Atom Transformation unit is configured as carrying out wave Atom Transformation to the initial value, obtains wave atom numeric field data;
Threshold process unit is configured as carrying out at threshold value the wave atom numeric field data by the default seismic data threshold value
Reason;
Wave atom inverse transformation unit is configured as carrying out wave atom inverse transformation to treated wave atom numeric field data, be handled
Data;
Data are implanted into unit, are configured as that the initial value of the first frequency component is implanted into the place by convex set projection formula
Reason data obtain the reconstruct data of the kth time iterative calculation,
The convex set projection formula are as follows:
The Dk(f, x, y) is the reconstruct data of the kth time iterative calculation, the D (f, x, y) is in the seismic data
The first frequency component, the S (x, y) are sample operator, the WA for identifying given data and missing datax,yFor positive wave
Atom Transformation operator, the Tk(f, x, y) is threshold operator, the Dk-1The reconstruct data iterated to calculate for kth -1 time, in institute
When stating k=1, D0=D (f, x, y).
6. device according to claim 5, which is characterized in that the data acquisition module is configured as:
The original earthquake data in acquisition time domain;
One-dimensional Fourier transform is done along time orientation to the original earthquake data of the time-domain, obtains the earthquake of frequency domain
Data.
7. device according to claim 5, which is characterized in that the threshold process unit is configured as:
The wave atom numeric field data is handled according to threshold formula,
The threshold formula are as follows:
The PkFor the threshold value of the kth time iterative calculation, the PkIncluded in the default seismic data threshold value P.
8. device according to claim 5, which is characterized in that described device further include:
Fourier inversion module is configured as doing the reconstruct seismic data of the first frequency component along frequency direction one-dimensional
Fourier inversion.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510885713.5A CN106842321B (en) | 2015-12-04 | 2015-12-04 | Reconstruction of seismic data method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510885713.5A CN106842321B (en) | 2015-12-04 | 2015-12-04 | Reconstruction of seismic data method and apparatus |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106842321A CN106842321A (en) | 2017-06-13 |
CN106842321B true CN106842321B (en) | 2019-02-15 |
Family
ID=59150672
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510885713.5A Active CN106842321B (en) | 2015-12-04 | 2015-12-04 | Reconstruction of seismic data method and apparatus |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106842321B (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107766293B (en) * | 2017-09-28 | 2021-04-23 | 中国计量科学研究院 | Signal spectrum analysis method and system when part of sampling data is missing regularly |
CN108387928B (en) * | 2018-02-11 | 2021-06-15 | 中国石油化工股份有限公司 | Data construction method based on seismic feature transformation space |
CN110764135A (en) * | 2018-07-26 | 2020-02-07 | 中国石油化工股份有限公司 | Irregular seismic data full-band reconstruction method |
CN110858000B (en) * | 2018-08-24 | 2021-07-02 | 中国石油天然气股份有限公司 | Seismic data reconstruction method and device, computer equipment and storage medium |
CN108983287B (en) * | 2018-10-23 | 2020-08-25 | 东华理工大学 | Curvelet transform anti-aliasing seismic data reconstruction method based on convex set projection algorithm |
CN112445649A (en) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | Seismic missing data recovery method, computer storage medium and computer system |
CN111007566B (en) * | 2019-12-27 | 2020-12-18 | 西南石油大学 | Curvature-driven diffusion full-convolution network seismic data bad channel reconstruction and denoising method |
CN114509805B (en) * | 2021-05-14 | 2023-02-28 | 中国地质大学(北京) | Vector convex set projection multi-component three-dimensional seismic data reconstruction method and device |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102436641A (en) * | 2011-10-10 | 2012-05-02 | 上海交通大学 | Automatic film restoration method based on wave atom transformation and nonparametric model |
CN104516014A (en) * | 2013-09-27 | 2015-04-15 | 中国石油天然气集团公司 | Fitting landform based wave field reconstruction method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9753166B2 (en) * | 2014-05-20 | 2017-09-05 | Cgg Services Sas | P-wave and S-wave separation of seismic data in the presence of statics and irregular geometry |
-
2015
- 2015-12-04 CN CN201510885713.5A patent/CN106842321B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102436641A (en) * | 2011-10-10 | 2012-05-02 | 上海交通大学 | Automatic film restoration method based on wave atom transformation and nonparametric model |
CN104516014A (en) * | 2013-09-27 | 2015-04-15 | 中国石油天然气集团公司 | Fitting landform based wave field reconstruction method |
Non-Patent Citations (4)
Title |
---|
Irregular seismic data reconstruction based on exponential threshold model of POCS method;Gao Jian-Jun et al.;《APPLIED GEOPHYSICS》;20100930;第7卷(第3期);第229-237页 |
基于Curvelet变换的缺失地震数据插值方法;刘国昌 等;《石油地球物理勘探》;20110430;第46卷(第2期);第237-246页 |
基于jitter采样和曲波变换的三维地震数据重建;张华 等;《地球物理学报》;20130531;第56卷(第5期);第1637-1649页 |
基于系数相关性阈值的波原子域叠前地震资料信噪分离方法;杨宁 等;《石油地球物理勘探》;20110228;第46卷(第1期);第53-57页 |
Also Published As
Publication number | Publication date |
---|---|
CN106842321A (en) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842321B (en) | Reconstruction of seismic data method and apparatus | |
CN104360393A (en) | Seismic data reconstruction method | |
CN105425301A (en) | Frequency domain three-dimensional irregular earthquake data reconstruction method | |
Zhang et al. | 3D simultaneous seismic data reconstruction and noise suppression based on the curvelet transform | |
CN107765308A (en) | Reconstruct low-frequency data frequency domain full waveform inversion method based on convolution thought Yu accurate focus | |
Chang et al. | Generative adversarial networks for seismic data interpolation | |
CN107589452B (en) | The data matching method and device of longitudinal wave and converted wave | |
CN110967734B (en) | Virtual source reconstruction method and system based on fast Fourier transform | |
Safarkhani et al. | Improving C 1 and C 3 empirical Green’s functions from ambient seismic noise in NW Iran using RMS ratio stacking method | |
López et al. | Spectral gap-based seismic survey design | |
CN110261912A (en) | The interpolation and denoising method and system of a kind of seismic data | |
CN110109179A (en) | Bandwidth compensation processing method, device and equipment | |
CN111123361A (en) | Method and device for regularly reconstructing seismic data of vertical seismic profile and storage medium | |
CN113311485B (en) | Seismic sediment characteristic enhanced filtering method and device | |
CN114488302A (en) | In-situ anisotropic ground stress field prediction method and system | |
KR101614138B1 (en) | Method for 3-dimensional moment tensor inversion by dividing modeling region | |
Pham et al. | Passive seismic signal denoising using convolutional neural network | |
US10338250B2 (en) | Method of removing incoherent noise | |
Peng et al. | A nonparametric method for automatic denoising of microseismic data | |
CN112379415A (en) | Multi-scale full-waveform inversion method and device for reconstructed low-frequency data based on downsampling | |
Al-Hashmi et al. | A wavelet transform method to detect P and S-phases in three component seismic data | |
Groos | Broadband seismic noise: classification and Green's function estimation | |
CN110361782A (en) | A kind of seismic waveform clustering method and device | |
CN108345034A (en) | A kind of seismic data rule method | |
CN115639596A (en) | Seismic data reconstruction method and device |
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 |