CN105425301A - Frequency domain three-dimensional irregular earthquake data reconstruction method - Google Patents
Frequency domain three-dimensional irregular earthquake data reconstruction method Download PDFInfo
- Publication number
- CN105425301A CN105425301A CN201610008843.5A CN201610008843A CN105425301A CN 105425301 A CN105425301 A CN 105425301A CN 201610008843 A CN201610008843 A CN 201610008843A CN 105425301 A CN105425301 A CN 105425301A
- Authority
- CN
- China
- Prior art keywords
- data
- reconstruction
- represent
- threshold
- frequency field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 230000001788 irregular Effects 0.000 title claims abstract description 20
- 230000008569 process Effects 0.000 claims abstract description 13
- 230000000694 effects Effects 0.000 claims description 15
- 230000008034 disappearance Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000012804 iterative process Methods 0.000 claims description 6
- 230000007812 deficiency Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000005428 wave function Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract 2
- 238000005070 sampling Methods 0.000 description 24
- 238000010586 diagram Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000010183 spectrum analysis Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000007115 recruitment 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
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 brings forward a frequency domain three-dimensional irregular earthquake data reconstruction method. The method is characterized in that first of all, three-dimensional earthquake data in a time domain is converted to a frequency domain by use of Fourier transform, and then, a projection onto convex set (POCS) algorithm is employed and curvelet transform capable of describing localized features of the earthquake data is introduced; and in an iteration process, a new threshold parameter attenuating according to an index rule is brought forward, and each frequency slice is individually reconstructed by use of a soft threshold operator, such that the iteration frequency is reduced, the reconstruction precision is improved, and the purpose of reconstructing the three-dimensional earthquake data is realized. According to the invention, the new threshold parameter attenuating according to the index rule is brought forward and each frequency slice is reconstructed in individually by use of the soft threshold operator, such that the disadvantage of quite slow convergence speed of a conventional threshold parameter is overcome, the calculation complexity of an algorithm is reduced, the calculation efficiency is substantially improved, and the operation time is reduced.
Description
Technical field
What the present invention relates to is a kind of method lacking Reconstruction of seismic data, specifically the three-dimensional irregular earthquake data re-establishing method of a kind of frequency field.
Technical background
In the wild in data acquisition, due to collecting device, the restriction of project cost and topographic condition, geological data carries out irregular sampling along direction in space being usually less than nyquist sampling theorem requirement, thus cause the geological data that collects irregular, imperfect, there is spatial aliasing, be difficult to the requirement meeting subsequent treatment, such as velocity analysis, Free Surface Multiple attenuation, migration etc., in this case, need to develop good Reconstruction of seismic data method, reconstruct the seismic trace of disappearance, be met the geological data of processing requirements, and the method also can instruct field data collection, compression earthquake data acquisition amount.
Data re-establishing method based on warp wavelet does not need the prior imformation of underground structure, and can reconstruct the seismic trace of rule disappearance and irregular disappearance, and computing velocity is fast, precision is high, is widely used.Xue Nian etc. adopt threshold value process of iteration, achieve two-dimension earthquake data reconstruction (the geological data interpolation converted based on Curvelet and the denoising [Master's thesis] based on warp wavelet, Southwest Jiaotong University, 2010), Liu Guochang etc. adopt convex set projection algorithm, bent wave zone achieve based on POCS algorithm two-dimension earthquake data reconstruction (based on curvelet conversion disappearance geological data interpolation method [J]. geophysical prospecting for oil, 2011,46 (2), 237-245.), the method more can rebuild non-linear lineups, more be applicable to the seismic wave field with anisotropic feature, but do not have to further investigate the three dimension data reconstruct method based on warp wavelet, and select also not study further to threshold parameter, only adopt index threshold parameter equation, its rebuild after effect or limited.And posttension China wait achieve on this basis based on the 3D seismic data of warp wavelet rebuilds (based on jitter sampling and warp wavelet 3D seismic data reconstruction [J]. Chinese Journal of Geophysics, 2013,56 (5): 1637-1649.), by successively rebuilding time slice, obtaining and rebuilding effect preferably, and inquired into threshold parameter selection strategy.Wang Benfeng etc. propose bent ripple method for reconstructing under different threshold value (POCS combines bent wave zone Reconstruction of seismic data method [J] of Jitter sampling theory of improvement. Chinese Journal of Geophysics, 2015,50 (1): 20-27.), but the method is rebuild in time domain, and only adopt hard-threshold parameter to process, when data reconstruction comprises high and low-yield amplitude inphase axle simultaneously, because POCS algorithm is in front iterative process several times, lower than the lineups coefficient below threshold parameter by filtering, therefore these methods are rebuilding the limited efficiency of short arc significant wave lineups, the low-yield lineups of seismic signal can not be recovered completely and improve the signal to noise ratio (S/N ratio) in local anomaly district, and iterations is many, reconstruction efficiency is low, constrains the further application of the method.
Summary of the invention
The object of the invention is in order to reconstructing of quick high accuracy road 3D seismic data can be lacked, and provide the three-dimensional irregular earthquake data re-establishing method of a kind of frequency field.
Technical solution of the present invention:
The present invention proposes the three-dimensional irregular earthquake data re-establishing method of frequency field, first utilize Fourier transform that time domain 3D seismic data is transformed to frequency field, then adopt convex set projection algorithm idea and introduce the warp wavelet can portraying geological data local characteristic.Propose the new threshold parameter of exponentially decaying in an iterative process, adopt soft-threshold operator, each frequency slice is rebuild separately, reduce iterations and improve reconstruction precision, thus reaching the object of the irregular geological data of reconstruction of three-dimensional.
The three-dimensional irregular earthquake data re-establishing method of a kind of frequency field, the Problems of Reconstruction of geological data is described as recovering complete data by one group of deficiency of data by the effect of linear operator, supposes as lower linear forward model:
y
obs=Mf
Here y
obs∈ R
nthe geological data that representative gathers; F ∈ R
n, and N>=n, represents to be reconstructed without alias data; M ∈ R
n × Nrepresent diagonal matrix, its element 1 and 0 represents known seismic trace and unknown seismic trace respectively, and suppose that coefficient x is the rarefaction representation of f in bent wave zone C, then above-mentioned equation is:
Y
obs=Ax and
Here subscript
hrepresent associate matrix, as the data y from collection
obswhen middle reconstruction is without alias data f, because x is sparse, thus Corresponding Sparse Algorithm can be adopted to solve this underdetermined equation.
After sparse promotion inverting, reconstruction signal by
determine, simultaneously
In this expression formula,
represent estimated value, L
1norm is defined as
x [i] is i-th element in vector x, and by solving above-mentioned equation, the original data without alias are just rebuild out;
In Reconstruction of seismic data process, adopt convex set projection algorithm, concrete reconstruction procedures is as follows:
Step 1: input has the 3D seismic data in disappearance road in the time domain, then adopts Fourier transform, data reconstruction is transformed to frequency field from time domain;
Step 2: adopt warp wavelet to carry out sparse transformation to frequency field 3D seismic data, obtain bent wave system number in bent wave zone, select suitable threshold parameter formula λ according to the size of bent wave system number
i(i=1,2,3, N, wherein N is iterations) and formula;
Step 3: in bent wave zone, adopts soft-threshold operator to process, is also namely greater than threshold value λ
ibent wave system number deduct the value of a threshold size, be less than threshold value-λ
ibent wave system number add the value of a threshold size, and other bent wave system number zero setting;
Step 4: the bent wave system number after thresholding is done contrary flexure wave conversion and obtains time domain geological data, and then with disappearance geological data do not lack seismic trace be filled into inverse transformation after geological data in go;
Step 5: finally the geological data that obtain is substituted into step (2), re-start iteration, until run N end, then does inverse-Fourier transform to the geological data after iteration N time and namely obtains final reconstructed results.
Further, described warp wavelet is defined as:
In formula: φ
j, l, krepresent bent wave function, j, l, k represent yardstick respectively, direction and location parameter, and f (x) is geological data, and its frequency field definition is:
The bent wave system number obtained after conversion, available C{j}{l} (k
1, k
2) represent its structure, wherein j represents yardstick, and l represents direction, (k
1, k
2) represent matrix coefficient on j yardstick l direction.
Further, described soft-threshold operator, its expression formula is as follows:
T
krepresent soft-threshold operator, S
k-1represent the bent wave system number of the data reconstruction that kth-1 iteration obtains, meet
f
tand F
-1 trepresent the positive inverse-Fourier transform about time variable t,
represent the time domain data reconstruction that kth-1 iteration obtains, wherein
λ represents that N ties up threshold value set, λ={ λ
1, λ
2, λ
n, and meet λ
1> λ
2> > λ
n, N represents maximum iteration time.
Further, threshold parameter formula, its expression formula is as follows:
Wherein Max is | the maximal value of CFy|, the i.e. maximal value of warp wavelet absolute coefficient.N is total iterations, and ε is the little value close to zero, relevant with the energy of noise in data.
Advantage of the present invention: time domain 3D seismic data is transformed to frequency field by utilizing Fourier transform by the present invention, and adopt convex set projection algorithm to carry out data reconstruction with the warp wavelet can portraying geological data local characteristic, thus avoid low-yield significant wave lineups coefficient by filtering, improve disappearance Reconstruction of seismic data precision, the new threshold parameter of exponentially decaying is proposed in an iterative process simultaneously, adopt soft-threshold operator, each frequency slice is rebuild separately, iterations can be reduced, reduce the computation complexity of algorithm, improve counting yield significantly, save operation time.
Accompanying drawing explanation
Fig. 1 represents the three-dimensional irregular earthquake data re-establishing method process flow diagram in embodiment of the present invention medium frequency territory.
The different threshold parameter of Fig. 2 rebuilds signal to noise ratio (S/N ratio) curve comparison figure.
Fig. 3 is two-dimensional random lack sampling and reconstructed results comparison diagram thereof.
Fig. 4 is for rebuilding front and back 2-d spectrum comparison diagram.
Fig. 5 is signal to noise ratio (S/N ratio) and sampling rate relation curve comparison diagram.
Fig. 6 is noisy data reconstruction processes figure.
Embodiment
Following case study on implementation for illustration of the present invention, but is not used for limiting the scope of the invention.
Embodiment 1
The step realizing the method mainly comprises, the structure of Reconstructed equation, and time domain transforms to frequency field, convex set projection reconstruction algorithm, threshold parameter process etc.Concrete steps are as follows:
Step 1: the structure of Reconstructed equation.The Problems of Reconstruction of geological data can be described as recovering complete data by one group of deficiency of data by the effect of linear operator, supposes as lower linear forward model
y
obs=Mf
Here y
obs∈ R
nthe geological data that representative gathers; F ∈ R
n, and N>=n, represents to be reconstructed without alias data; M ∈ R
n × Nrepresent diagonal matrix, its element 1 and 0 represents known seismic trace and unknown seismic trace respectively, and suppose that coefficient x is the rarefaction representation of f in bent wave zone C, then above-mentioned equation can be write as
Y
obs=Ax and
Here subscript
hrepresent associate matrix, when from the data y collected
obswhen middle reconstruction is without alias data f, because x is sparse, thus this underdetermined equation is made to have solution.
After sparse promotion inverting, reconstruction signal by
determine, simultaneously
In this expression formula,
represent estimated value, L
1norm is defined as
x [i] is i-th element in vector x, and by solving above-mentioned underdetermined equation, the original complete and geological data of rule just can be rebuild out.
Step 2: time domain transforms to frequency field.The present invention adopts Fourier transform to calculate time domain geological data, and its discrete Fourier transform (DFT) formula is as follows:
In formula: m, n are integer.m=0,1,…M-1;n=0,1,…N-1。N △ t=T is the window length of spectrum analysis.
Step 3: convex set projection reconstruction algorithm.The shortcoming low for traditional reconstruction algorithm reconstruction precision and computing velocity is slow, propose the algorithm carrying out rebuilding in frequency field, and adopt soft-threshold operator and the new threshold parameter formula of exponentially decaying, its algorithm iteration expression formula is:
Wherein, f
krepresent the time domain data reconstruction that kth time iteration obtains,
represent that acquired original is to data y
obs(t, x, y), meets
f
tand F
-1 trepresent the positive inverse-Fourier transform about time variable t, C and C
-1represent positive and negative warp wavelet, its warp wavelet is defined as:
In formula: φ
j, l, krepresent bent wave function, j, l, k represent yardstick respectively, direction and location parameter, and f (x) is geological data, and its frequency field definition is:
The bent wave system number obtained after conversion, available C{j}{l} (k
1, k
2) represent its structure, wherein j represents yardstick, and l represents direction, (k
1, k
2) represent matrix coefficient on j yardstick l direction.
T
krepresent soft-threshold operator, its element meets:
S
k-1represent the bent wave system number of the data reconstruction that kth-1 iteration obtains, meet
λ represents that N ties up threshold value set, λ={ λ
1, λ
2, λ
n, and meet λ
1> λ
2> > λ
n, N represents maximum iteration time.
Step 4: threshold parameter formula.In an iterative process, its thresholding parameter value generally changes from big to small, also namely most of bent wave system number is removed in front iteration several times, only leave the lineups that energy is stronger, and in the end several times in iteration, the bent wave system number of the overwhelming majority retains, and only removes the random noise that some energy are less, i.e. threshold parameter λ
imeet:
||CFy||
∞=λ
1>λ
2>···>ε
In formula, ε is the little value close to zero, relevant with the energy of noise in data, different pieces of information ε value difference to some extent.The present invention on the basis of linear threshold parameter equation and index threshold parameter equation in the past, in conjunction with frequency field warp wavelet feature, propose to select by
the threshold parameter of rule decay, under the prerequisite ensureing precision, can improve speed of convergence sooner, this threshold parameter formula is:
Wherein Max is | the maximal value of CFy|, the i.e. maximal value of warp wavelet absolute coefficient.N is total iterations, and ε is the minimal value close to zero, relevant with the energy of noise in data.
Realizing the method concrete operations is:
(1) threshold parameter comparative analysis
In order to correlation data rebuilds effect better, definition signal to noise ratio snr=20log
10|| f
0||
2/ || f-f
0||
2, unit is dB, wherein f
0represent primary model data, f represents reconstructed results, and signal to noise ratio (S/N ratio) is higher, represent reconstructed results and model data is more close, effect is unreasonable to be thought, the present invention is simultaneously in 3D seismic data process of reconstruction, the scale parameter of warp wavelet is 5, and the angle number on the thickest yardstick is 8.
Then the threshold parameter formula adopting conventional linear threshold value, index threshold and the present invention to propose carries out data reconstruction to the data (Fig. 3 (c) represents the time slice two-dimensional random lack sampling figure of 3D seismic data) after theoretical model 50% lack sampling.Reconstructed results is as shown in Figure 2 a (when Fig. 2 (a) represents iterations 50, each iterations and signal to noise ratio (S/N ratio) figure), Fig. 2 (a) for maximum iteration time be 50 times time, each iterations threshold parameter different from three kinds rebuilds rear Between Signal To Noise Ratio curve map, after can finding out the threshold parameter reconstruction that in each iterative process, the present invention proposes, signal to noise ratio (S/N ratio) is the highest, next is index threshold parameter, is finally only linear threshold parameter.Maximum iteration time is changed from 5 ~ 100 times simultaneously, Between Signal To Noise Ratio curve map after calculating each maximum iteration time and rebuilding, as Suo Shi Fig. 2 (b) (Fig. 2 (b) represents maximum iteration time and its signal to noise ratio (S/N ratio) figure of each computing), therefrom also can find out in different maximum iteration time, it is maximum that the threshold parameter that the present invention proposes rebuilds rear signal to noise ratio (S/N ratio), and can find out, wanting to obtain signal to noise ratio (S/N ratio) is the reconstructed results of 20dB, linear threshold needs iteration about 70 times, index threshold needs iteration about 17 times, and the threshold parameter that the present invention proposes only needs iteration about 13 times, certain computing time can be saved.And for the index square root threshold parameter that index threshold and the present invention propose, after iterations is greater than 30 times, signal to noise ratio (S/N ratio) recruitment is relatively less with iterations increase, and therefore, the three-dimensional irregular earthquake data re-establishing method of frequency field of the present invention all selects 30 iteration.
(2) data reconstruction processes
Adopt sound wave finite difference method, simulate the Seismic forward record under uniform grid sampling, and just drill geological data by wave detector to obtained, shot point and time carry out being arranged in 3-D data volume.Ideal data are (Fig. 3 (c) represents the time slice two-dimensional random lack sampling figure of 3D seismic data) as Suo Shi Fig. 3 (c), wherein time slice is 0.48s, shot record migration respective distances is 1512m (128 big gun), and the distance that common receiver is corresponding is 1512m (128 road).Then the three-dimensional irregular earthquake data re-establishing method of frequency field of the present invention is adopted to rebuild it, rebuild as shown in Figure 3 d (Fig. 3 (d) represents the two-dimensional random lack sampling reconstructed results of cutting into slices sometime), signal to noise ratio (S/N ratio) after reconstruction is 31.67dB, and reconstruction precision is higher.From the analysis of reconstruction surrounding time section 2-d spectrum, (the 2-d spectrum analysis chart of a certain section of Fig. 4 (e) original earthquake data as shown in Figure 4; Fig. 4 (f) adopts the 2-d spectrum analysis chart of the inventive method reconstructed results), 2-d spectrum after the inventive method reconstruction is closest to the 2-d spectrum of raw data, this also illustrates the inventive method can reflect seismic event before variation characteristic, remain weak significant wave lineups, improve the signal to noise ratio (S/N ratio) in local anomaly district.
(3) different sampling rate compares
In order to contrast the reconstruction effect after a peacekeeping two dimension lack sampling, random lack sampling rate is kept to increase progressively from 10 ~ 80%, and then rebuild, signal to noise ratio (S/N ratio) after record reconstruction and the relation between sampling rate, as Suo Shi Fig. 5 (g) (Fig. 5 (g) represents the sampling rate after adopting warp wavelet to rebuild and Between Signal To Noise Ratio figure), therefrom can find out that sampling rate is fewer, effect after reconstruction is poor, but comparatively speaking, the effect that two dimension lack sampling is rebuild rebuilds better effects if than one dimension lack sampling, therefore need to adopt three-dimensional method for reconstructing, increase the information in other directions, to make up the deficiency of low-dimensional method for reconstructing.Simultaneously in order to the superiority of outstanding warp wavelet, itself and the three dimension data reconstruct algorithm based on Fourier transform are compared, Fourier transform also adopts 30 iteration when rebuilding, as Suo Shi Fig. 5 (h) (Fig. 5 (h) represents the sampling rate after adopting Fourier transform to rebuild and Between Signal To Noise Ratio figure), can find out, when adopting identical sampling rate, the effect that warp wavelet is rebuild is far superior to Fourier basis, this also illustrates warp wavelet more can reflect seismic event before variation characteristic, improve the reconstruction precision in local anomaly district.Also can find out, no matter be warp wavelet or Fourier transform, the reconstruction effect of its 3D seismic data is all better than two-dimension earthquake data reconstruction effect simultaneously.
(4) noisy data reconstruction
Because geological data all contains noise usually, so the noise resisting ability of the inventive method under needing to check two-dimensional random lack sampling, gaussian random noise is added in Fig. 3 (c) time slice, as Suo Shi Fig. 6 (i) (a certain noisy slice map that Fig. 6 (i) is 3D seismic data), then carry out anti-noise and rebuild simulated experiment, two-dimensional random lack sampling result is (Fig. 6 (j) is two-dimensional random lack sampling slice map sometime) as Suo Shi Fig. 6 (j), then the inventive method is adopted to rebuild it, its reconstructed results is (Fig. 6 (k) adopts the inventive method reconstructed results figure) as Suo Shi Fig. 6 (k), signal to noise ratio (S/N ratio) 10.92dB, it rebuilds front and back error section (Fig. 6 (l) rebuilds front and back Error Graph) as Suo Shi Fig. 6 (l), can be found out by contrast, effective information after reconstruction is almost constant, data reconstruction effect is better, this also illustrates that the inventive method has good noise resisting ability in Reconstruction of seismic data, can apply in the process of real data.
Claims (6)
1. the three-dimensional irregular earthquake data re-establishing method of frequency field, it is characterized in that, first utilize Fourier transform that time domain 3D seismic data is transformed to frequency field, then adopt convex set projection algorithm idea and introduce the warp wavelet can portraying geological data local characteristic; Propose the new threshold parameter of exponentially decaying in an iterative process, adopt soft-threshold operator, each frequency slice is rebuild separately, thus reduce iterations and improve reconstruction precision, reach the object of reconstruction of three-dimensional geological data.
2. the three-dimensional irregular earthquake data re-establishing method of a kind of frequency field according to claim 1, it is characterized in that, the Problems of Reconstruction of geological data is described as recovering complete data by one group of deficiency of data by the effect of linear operator, supposes as lower linear forward model:
y
obs=Mf
Here y
obs∈ R
nthe geological data that representative gathers; F ∈ R
n, and N>=n, represents to be reconstructed without alias data; M ∈ R
n × Nrepresent diagonal matrix, its element 1 and 0 represents known seismic trace and unknown seismic trace respectively, and suppose that coefficient x is the rarefaction representation of f in bent wave zone C, then above-mentioned equation is:
Y
obs=Ax and
Here subscript
hrepresent associate matrix, as the data y from collection
obswhen middle reconstruction is without alias data f, because x is sparse, thus Corresponding Sparse Algorithm can be adopted to solve this underdetermined equation;
After sparse promotion inverting, reconstruction signal by
determine, simultaneously
In this expression formula,
represent estimated value, L
1norm is defined as
x [i] is i-th element in vector x, and by solving above-mentioned equation, the original data without alias are just rebuild out.
3. the three-dimensional irregular earthquake data re-establishing method of a kind of frequency field according to claim 1 and 2, is characterized in that, in Reconstruction of seismic data process, adopt convex set projection algorithm, concrete reconstruction procedures is as follows:
Step 1: input has the 3D seismic data in disappearance road in the time domain, then adopts Fourier transform, data reconstruction is transformed to frequency field from time domain;
Step 2: adopt warp wavelet to carry out sparse transformation to frequency field 3D seismic data, obtain bent wave system number in bent wave zone, select suitable threshold parameter formula λ according to the size of bent wave system number
i(i=1,2,3 ..., N, wherein N is iterations);
Step 3: in bent wave zone, adopts soft-threshold operator to process, is also namely greater than threshold value λ
ibent wave system number deduct the value of a threshold size, be less than threshold value-λ
ibent wave system number add the value of a threshold size, and other bent wave system number zero setting;
Step 4: the bent wave system number after thresholding is done contrary flexure wave conversion and obtains time domain geological data, and then with disappearance geological data do not lack seismic trace be filled into inverse transformation after geological data in go;
Step 5: finally the geological data that obtain is substituted into step (2), re-start iteration, until run N end, then does inverse-Fourier transform to the geological data after iteration N time and namely obtains final reconstructed results.
4. the three-dimensional irregular earthquake data re-establishing method of a kind of frequency field according to claim 3, it is characterized in that, described warp wavelet is defined as:
In formula: φ
j, l, krepresent bent wave function, j, l, k represent yardstick respectively, direction and location parameter, and f (x) is geological data, and its frequency field definition is:
The bent wave system number obtained after conversion, available C{j}{l} (k
1, k
2) represent its structure, wherein j represents yardstick, and l represents direction, (k
1, k
2) represent matrix coefficient on j yardstick l direction.
5. the three-dimensional irregular earthquake data re-establishing method of a kind of frequency field according to claim 3, it is characterized in that, described soft-threshold operator, its expression formula is as follows:
T
krepresent soft-threshold operator, S
k-1represent the bent wave system number of the data reconstruction that kth-1 iteration obtains, meet
f
tand F
-1 trepresent the positive inverse-Fourier transform about time variable t,
represent the time domain data reconstruction that kth-1 iteration obtains, wherein
λ represents that N ties up threshold value set, λ={ λ
1, λ
2..., λ
n, and meet λ
1> λ
2> ... > λ
n, N represents maximum iteration time.
6. the three-dimensional irregular earthquake data re-establishing method of a kind of frequency field according to claim 3, it is characterized in that, threshold parameter formula, its expression formula is as follows:
Wherein Max is the maximal value of warp wavelet absolute coefficient; N is total iterations, and ε is the little value close to zero, relevant with the energy of noise in data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610008843.5A CN105425301B (en) | 2016-01-08 | 2016-01-08 | A kind of three-dimensional irregular earthquake data re-establishing method of frequency domain |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610008843.5A CN105425301B (en) | 2016-01-08 | 2016-01-08 | A kind of three-dimensional irregular earthquake data re-establishing method of frequency domain |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105425301A true CN105425301A (en) | 2016-03-23 |
CN105425301B CN105425301B (en) | 2018-10-19 |
Family
ID=55503625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610008843.5A Expired - Fee Related CN105425301B (en) | 2016-01-08 | 2016-01-08 | A kind of three-dimensional irregular earthquake data re-establishing method of frequency domain |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105425301B (en) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105974468A (en) * | 2016-05-04 | 2016-09-28 | 东华理工大学 | Method of simultaneously carrying out five-dimensional seismic data reconstruction and noise suppression |
CN106199701A (en) * | 2016-07-15 | 2016-12-07 | 中国石油天然气集团公司 | The method for reconstructing of irregular geological data and device |
CN106249291A (en) * | 2016-09-26 | 2016-12-21 | 东华理工大学 | A kind of high precision seismic data re-establishing method based on two-dimentional non-homogeneous warp wavelet |
CN106371138A (en) * | 2016-08-17 | 2017-02-01 | 中国石油天然气集团公司 | Seismic data reconstruction method and apparatus |
CN106646612A (en) * | 2016-12-20 | 2017-05-10 | 中国地质大学(北京) | Seismic data reconstruction method based on matrix reduced rank |
CN106970419A (en) * | 2017-04-10 | 2017-07-21 | 东华理工大学 | A kind of non-homogeneous bent ripple 3D seismic data method for reconstructing based on linear Bregman algorithms |
CN107121701A (en) * | 2017-05-05 | 2017-09-01 | 吉林大学 | The multi-component earthquake data Corssline directions wave field method for reconstructing converted based on Shearlet |
CN107272057A (en) * | 2017-05-16 | 2017-10-20 | 中国石油天然气股份有限公司 | Method and device for processing three-dimensional seismic data |
CN108387928A (en) * | 2018-02-11 | 2018-08-10 | 中国石油化工股份有限公司 | Data construction method based on seismic signature transformation space |
CN108519619A (en) * | 2018-04-11 | 2018-09-11 | 中国石油大学(华东) | Based on the time-frequency analysis technology for becoming privileged warp wavelet |
CN108919357A (en) * | 2018-05-16 | 2018-11-30 | 中国海洋石油集团有限公司 | A kind of ghost reflection drawing method based on frequency spectrum reconfiguration |
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 |
CN111325834A (en) * | 2020-03-24 | 2020-06-23 | 中海石油(中国)有限公司 | Modeling method and system based on digital image processing |
CN111337973A (en) * | 2018-12-19 | 2020-06-26 | 中国石油天然气集团有限公司 | Seismic data reconstruction method and system |
CN112241020A (en) * | 2019-07-16 | 2021-01-19 | 中国石油天然气集团有限公司 | Method and apparatus for determining undersampling rate in sparse seismic data acquisition |
CN112445649A (en) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | Seismic missing data recovery method, computer storage medium and computer system |
WO2021126360A1 (en) * | 2019-12-18 | 2021-06-24 | Bp Corporation North America Inc. | Enhanced projection on convex sets for interpolation and deblending |
CN114509805A (en) * | 2021-05-14 | 2022-05-17 | 中国地质大学(北京) | Vector convex set projection multi-component three-dimensional seismic data reconstruction method and device |
CN117148432A (en) * | 2023-10-27 | 2023-12-01 | 胜利信科(山东)勘察测绘有限公司 | Shallow profile data space interpolation method based on multi-scale component extraction |
CN118211067A (en) * | 2024-03-19 | 2024-06-18 | 哈尔滨工业大学 | Seismic data interpolation method based on self-supervision migration learning convex set projection network |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2460175A (en) * | 2008-05-21 | 2009-11-25 | Bp Corp North America Inc | Interpolating seismic data by projection onto convex sets |
CN104007469A (en) * | 2014-05-24 | 2014-08-27 | 长江大学 | Weak seismic signal reconstruction method based on curvelet transform |
CN104536044A (en) * | 2015-01-16 | 2015-04-22 | 中国石油大学(北京) | Interpolation and denoising method and system for seismic data |
-
2016
- 2016-01-08 CN CN201610008843.5A patent/CN105425301B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2460175A (en) * | 2008-05-21 | 2009-11-25 | Bp Corp North America Inc | Interpolating seismic data by projection onto convex sets |
CN104007469A (en) * | 2014-05-24 | 2014-08-27 | 长江大学 | Weak seismic signal reconstruction method based on curvelet transform |
CN104536044A (en) * | 2015-01-16 | 2015-04-22 | 中国石油大学(北京) | Interpolation and denoising method and system for seismic data |
Non-Patent Citations (5)
Title |
---|
HUA ZHANG,ET AL.: "3D seismic data reconstruction based on complex-valued curvelet transform in frequency domain", 《JOURNAL OF APPLIED GEOPHYSICS》 * |
WANG BEN-FENG,ET AL.: "CURVELET-BASED 3D RECONSTRUCTION OF DIGITAL CORES USING THE POCS METHOD", 《CHINESE JOURNAL OF GEOPHYSICS》 * |
ZHANG HUA,ET AL.: "Seismic data reconstruction based on CS and Fourier theory", 《APPLIED GEOPHYSICS》 * |
刘国昌等: "基于Curvelet变换的缺失地震数据插值方法", 《石油地球物理勘探》 * |
王本锋等: "POCS联合改进的Jitter采样理论曲波域地震数据重建", 《石油地球物理勘探》 * |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105974468A (en) * | 2016-05-04 | 2016-09-28 | 东华理工大学 | Method of simultaneously carrying out five-dimensional seismic data reconstruction and noise suppression |
CN105974468B (en) * | 2016-05-04 | 2018-10-19 | 东华理工大学 | A kind of method that can be carried out at the same time five dimension Reconstruction of seismic data and noise compacting |
CN106199701A (en) * | 2016-07-15 | 2016-12-07 | 中国石油天然气集团公司 | The method for reconstructing of irregular geological data and device |
CN106199701B (en) * | 2016-07-15 | 2018-06-01 | 中国石油天然气集团公司 | The method for reconstructing and device of irregular seismic data |
CN106371138B (en) * | 2016-08-17 | 2018-10-16 | 中国石油天然气集团公司 | Reconstruction of seismic data method and apparatus |
CN106371138A (en) * | 2016-08-17 | 2017-02-01 | 中国石油天然气集团公司 | Seismic data reconstruction method and apparatus |
CN106249291A (en) * | 2016-09-26 | 2016-12-21 | 东华理工大学 | A kind of high precision seismic data re-establishing method based on two-dimentional non-homogeneous warp wavelet |
CN106646612A (en) * | 2016-12-20 | 2017-05-10 | 中国地质大学(北京) | Seismic data reconstruction method based on matrix reduced rank |
CN106970419A (en) * | 2017-04-10 | 2017-07-21 | 东华理工大学 | A kind of non-homogeneous bent ripple 3D seismic data method for reconstructing based on linear Bregman algorithms |
CN106970419B (en) * | 2017-04-10 | 2019-01-25 | 东华理工大学 | A kind of non-homogeneous bent wave 3D seismic data method for reconstructing based on linear Bregman algorithm |
CN107121701A (en) * | 2017-05-05 | 2017-09-01 | 吉林大学 | The multi-component earthquake data Corssline directions wave field method for reconstructing converted based on Shearlet |
CN107272057A (en) * | 2017-05-16 | 2017-10-20 | 中国石油天然气股份有限公司 | Method and device for processing three-dimensional seismic data |
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 |
CN108519619A (en) * | 2018-04-11 | 2018-09-11 | 中国石油大学(华东) | Based on the time-frequency analysis technology for becoming privileged warp wavelet |
CN108919357B (en) * | 2018-05-16 | 2019-10-11 | 中国海洋石油集团有限公司 | A kind of ghost reflection drawing method based on frequency spectrum reconfiguration |
CN108919357A (en) * | 2018-05-16 | 2018-11-30 | 中国海洋石油集团有限公司 | A kind of ghost reflection drawing method based on frequency spectrum reconfiguration |
CN110764135A (en) * | 2018-07-26 | 2020-02-07 | 中国石油化工股份有限公司 | Irregular seismic data full-band reconstruction method |
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 |
CN111337973B (en) * | 2018-12-19 | 2022-08-05 | 中国石油天然气集团有限公司 | Seismic data reconstruction method and system |
CN111337973A (en) * | 2018-12-19 | 2020-06-26 | 中国石油天然气集团有限公司 | Seismic data reconstruction method and system |
CN112241020B (en) * | 2019-07-16 | 2024-04-30 | 中国石油天然气集团有限公司 | Method and device for determining undersampling rate in sparse seismic data acquisition |
CN112241020A (en) * | 2019-07-16 | 2021-01-19 | 中国石油天然气集团有限公司 | Method and apparatus for determining undersampling rate in sparse seismic data acquisition |
CN112445649A (en) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | Seismic missing data recovery method, computer storage medium and computer system |
WO2021126360A1 (en) * | 2019-12-18 | 2021-06-24 | Bp Corporation North America Inc. | Enhanced projection on convex sets for interpolation and deblending |
CN111325834B (en) * | 2020-03-24 | 2023-12-26 | 中海石油(中国)有限公司 | Modeling method and system based on digital image processing |
CN111325834A (en) * | 2020-03-24 | 2020-06-23 | 中海石油(中国)有限公司 | Modeling method and system based on digital image processing |
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 |
CN117148432A (en) * | 2023-10-27 | 2023-12-01 | 胜利信科(山东)勘察测绘有限公司 | Shallow profile data space interpolation method based on multi-scale component extraction |
CN117148432B (en) * | 2023-10-27 | 2024-03-19 | 胜利信科(山东)勘察测绘有限公司 | Shallow profile data space interpolation method based on multi-scale component extraction |
CN118211067A (en) * | 2024-03-19 | 2024-06-18 | 哈尔滨工业大学 | Seismic data interpolation method based on self-supervision migration learning convex set projection network |
Also Published As
Publication number | Publication date |
---|---|
CN105425301B (en) | 2018-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105425301A (en) | Frequency domain three-dimensional irregular earthquake data reconstruction method | |
CN107728211B (en) | Seismic signal recovery algorithm based on tensor nuclear norm regularization | |
CN106646612B (en) | Reconstruction of seismic data method based on matrix contraction | |
CN106970419B (en) | A kind of non-homogeneous bent wave 3D seismic data method for reconstructing based on linear Bregman algorithm | |
CN105974468B (en) | A kind of method that can be carried out at the same time five dimension Reconstruction of seismic data and noise compacting | |
CN107817527B (en) | Seismic exploration in desert stochastic noise suppression method based on the sparse compressed sensing of block | |
Naghizadeh et al. | Multidimensional de-aliased Cadzow reconstruction of seismic records | |
CN111158049B (en) | Seismic reverse time migration imaging method based on scattering integration method | |
CN103293551A (en) | Model constraint based impedance inversion method and model constraint based impedance inversion system | |
CN106842321B (en) | seismic data reconstruction method and device | |
CN107765308B (en) | Reconstruct low-frequency data frequency domain full waveform inversion method based on convolution thought Yu accurate focus | |
CN105467442A (en) | A globally optimized time-varying sparse deconvolution method and apparatus | |
CN105911584B (en) | Implicit staggered-grid finite difference elastic wave numerical simulation method and device | |
Zhang et al. | Seismic data reconstruction based on CS and Fourier theory | |
CA2800425A1 (en) | Surface-consistent amplitude and deconvolution simultaneous joined inversion | |
CN104730576A (en) | Curvelet transform-based denoising method of seismic signals | |
CN110208856A (en) | A kind of desert Complex Noise drawing method based on manifold subregion 2D-VMD | |
CN109738950A (en) | The noisy-type data primary wave inversion method of domain inverting is focused based on sparse 3 D | |
CN104391324A (en) | Seismic trace set dynamic correction stretching correction pre-processing technology before AVO inversion depending on frequency | |
CN113077386A (en) | Seismic data high-resolution processing method based on dictionary learning and sparse representation | |
CN105319593A (en) | Combined denoising method based on curvelet transform and singular value decomposition | |
CN105676292A (en) | 3D earthquake data de-noising method based on 2D curvelet transform | |
Wang et al. | Fast 3D time-domain airborne EM forward modeling using random under-sampling | |
CN109782346B (en) | Acquisition footprint pressing method based on morphological component analysis | |
Saengduean et al. | Multi-source wavefield reconstruction combining interferometry and compressive sensing: application to a linear receiver array |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181019 Termination date: 20200108 |