CN106842321A - Reconstruction of seismic data method and apparatus - Google Patents
Reconstruction of seismic data method and apparatus Download PDFInfo
- Publication number
- CN106842321A CN106842321A CN201510885713.5A CN201510885713A CN106842321A CN 106842321 A CN106842321 A CN 106842321A CN 201510885713 A CN201510885713 A CN 201510885713A CN 106842321 A CN106842321 A CN 106842321A
- Authority
- CN
- China
- Prior art keywords
- data
- reconstruct
- geological data
- geological
- 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.)
- Granted
Links
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. analysis, for interpretation, for correction
- G01V1/32—Transforming one recording into another or one representation into another
Abstract
The invention discloses a kind of Reconstruction of seismic data method and apparatus, belong to field of seismic exploration.Methods described includes：Collection geological data；Obtain default geological data threshold value；The geological data is reconstructed by convex set projection algorithm in ripple atom domain according to the default geological data threshold value, the geological data after reconstruct is obtained.The present invention is reconstructed with convex set projection algorithm by ripple atom domain come the geological data to irregular missing, solving when original earthquake data is reconstructed by warp wavelet in correlation technique, because the anisotropy of warp wavelet makes its amount of calculation larger, the less efficient problem of reconstruction of geological data.On the premise of having reached the precision of the geological data for ensureing reconstruct, the effect of computational efficiency higher is possessed when geological data is reconstructed.
Description
Technical field
The present invention relates to field of seismic exploration, more particularly to a kind of Reconstruction of seismic data method and apparatus.
Background technology
Seismic prospecting is the important technology of exploration oil and gas resource before probing, and the technology is in coalfield and work
The aspects such as journey geologic prospect, Study on Regional Geology and the earth's crust research, are also used widely.And due to environment
The limitation such as complexity, the error of instrument and financial cost, the geological data for collecting often lack sampling,
And the geological data of lack sampling can cause information irregular and missing.Irregular and missing original earthquake data
Had a strong impact on to comprehending generation at data.Accordingly, it would be desirable to be rebuild to original earthquake data.
In correlation technique when being rebuild to original earthquake data, usually using the earthquake based on sparse transformation
Data re-establishing method, in the method, is become by the Qu Bo (curvelet) in sparse transformation and brought to original
Geological data is rebuild, and warp wavelet, can due to its good directionality, locality and anisotropy
Obtain and preferably rebuild effect.
Realize it is of the invention during, inventor find prior art at least there is problems with：Above-mentioned side
Method when being reconstructed to original earthquake data by warp wavelet, because the anisotropy of warp wavelet makes it
Amount of calculation is larger, and the reconstruction of geological data is less efficient.
The content of the invention
During in order to solve in the prior art to be reconstructed original earthquake data by warp wavelet, due to Qu Bo
The anisotropy of conversion makes its amount of calculation larger, and the less efficient problem of reconstruction of geological data, the present invention is carried
A kind of Reconstruction of seismic data method and apparatus are supplied.The technical scheme is as follows：
According to the first aspect of the invention, there is provided a kind of Reconstruction of seismic data method, methods described includes：
Collection geological data；
Obtain default geological data threshold value；
Convex set projection algorithm is passed through to the geological data in ripple atom domain according to the default geological data threshold value
It is reconstructed, obtains the geological data after reconstruct.
Alternatively, the collection geological data, including：
The original earthquake data in acquisition time domain；
Original earthquake data to the time-domain does one-dimensional Fourier transform along time orientation, obtains frequency domain
The geological data.
Alternatively, the geological data includes multiple frequency components, described according to the default geological data
Threshold value is reconstructed to the geological data in ripple atom domain by convex set projection algorithm, obtains the ground after reconstruct
Shake data, including：
The first frequency component in the geological data is carried out by the convex set projection algorithm in ripple atom domain
The iterative calculation of preset times, the first frequency component is any frequency component in the geological data；
Detection iterates to calculate the signal to noise ratio of the reconstruct data for drawing every time, until there is signal to noise ratio more than default letter
During the reconstruct data for making an uproar than threshold value, stop iterative calculation；
The signal to noise ratio is more than the reconstruct data of default snr threshold as the weight of the first frequency component
Structure geological data.
Alternatively, it is described to pass through the convex set projection algorithm to first in the geological data in ripple atom domain
Frequency component carries out the iterative calculation of preset times, including：
The initial value of kth time iterative calculation is obtained, when the k is equal to 1, the initial value is for describedly
The first frequency component in shake data, when the k is more than 1, the initial value is -1 iteration of kth
The reconstruct data being calculated；
Traveling wave Atom Transformation is entered to the initial value, ripple atom numeric field data is obtained；
Threshold process is carried out to the ripple atom numeric field data by the default geological data threshold value；
Traveling wave atom inverse transformation is entered to the ripple atom numeric field data after treatment, processing data is obtained；
The initial value of the first frequency component is implanted into by the processing data by convex set projection formula and obtains institute
The reconstruct data of kth time iterative calculation are stated,
The convex set projection formula is：
The Dk(f, x, y) is the reconstruct data of the kth time iterative calculation, the D (f, x, y) for describedly
The first frequency component, the S (x, y) in shake data are that the sampling for recognizing given data and missing data is calculated
Sub, described WAx,yIt is positive wave Atom Transformation operator, the Tk(f, x, y) is threshold operator, the Dk-1It is kth -1
The reconstruct data of secondary iterative calculation, in the k=1, D0=D (f, x, y).
Alternatively, it is described the ripple atom numeric field data is carried out at threshold value by the default geological data threshold value
Reason, including：
The ripple atom numeric field data is processed according to threshold formula,
The threshold formula is：
The PkIt is the threshold value of the kth time iterative calculation, the PkIt is included in the default geological data threshold value
In P.
Alternatively, the reconstruct data that the signal to noise ratio is more than default snr threshold are used as described first
After the reconstruct geological data of frequency component, methods described also includes：
Reconstruct geological data to the first frequency component does one-dimensional Fourier inversion along frequency direction.
According to the second aspect of the invention, there is provided a kind of Reconstruction of seismic data device, described device includes：
Data acquisition module, is configured as gathering geological data；
Threshold value acquisition module, is configured as obtaining default geological data threshold value；
Data reconstruction module, is configured as according to the default geological data threshold value in ripple atom domain to described
Shake data are reconstructed by convex set projection algorithm, obtain the geological data after reconstruct.
Alternatively, the data acquisition module, is configured as：
The original earthquake data in acquisition time domain；
Original earthquake data to the time-domain does one-dimensional Fourier transform along time orientation, obtains frequency domain
The geological data.
Alternatively, the geological data includes multiple frequency components,
The data reconstruction module, including：
Iterative calculation submodule, is configured as in ripple atom domain by the convex set projection algorithm to the earthquake
First frequency component in data carries out the iterative calculation of preset times, and the first frequency component is for describedly
Any frequency component in shake data；
Signal-to-noise ratio computation submodule, is configured as detecting the signal to noise ratio for iterating to calculate the reconstruct data for drawing every time,
Until when there is signal to noise ratio more than the reconstruct data for presetting snr threshold, stopping iterative calculation；
Reconstruct data determination sub-module, is configured as being more than the signal to noise ratio reconstruct of default snr threshold
Data as the first frequency component reconstruct geological data.
Alternatively, the iterative calculation submodule, including：
Initial value acquiring unit, is configured as obtaining the initial value of kth time iterative calculation, and 1 is equal in the k
When, the initial value is the first frequency component in the geological data, when the k is more than 1,
The initial value is that kth iterates to calculate the reconstruct data for obtaining -1 time；
Ripple Atom Transformation unit, is configured as entering traveling wave Atom Transformation to the initial value, obtains ripple atom domain
Data；
Threshold process unit, is configured to the default geological data threshold value to the ripple atom numeric field data
Carry out threshold process；
Ripple atom inverse transformation unit, is configured as entering traveling wave atom inverse transformation to the ripple atom numeric field data after treatment,
Obtain processing data；
Data are implanted into unit, are configured to convex set projection formula by the initial value of the first frequency component
The reconstruct data that the processing data obtains the kth time iterative calculation are implanted into,
The convex set projection formula is：
The Dk(f, x, y) is the reconstruct data of the kth time iterative calculation, the D (f, x, y) for describedly
The first frequency component, the S (x, y) in shake data are that the sampling for recognizing given data and missing data is calculated
Sub, described WAx,yIt is positive wave Atom Transformation operator, the Tk(f, x, y) is threshold operator, the Dk-1It is kth -1
The reconstruct data of secondary iterative calculation, in the k=1, D0=D (f, x, y).
Alternatively, the threshold process unit, is configured as：
The ripple atom numeric field data is processed according to threshold formula,
The threshold formula is：
The PkIt is the threshold value of the kth time iterative calculation, the PkIt is included in the default geological data threshold value
In P.
Alternatively, described device also includes：
Fourier inversion module, is configured as to the reconstruct geological data of the first frequency component along frequency
Do one-dimensional Fourier inversion in direction.
The beneficial effect brought of technical scheme that the present invention is provided is：
It is reconstructed come the geological data to irregular missing with convex set projection algorithm by ripple atom domain, is solved
It is each due to warp wavelet when being reconstructed to original earthquake data by warp wavelet in correlation technique of having determined
Anisotropy makes its amount of calculation larger, the less efficient problem of reconstruction of geological data.Reach and ensured reconstruct
On the premise of the precision of geological data, the effect of computational efficiency higher is possessed when geological data is reconstructed.
Brief description of the drawings
Technical scheme in order to illustrate more clearly the embodiments of the present invention, institute in being described to embodiment below
The accompanying drawing for needing to use is briefly described, it should be apparent that, drawings in the following description are only the present invention
Some embodiments, for those of ordinary skill in the art, on the premise of not paying creative work,
Other accompanying drawings can also be obtained according to these accompanying drawings.
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 of iterative calculation in Fig. 2-1 illustrated embodiments；
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 embodiments；
Fig. 3-3 is the block diagram of iterative calculation submodule in Fig. 3-1 illustrated embodiments；
Fig. 3-4 is the block diagram of another Reconstruction of seismic data device provided in an embodiment of the present invention.
By above-mentioned accompanying drawing, it has been shown that the clear and definite embodiment of the disclosure, will hereinafter be described in more detail.
These accompanying drawings and word description are not intended to limit the scope that the disclosure is conceived by any mode, but logical
Cross with reference to the concept that specific embodiment is those skilled in the art's explanation disclosure.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with accompanying drawing to the present invention
Implementation method 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, and the method can be with
Including following several steps：
In a step 101, geological data is gathered.
In a step 102, default geological data threshold value is obtained.
In step 103, geological data is thrown by convex set in ripple atom domain according to default geological data threshold value
Shadow algorithm is reconstructed, and obtains the geological data after reconstruct.
In sum, Reconstruction of seismic data method provided in an embodiment of the present invention, by ripple atom domain with convex
Collection projection algorithm is reconstructed come the geological data to irregular missing, solves in correlation technique by Qu Bo
When conversion is reconstructed to original earthquake data, because the anisotropy of warp wavelet makes its amount of calculation larger,
The less efficient problem of reconstruction of geological data.On the premise of having reached the precision of the geological data for ensureing reconstruct,
Possess the effect of computational efficiency higher when geological data is reconstructed.
Fig. 2-1 is the flow chart of another Reconstruction of seismic data method provided in an embodiment of the present invention, and the method can
With including following several steps：
In step 201, the original earthquake data in acquisition time domain.
When using Reconstruction of seismic data method provided in an embodiment of the present invention, can acquisition time domain first
Original earthquake data.The method of the original earthquake data in acquisition time domain may be referred to correlation technique, exemplary
, the original earthquake data of the time-domain of collection can be the prestack three-dimensional irregular missing earthquake in a certain work area
Data, the record length of interception can be 1 second, and the sampling interval can be 1 millisecond, totally 1000 sampled points.
In step 202., the original earthquake data to time-domain does one-dimensional Fourier transform along time orientation,
Obtain the geological data of frequency domain.
After the original earthquake data for acquiring time-domain, can be to the original earthquake data of time-domain along the time
One-dimensional Fourier transform is done in direction, obtains the geological data of frequency domain.Original earthquake data is transformed into frequency
Behind domain, due to need to only calculate the data of positive frequency, thus can reduce by 50% amount of calculation in subsequent treatment.
In step 203, default geological data threshold value is obtained.
After the geological data for obtaining frequency domain, default geological data threshold value can be obtained, the default earthquake number
Can be configured in advance according to threshold value, specific set-up mode may be referred to correlation technique, the present embodiment is no longer gone to live in the household of one's in-laws on getting married
State.
In step 204, the first frequency in geological data is divided by convex set projection algorithm in ripple atom domain
Amount carries out the iterative calculation of preset times.
After the geological data for obtaining frequency domain, convex set projection algorithm can be passed through over the ground in ripple atom domain
First frequency component in shake data carries out the iterative calculation of preset times.Wherein, first frequency component is ground
Any frequency component in shake data.
As shown in Fig. 2-2, this step can include following several sub-steps：
In sub-step 2041, the initial value of kth time iterative calculation is obtained.
Preset times are carried out to the first frequency component in geological data by convex set projection algorithm in ripple atom domain
Iterative calculation when, the initial value of kth time iterative calculation can be obtained first, when k is equal to 1, initially
It is worth the first frequency component to be record in geological data, when k is more than 1, initial value is -1 iteration of kth
The reconstruct data being calculated.
In sub-step 2042, traveling wave Atom Transformation is entered to initial value, obtain ripple atom numeric field data.
When the initial value of kth time iterative calculation is obtained, traveling wave Atom Transformation can be entered to initial value, obtained
To ripple atom numeric field data.Ripple Atom Transformation is the Functional Analysis side with multi-resolution analysis, locality and directionality
Method, compared to warp wavelet, ripple Atom Transformation more can sparsely represent geological data, and more preferable Approaching earthquake is same
Phase axle.Exemplary, the mathematic(al) representation of traveling wave Atom Transformation is entered to initial value can be：
Dk(f, x, y)=WAx,yDk-1。
Wherein, DkThe ripple atom numeric field data that (f, x, y) is obtained when being kth time iterative calculation, WAx,yIt is for sky
Between variable x, y positive wave Atom Transformation operator, Dk-1It is the initial value of kth subwave Atom Transformation, in k=1
When, D0=D (f, x, y), i.e., the first frequency component that record in geological data.
In sub-step 2043, threshold process is carried out to ripple atom numeric field data by default geological data threshold value.
After ripple atom numeric field data is obtained, ripple atom numeric field data can be carried out by default geological data threshold value
Threshold process,.Optionally, ripple atom numeric field data can be processed according to threshold formula.
Threshold formula is：
Wherein, Tk(f, x, y) is threshold operator, PkIt is the threshold value of kth time iterative calculation, PkIt is included in default ground
In shake data threshold P.
In sub-step 2044, traveling wave atom inverse transformation is entered to the ripple atom numeric field data after treatment, processed
Data.
After threshold process is carried out to ripple atom numeric field data, traveling wave can be entered to the ripple atom numeric field data after treatment
Atom inverse transformation, obtains processing data.Can by the ripple atom numeric field data after treatment with for space variable x,
The ripple atom inverse transformation operator of y is multiplied to obtain processing data, exemplary, and ripple atom inverse transformation operator can be with
It is WAx,y -1。
In sub-step 2045, the initial value of first frequency component is implanted into by treatment number by convex set projection formula
According to the reconstruct data for obtaining kth time iterative calculation.
After processing data is obtained, can be by convex set projection formula by the initial value of first frequency component
Implantation processing data obtains the reconstruct data of kth time iterative calculation.And convex set projection formula can be：
Wherein, Dk(f, x, y) is the reconstruct data of kth time iterative calculation, D (f, x, y) in geological data
First frequency component, S (x, y) are sample operator, the WA for recognizing given data and missing datax,yIt is positive wave atom
Transformation operator, Tk(f, x, y) is threshold operator, Dk-1It is -1 reconstruct data of iterative calculation of kth, in k=1
When, D0=D (f, x, y).
In step 205, detection iterates to calculate the signal to noise ratio of the reconstruct data for drawing every time, until there is letter
Make an uproar during than reconstruct data more than default snr threshold, stop iterative calculation.
When calculating is iterated, the signal to noise ratio for iterating to calculate the reconstruct data for drawing every time can be detected, directly
When there is signal to noise ratio more than the reconstruct data for presetting snr threshold, stop iterative calculation.
Exemplary, signal-to-noise ratio computation formula can be：
Wherein, SNR is signal to noise ratio, Dobs(f, x, y) is the irregular missing earthquake of observation data, i.e. frequency domain
Data, Drec(f, x, y) is the reconstruction data observed after data implantation, Dcal(f, x, y) is before observation data are implanted into
Rebuild data.
It should be noted that Reconstruction of seismic data method provided in an embodiment of the present invention, compared to correlation technique
In reconstruction geological data, signal to noise ratio can lift 1 to 2 times.
In step 206, the reconstruct data for signal to noise ratio being more than default snr threshold are divided as first frequency
The reconstruct geological data of amount.
When there is signal to noise ratio more than the reconstruct data for presetting snr threshold, can be by signal to noise ratio more than default
The reconstruct data of snr threshold as first frequency component reconstruct geological data.So far, first has been obtained
The reconstruct geological data of frequency component.
In step 207, the reconstruct geological data to first frequency component is one-dimensional Fourier along frequency direction
Inverse transformation.
After the reconstruct earthquake number for obtaining first frequency component, can be to the reconstruct earthquake of first frequency component
Data do one-dimensional Fourier inversion along frequency direction, convert it to time-domain, follow-up to earthquake to facilitate
Data carry out correlation analysis.
It should be noted that can be by step 204 to step 207, to each frequency component of frequency domain
Geological data is rebuild, and obtains the reconstruct geological data of time-domain.In addition it is also possible to first pass through step 204
To step 206 obtain each frequency domain reconstruct geological data, then the frequency domain to all frequency components weight
The geological data of structure does one-dimensional Fourier inversion along frequency direction, and the embodiment of the present invention is not restricted.
It should be added that, Reconstruction of seismic data method provided in an embodiment of the present invention, by not
Before the geological data of rule missing enters traveling wave Atom Transformation, frequency domain is first converted it to, reached reduction meter
Calculation amount, lifts the effect of Reconstruction of seismic data efficiency.
In sum, Reconstruction of seismic data method provided in an embodiment of the present invention, by ripple atom domain with convex
Collection projection algorithm is reconstructed come the geological data to irregular missing, solves in correlation technique by Qu Bo
When conversion is reconstructed to original earthquake data, because the anisotropy of warp wavelet makes its amount of calculation larger,
The less efficient problem of reconstruction of geological data.On the premise of having reached the precision of the geological data for ensureing reconstruct,
Possess the effect of computational efficiency higher when geological data is reconstructed.
Following is apparatus of the present invention embodiment, can be used for performing the inventive method embodiment.For the present invention
The details not disclosed in device embodiment, refer to the inventive method embodiment.
Fig. 3-1 is a kind of block diagram of Reconstruction of seismic data device provided in an embodiment of the present invention, the geological data weight
Building device includes：
Data acquisition module 310, is configured as gathering geological data.
Threshold value acquisition module 320, is configured as obtaining default geological data threshold value.
Data reconstruction module 330, is configured as according to default geological data threshold value in ripple atom domain to geological data
It is reconstructed by convex set projection algorithm, obtains the geological data after reconstruct.
In sum, Reconstruction of seismic data device provided in an embodiment of the present invention, by ripple atom domain with convex
Collection projection algorithm is reconstructed come the geological data to irregular missing, solves in correlation technique by Qu Bo
When conversion is reconstructed to original earthquake data, because the anisotropy of warp wavelet makes its amount of calculation larger,
The less efficient problem of reconstruction of geological data.On the premise of having reached the precision of the geological data for ensureing reconstruct,
Possess the effect of computational efficiency higher when geological data is reconstructed.
Alternatively, data acquisition module 310, are configured as：
The original earthquake data in acquisition time domain；
Original earthquake data to time-domain does one-dimensional Fourier transform along time orientation, obtains the ground of frequency domain
Shake data.
Alternatively, geological data includes multiple frequency components, as shown in figure 3-2, data reconstruction module 330,
Including：
Iterative calculation submodule 331, is configured as in ripple atom domain by convex set projection algorithm in geological data
First frequency component carry out the iterative calculation of preset times, first frequency component is any in geological data
Frequency component；
Signal-to-noise ratio computation submodule 332, is configured as detecting the noise for iterating to calculate the reconstruct data for drawing every time
Than until when there is signal to noise ratio more than the reconstruct data for presetting snr threshold, stopping iterative calculation；
Reconstruct data determination sub-module 333, is configured as being more than signal to noise ratio the reconstruct number of default snr threshold
According to the reconstruct geological data as first frequency component.
Alternatively, as shown in Fig. 3-3, submodule 331 is iterated to calculate, including：
Initial value acquiring unit 331a, is configured as obtaining the initial value of kth time iterative calculation, is equal in k
When 1, initial value is the first frequency component in geological data, and when k is more than 1, initial value is kth -1 time
The reconstruct data that iterative calculation is obtained；
Ripple Atom Transformation unit 331b, is configured as entering initial value traveling wave Atom Transformation, obtains ripple atom domain
Data；
Threshold process unit 331c, be configured to default geological data threshold value is carried out to ripple atom numeric field data
Threshold process；
Ripple atom inverse transformation unit 331d, is configured as entering traveling wave atom to the ripple atom numeric field data after treatment anti-
Conversion, obtains processing data；
Data are implanted into unit 331e, are configured to convex set projection formula by the initial value of first frequency component
Implantation processing data obtains the reconstruct data of kth time iterative calculation, and convex set projection formula is：
Dk(f, x, y) is the first frequency in geological data for the secondary reconstruct data for iterating to calculate of kth, D (f, x, y)
Rate component, S (x, y) are sample operator, the WA for recognizing given data and missing datax,yFor positive wave Atom Transformation is calculated
Son, Tk(f, x, y) is threshold operator, Dk-1It is -1 reconstruct data of iterative calculation of kth, in k=1,
D0=D (f, x, y).
Alternatively, threshold process unit 331c, is configured as：
Ripple atom numeric field data is processed according to threshold formula,
Threshold formula is：
PkIt is the threshold value of kth time iterative calculation, PkIt is included in default geological data threshold value.
Alternatively, as shown in Figure 3-4, the device also includes：
Fourier inversion module 340, is configured as to the reconstruct geological data of first frequency component along frequency side
To doing one-dimensional Fourier inversion.
It should be added that, Reconstruction of seismic data method provided in an embodiment of the present invention, by not
Before the geological data of rule missing enters traveling wave Atom Transformation, frequency domain is first converted it to, reached reduction meter
Calculation amount, lifts the effect of Reconstruction of seismic data efficiency.
In sum, Reconstruction of seismic data device provided in an embodiment of the present invention, by ripple atom domain with convex
Collection projection algorithm is reconstructed come the geological data to irregular missing, solves in correlation technique by Qu Bo
When conversion is reconstructed to original earthquake data, because the anisotropy of warp wavelet makes its amount of calculation larger,
The less efficient problem of reconstruction of geological data.On the premise of having reached the precision of the geological data for ensureing reconstruct,
Possess the effect of computational efficiency higher when geological data is reconstructed.
One of ordinary skill in the art will appreciate that realize all or part of step of above-described embodiment can pass through
Hardware is completed, it is also possible to instruct the hardware of correlation to complete by program, described program can be stored in
In a kind of computer-readable recording medium, storage medium mentioned above can be read-only storage, disk or
CD etc..
Presently preferred embodiments of the present invention is the foregoing is only, is not intended to limit the invention, it is all of the invention
Within spirit and principle, any modification, equivalent substitution and improvements made etc. should be included in of the invention
Within protection domain.
Claims (12)
1. a kind of Reconstruction of seismic data method, it is characterised in that methods described includes：
Collection geological data；
Obtain default geological data threshold value；
Convex set projection algorithm is passed through to the geological data in ripple atom domain according to the default geological data threshold value
It is reconstructed, obtains the geological data after reconstruct.
2. method according to claim 1, it is characterised in that the collection geological data, including：
The original earthquake data in acquisition time domain；
Original earthquake data to the time-domain does one-dimensional Fourier transform along time orientation, obtains frequency domain
The geological data.
3. method according to claim 2, it is characterised in that the geological data includes multiple frequencies
Rate component, it is described that convex set is passed through to the geological data in ripple atom domain according to the default geological data threshold value
Projection algorithm is reconstructed, and obtains the geological data after reconstruct, including：
The first frequency component in the geological data is carried out by the convex set projection algorithm in ripple atom domain
The iterative calculation of preset times, the first frequency component is any frequency component in the geological data；
Detection iterates to calculate the signal to noise ratio of the reconstruct data for drawing every time, until there is signal to noise ratio more than default letter
During the reconstruct data for making an uproar than threshold value, stop iterative calculation；
The signal to noise ratio is more than the reconstruct data of default snr threshold as the weight of the first frequency component
Structure geological data.
4. method according to claim 3, it is characterised in that it is described in ripple atom domain by described convex
Collection projection algorithm carries out the iterative calculation of preset times to the first frequency component in the geological data, including：
The initial value of kth time iterative calculation is obtained, when the k is equal to 1, the initial value is for describedly
The first frequency component in shake data, when the k is more than 1, the initial value is -1 iteration of kth
The reconstruct data being calculated；
Traveling wave Atom Transformation is entered to the initial value, ripple atom numeric field data is obtained；
Threshold process is carried out to the ripple atom numeric field data by the default geological data threshold value；
Traveling wave atom inverse transformation is entered to the ripple atom numeric field data after treatment, processing data is obtained；
The initial value of the first frequency component is implanted into by the processing data by convex set projection formula and obtains institute
The reconstruct data of kth time iterative calculation are stated,
The convex set projection formula is：
The Dk(f, x, y) is the reconstruct data of the kth time iterative calculation, the D (f, x, y) for describedly
The first frequency component, the S (x, y) in shake data are that the sampling for recognizing given data and missing data is calculated
Sub, described WAx,yIt is positive wave Atom Transformation operator, the Tk(f, x, y) is threshold operator, the Dk-1It is kth -1
The reconstruct data of secondary iterative calculation, in the k=1, D0=D (f, x, y).
5. method according to claim 4, it is characterised in that described by the default geological data
Threshold value carries out threshold process to the ripple atom numeric field data, including：
The ripple atom numeric field data is processed according to threshold formula,
The threshold formula is：
The PkIt is the threshold value of the kth time iterative calculation, the PkIt is included in the default geological data threshold value
In P.
6. method according to claim 3, it is characterised in that it is described by the signal to noise ratio more than default
After reconstruct geological data of the reconstruct data of snr threshold as the first frequency component, methods described
Also include：
Reconstruct geological data to the first frequency component does one-dimensional Fourier inversion along frequency direction.
7. a kind of Reconstruction of seismic data device, it is characterised in that described device includes：
Data acquisition module, is configured as gathering geological data；
Threshold value acquisition module, is configured as obtaining default geological data threshold value；
Data reconstruction module, is configured as according to the default geological data threshold value in ripple atom domain to described
Shake data are reconstructed by convex set projection algorithm, obtain the geological data after reconstruct.
8. device according to claim 7, it is characterised in that the data acquisition module, is configured
For：
The original earthquake data in acquisition time domain；
Original earthquake data to the time-domain does one-dimensional Fourier transform along time orientation, obtains frequency domain
The geological data.
9. device according to claim 8, it is characterised in that the geological data includes multiple frequencies
Rate component,
The data reconstruction module, including：
Iterative calculation submodule, is configured as in ripple atom domain by the convex set projection algorithm to the earthquake
First frequency component in data carries out the iterative calculation of preset times, and the first frequency component is for describedly
Any frequency component in shake data；
Signal-to-noise ratio computation submodule, is configured as detecting the signal to noise ratio for iterating to calculate the reconstruct data for drawing every time,
Until when there is signal to noise ratio more than the reconstruct data for presetting snr threshold, stopping iterative calculation；
Reconstruct data determination sub-module, is configured as being more than the signal to noise ratio reconstruct of default snr threshold
Data as the first frequency component reconstruct geological data.
10. device according to claim 9, it is characterised in that the iterative calculation submodule, bag
Include：
Initial value acquiring unit, is configured as obtaining the initial value of kth time iterative calculation, and 1 is equal in the k
When, the initial value is the first frequency component in the geological data, when the k is more than 1,
The initial value is that kth iterates to calculate the reconstruct data for obtaining -1 time；
Ripple Atom Transformation unit, is configured as entering traveling wave Atom Transformation to the initial value, obtains ripple atom domain
Data；
Threshold process unit, is configured to the default geological data threshold value to the ripple atom numeric field data
Carry out threshold process；
Ripple atom inverse transformation unit, is configured as entering traveling wave atom inverse transformation to the ripple atom numeric field data after treatment,
Obtain processing data；
Data are implanted into unit, are configured to convex set projection formula by the initial value of the first frequency component
The reconstruct data that the processing data obtains the kth time iterative calculation are implanted into,
The convex set projection formula is：
The Dk(f, x, y) is the reconstruct data of the kth time iterative calculation, the D (f, x, y) for describedly
The first frequency component, the S (x, y) in shake data are that the sampling for recognizing given data and missing data is calculated
Sub, described WAx,yIt is positive wave Atom Transformation operator, the Tk(f, x, y) is threshold operator, the Dk-1It is kth -1
The reconstruct data of secondary iterative calculation, in the k=1, D0=D (f, x, y).
11. devices according to claim 10, it is characterised in that the threshold process unit, are matched somebody with somebody
It is set to：
The ripple atom numeric field data is processed according to threshold formula,
The threshold formula is：
The PkIt is the threshold value of the kth time iterative calculation, the PkIt is included in the default geological data threshold value
In P.
12. devices according to claim 9, it is characterised in that described device also includes：
Fourier inversion module, is configured as to the reconstruct geological data of the first frequency component along frequency
Do one-dimensional Fourier inversion in direction.
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 true CN106842321A (en) | 2017-06-13 |
| CN106842321B 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) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107766293A (en) * | 2017-09-28 | 2018-03-06 | 中国计量科学研究院 | Signal spectral analysis method and system when fractional-sample data rule lacks |
| CN108387928A (en) * | 2018-02-11 | 2018-08-10 | 中国石油化工股份有限公司 | Data construction method based on seismic signature transformation space |
| CN108983287A (en) * | 2018-10-23 | 2018-12-11 | 东华理工大学 | A kind of anti-alias Reconstruction of seismic data method of warp wavelet based on convex set projection algorithm |
| CN110764135A (en) * | 2018-07-26 | 2020-02-07 | 中国石油化工股份有限公司 | Irregular seismic data full-band reconstruction method |
| CN110858000A (en) * | 2018-08-24 | 2020-03-03 | 中国石油天然气股份有限公司 | Seismic data reconstruction method and device, computer equipment and storage medium |
| CN111007566A (en) * | 2019-12-27 | 2020-04-14 | 西南石油大学 | Curvature-driven diffusion full-convolution network seismic data bad channel reconstruction and denoising method |
| CN112445649A (en) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | Seismic missing data recovery method, computer storage medium and computer system |
| CN114509805A (en) * | 2021-05-14 | 2022-05-17 | 中国地质大学（北京） | Vector convex set projection multi-component three-dimensional seismic data reconstruction method and device |
Citations (3)
| 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 |
| US20150338536A1 (en) * | 2014-05-20 | 2015-11-26 | Cgg Services Sa | 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 (3)
| 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 |
| US20150338536A1 (en) * | 2014-05-20 | 2015-11-26 | Cgg Services Sa | P-wave and s-wave separation of seismic data in the presence of statics and irregular geometry |
Non-Patent Citations (4)
| Title |
|---|
| GAO JIAN-JUN ET AL.: "Irregular seismic data reconstruction based on exponential threshold model of POCS method", 《APPLIED GEOPHYSICS》 * |
| 刘国昌 等: "基于Curvelet变换的缺失地震数据插值方法", 《石油地球物理勘探》 * |
| 张华 等: "基于jitter采样和曲波变换的三维地震数据重建", 《地球物理学报》 * |
| 杨宁 等: "基于系数相关性阈值的波原子域叠前地震资料信噪分离方法", 《石油地球物理勘探》 * |
Cited By (13)
| 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 |
| CN107766293A (en) * | 2017-09-28 | 2018-03-06 | 中国计量科学研究院 | Signal spectral analysis method and system when fractional-sample data rule lacks |
| CN108387928A (en) * | 2018-02-11 | 2018-08-10 | 中国石油化工股份有限公司 | Data construction method based on seismic signature transformation space |
| 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 |
| CN110858000A (en) * | 2018-08-24 | 2020-03-03 | 中国石油天然气股份有限公司 | Seismic data reconstruction method and device, computer equipment and storage medium |
| 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 |
| CN108983287A (en) * | 2018-10-23 | 2018-12-11 | 东华理工大学 | A kind of anti-alias Reconstruction of seismic data method of warp wavelet 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 |
| CN111007566A (en) * | 2019-12-27 | 2020-04-14 | 西南石油大学 | Curvature-driven diffusion full-convolution network seismic data bad channel reconstruction and denoising method |
| CN114509805A (en) * | 2021-05-14 | 2022-05-17 | 中国地质大学（北京） | Vector convex set projection multi-component three-dimensional seismic data reconstruction method and device |
| CN114509805B (en) * | 2021-05-14 | 2023-02-28 | 中国地质大学（北京） | Vector convex set projection multi-component three-dimensional seismic data reconstruction method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106842321B (en) | 2019-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106842321A (en) | Reconstruction of seismic data method and apparatus | |
| CN103293551B (en) | A kind of based on model constrained impedance inversion approach and system | |
| CN107422379A (en) | Multiple dimensioned seismic full-field shape inversion method based on local auto-adaptive convexification method | |
| CN103869364B (en) | Multiple wave suppression method based on dual parabolic Radon transformation | |
| Mborah et al. | Enhancing manual P-phase arrival detection and automatic onset time picking in a noisy microseismic data in underground mines | |
| CN103163554A (en) | Self-adapting wave form retrieval method through utilization of zero offset vertical seismic profile (VSP) data to estimate speed and Q value | |
| CN106707334B (en) | A method of improving seismic data resolution | |
| CN103364832A (en) | Seismic attenuation qualitative estimation method based on self-adaptive optimal kernel time frequency distribution | |
| Sabbione et al. | Microseismic data denoising via an apex-shifted hyperbolic Radon transform | |
| CN107179550B (en) | A kind of seismic signal zero phase deconvolution method of data-driven | |
| CN103149592A (en) | Method for separating variable offset vertical seismic profile (VSP) wave fields | |
| CN104932010A (en) | Diffraction wave separating method based on shortcut edging sparse Radon transformation | |
| CN106680874A (en) | Harmonic noise suppression method based on waveform morphology sparse modeling | |
| CN107678064A (en) | A kind of interval transit time real time extracting method | |
| CN107356964A (en) | Q value estimation and compensation method of the S-transformation domain based on variation principle | |
| CN104614769A (en) | Beam-forming filtering method for suppressing seismic surface waves | |
| CN104730576A (en) | Curvelet transform-based denoising method of seismic signals | |
| CN107402405B (en) | Quiet phase virtual source trace gather construction method | |
| CN105259579A (en) | A high-amplitude shielding layer rejecting method based on seismic data instantaneous attributes | |
| CN106526678A (en) | Reflection acoustic logging wave field separation method and device | |
| CN104635264B (en) | The processing method of earthquake data before superposition and equipment | |
| CN104391324A (en) | Seismic trace set dynamic correction stretching correction pre-processing technology before AVO inversion depending on frequency | |
| CN102338886B (en) | Polarization filtering method for effectively attenuating surface waves in three-component seismic record | |
| KR101923166B1 (en) | Method for correcting amplitude magnitude of seismic signal extracted from seismic ambient noise | |
| Rosa-Cintas et al. | Polarization analysis in the stationary wavelet packet domain: Application to HVSR method |
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 |