CN107678065B - The guarantor for improving seismic resolution constructs well control space the Method of Deconvolution and device - Google Patents
The guarantor for improving seismic resolution constructs well control space the Method of Deconvolution and device Download PDFInfo
- Publication number
- CN107678065B CN107678065B CN201610624565.6A CN201610624565A CN107678065B CN 107678065 B CN107678065 B CN 107678065B CN 201610624565 A CN201610624565 A CN 201610624565A CN 107678065 B CN107678065 B CN 107678065B
- Authority
- CN
- China
- Prior art keywords
- reflection coefficient
- coefficient matrix
- matrix
- regularization parameter
- deconvolution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
Abstract
The invention discloses the guarantors for improving seismic resolution to construct well control space the Method of Deconvolution and device.This method, using along direction in space Regularization Technique, realizes multiple tracks deconvolution simultaneously on the basis of conventional deconvolution model；Well information has been introduced into inversion result by initial model, has been used restraint and Quality Control, has kept result more accurate, has avoided space wild effect, while improving longitudinal resolution, has kept the lateral continuity of high-resolution result.In practical application, this method can effectively improve the recognition accuracy to thin sand body, also have highly important practical significance to subsequent poststack seismic attributes analysis.
Description
Technical field
The present invention relates to oil gas technical field of physical geography, more particularly, to a kind of guarantor's structure for improving seismic resolution
It makes well control space the Method of Deconvolution and a kind of guarantor for improving seismic resolution constructs well control space deconvolution device.
Background technique
Poststack seismic data can regard the convolution of wavelet and reflection coefficient as, and which results in bands while wavelet and data
Limit.This means that wavelet has filtered the out-of-band important geological information of some seismic datas.High-Resolution Seismic Data knot
Fruit, for thin reservoir, the identification of the small scale geologic body such as construction, craven fault a little, straticulate exploration subdivision, oil sheet identification, well
Between connectivity analysis, and to establish fine three-dimensional oil reservoir geological model most important.And deconvolution is to improve seismic data point
Resolution and a kind of important means for expanding frequency band.Current the Method of Deconvolution can generally be divided into two classes, and one kind is to predict
Deconvolution is the Method of Deconvolution of representative；Another kind of is the Method of Deconvolution that any hypothesis is not done to the phase of seismic wavelet, such as
Homomorphic deconvolution, minimum entropy deconvolution etc..These common methods have its respective excellent, disadvantage.
Predictive deconvolution is applied in surface seismic exploration outside multiple suppression, is also used to eliminate the sides such as marine singing effect
The implementation method and application petroleum instrument .2009 of the predictive deconvolution operator extraction of face (Hu Liu multitrace seismogram, 23 (5):
55-57).It is the least-square filtering based on lienar prediction theory, i.e. Wiener filtering.Any signal itself all exists random dry
It disturbs, the purpose for doing signal processing seeks to the determination ingredient extracted include in random signal.In this way Wiener Filtering is exactly
A kind of linear filter, can be by output end signal performance reported as precisely as possible when signal and random noise input the filter simultaneously
Out.It has there are two basic assumption condition, first is that seismic wavelet is minimum phase, second is that reflection coefficient sequence is mistake at random
Journey.But actual seismic data is difficult to meet the two assumed conditions.
Homomorphic deconvolution is a kind of non-linear filtering method.It is theoretical using homomorphic filtering, by a kind of Function Mapping to another
One function field is made the signal formed with convolution relationship first be transformed to the form being multiplied, then the signal of multiplication relationship composition is converted
For the form (multiple plug spectrum) of addition.It is after being filtered operation in multiple plug spectral domain, to be transformed into original function by inverse function
The broad sense in domain is, linear filtering system.But homomorphic deconvolution requires wavelet based on low frequency, reflection coefficient sequence is with high frequency
It is main, therefore the two can separate to (in fact the two partially overlaps, thus cannot divide completely the two intermediary heat spectrum is upper
From), to obtain seismic wavelet.Wang Jun etc. (improvement of homomorphic deconvolution and application geophysical prospecting for oil .2003,38
Volume: 27-30) use minimum entropy deconvolution MED, L mould and PARSIMONY for objective function, it is obtained from superposed average logarithmic spectrum
Best wavelet and anti-wavelet are automatically determined in seismic wavelet and anti-wavelet, reduce human factor.Then seismic channel is mixed
Phase deconvolution thus combines homomorphic deconvolution and minimum entropy deconvolution, L mould deconvolution, mutually learns from other's strong points to offset one's weaknesses.
Week, superfine (well control the Method of Deconvolution lithologic deposit .2010,7:72-76 under signal to noise spectrum constraint) was by drawing
The high-frequency information for entering well-log information improves conventional the Method of Deconvolution, realizes well control deconvolution；Using signal to noise spectrum as matter
Constraint well control deconvolution parameter is chosen in amount control, it is ensured that deconvolution result keeps high s/n ratio while improving resolution ratio.The party
Method can make interval of interest seismic data frequency spectrum obtain certain expansion, and signal-to-noise ratio is held essentially constant in effective bandwidth, and slightly
It improves.It is that joined well information during deconvolution carries out Quality Control with similarity of the invention, still, when the letter of data
When relatively low or seismic data of making an uproar includes high wave number noise and error, this method will lead to the reflection coefficient or impedance transverse direction of estimation
It is unstable.
Li Zhina etc. (number of patent application CN104914466A, 2015, September) proposes a kind of raising seismic data resolution
Method.The technology is to be related to the pre-stack processing of multiple wave development seismic data.It is primarily based on SRME and focusing transform theory pushes away
Lead the multidimensional wavelet deconvolution model obtained using multiple wave；Then it is returned using the unstable state for introducing shaping regularization adaptive
Matched filtering method will be focused in domain and be separated by the high-resolution data that multiple wave is restored, and realize the high-resolution of initial data
Rate conversion.The technology, which is realized, can be improved while improving longitudinal resolution based on wave theory from spatially using restraint
The lateral resolution of data.It is during deconvolution with similarity of the invention from spatially using restraint, but two
Person based on principle and realization method it is not identical, what is finally wanted is a little that the patent is not added the information of well and is subject to
Constraint and Quality Control.
In short, conventional the Method of Deconvolution mostly application at present adds a regularization term along time orientation, and pass through
Inverting mode obtains a stabilization and optimal solution.But these conventional methods are mostly using a lot of processing mode.When
When the noise of data is relatively low or seismic data includes high wave number noise and error, single track regularization method often will lead to estimation
Reflection coefficient or impedance it is laterally unstable.In addition, conventional deconvolution processing has ignored the spatial relationship between the road Dao Yu, and
And the information that well is not added uses restraint and Quality Control.
Summary of the invention
The invention proposes a kind of guarantors for improving seismic resolution to construct well control space the Method of Deconvolution, and the present invention also proposes
Corresponding device.
According to an aspect of the invention, it is proposed that a kind of guarantor for improving seismic resolution constructs well control space deconvolution side
Method, this method comprises:
Step 101, the first reflection coefficient matrix is obtained based on log；
Step 102, initial spatial regularization parameter sequence λ is providedi, i=1,2 ..., N, N are original multiple tracks earthquake
The total road number for including in record, enters step, 103；
Step 103, using the regularization along direction in space, based on current spatial regularization parameter sequence to the multiple tracks
Earthquake record carries out space deconvolution, obtains the second reflection coefficient matrix, enters step 104；
Step 104, if the error between obtained the first reflection coefficient matrix and the second reflection coefficient matrix meet it is predetermined
It is expected that then entering step 106, otherwise, 105 are entered step；
Step 105, current spatial regularization parameter sequence is adjusted, step 103 is then return to；
It step 106, will be wherein each using the second currently available reflection coefficient matrix as final reflection coefficient matrix
The wavelet in the corresponding reflection coefficient sequence in road and the road carries out convolution, and is superimposed the convolution result of multiple tracks, to obtain required fold
Seismic data cube afterwards.
According to another aspect of the present invention, a kind of guarantor's construction well control space deconvolution dress for improving seismic resolution is disclosed
It sets, which includes the first reflection coefficient matrix determination unit, initialization unit, the second reflection coefficient matrix determination unit, misses
Poor judging unit, spatial regularization parameter adjustment unit, convolution superpositing unit, in which:
First reflection coefficient matrix determination unit, for obtaining the first reflection coefficient matrix based on log；
Initialization unit, for providing initial spatial regularization parameter sequence λi, i=1,2 ..., N, wherein N is original
The total road number for including in the multitrace seismogram of beginning；
Second reflection coefficient matrix determination unit, using the regularization along direction in space, for based on current space just
Then change argument sequence and space deconvolution is carried out to obtain the second reflection coefficient matrix to the multitrace seismogram, starts Shi Congchu
The case where beginningization unit reception space regularization parameter sequence, the error described in error judgment unit judges does not meet preset expected
Under, from spatial regularization parameter adjustment unit reception space regularization parameter sequence；
Error judgment unit, respectively from the first reflection coefficient matrix determination unit and the second reflection coefficient matrix determination unit
The first reflection coefficient matrix and the second reflection coefficient matrix are received, for judging the first reflection coefficient matrix and the second reflection coefficient
Whether the error between matrix meets preset expected；
The case where spatial regularization parameter adjustment unit, the error described in error judgment unit judges does not meet preset expected
Under, for adjusting current spatial regularization parameter sequence；
Convolution superpositing unit, in the case that the error described in error judgment unit judges meets preset expected, for by the
The second currently available reflection coefficient matrix of two reflection coefficient matrix determination units as final reflection coefficient matrix, and by its
In convolution, and the convolution result of superposition multiple tracks are carried out per the wavelet in corresponding reflection coefficient sequence and the road together, to obtain
Required post-stack seismic data body.
In the present invention, the deconvolution of well control space is constructed by protecting, is solved caused by single-channel seismic data regularization method
Reflection coefficient or the transversely unstable problem of wave impedance；The information of well is introduced into inversion result by initial model simultaneously,
It uses restraint and Quality Control, achievees the purpose that improve seismic data resolution.The present invention can be greatly improved earthquake deconvolution letter
Number processing accuracy, have high industrial utility value and popularization and application foreground.
Detailed description of the invention
Exemplary embodiment of the invention is described in more detail in conjunction with the accompanying drawings, it is of the invention above-mentioned and its
Its purpose, feature and advantage will be apparent, wherein in exemplary embodiment of the invention, identical reference label
Typically represent same parts.
Fig. 1 shows the flow chart according to an embodiment of the invention for protecting construction well control space the Method of Deconvolution.
Fig. 2 shows the realization schematic diagrames of a specific application example according to the present invention.
Fig. 3 (a), (b), (c) show certain original noisy seismic data cube figure, without regularization deconvolution post-stack seismic data
Body figure, the post-stack seismic data body figure obtained using the present invention.
Fig. 4 shows A well, B well, the well control multichannel deconvolution result of C well and the well filtering reflective obtained using the present invention
Index contrast figure.
The deconvolution achievement Profile Correlation figure that Fig. 5 (a), (b) show certain original seismic profile and obtain using the present invention.
Specific embodiment
The preferred embodiment that the present invention will be described in more detail below with reference to accompanying drawings.Although showing the present invention in attached drawing
Preferred embodiment, however, it is to be appreciated that may be realized in various forms the present invention without the embodiment party that should be illustrated here
Formula is limited.On the contrary, thesing embodiments are provided so that the present invention is more thorough and complete, and can will be of the invention
Range is completely communicated to those skilled in the art.
Embodiment 1
Fig. 1 shows the flow chart according to an embodiment of the invention for protecting construction well control space the Method of Deconvolution.According to
The present embodiment, this method comprises:
Step 101, the first reflection coefficient matrix is obtained based on log, for example, can be bent to collected well logging
Line is filtered the reflection coefficient sequence of extraction well after seismic data processing operation；
Step 102, initial spatial regularization parameter sequence λ is providedi, i=1,2 ..., N, N are original multiple tracks earthquake
The total road number for including in record, enters step, 103；
Step 103, using the regularization along direction in space, based on current spatial regularization parameter sequence to the multiple tracks
Earthquake record carries out space deconvolution, obtains the second reflection coefficient matrix, enters step 104；
Step 104, if the error between obtained the first reflection coefficient matrix and the second reflection coefficient matrix meet it is predetermined
It is expected that then entering step 106, otherwise, 105 are entered step；
Step 105, current spatial regularization parameter sequence is adjusted, step 103 is then return to；
It step 106, will be wherein each using the second currently available reflection coefficient matrix as final reflection coefficient matrix
The wavelet in the corresponding reflection coefficient sequence in road and the road carries out convolution, and is superimposed the convolution result of multiple tracks, to obtain required fold
Seismic data cube afterwards.
Above-described embodiment at least has following the utility model has the advantages that 1. realize that multichannel seismic data is same using spatial regularization technology
Shi Jinhang deconvolution seismic data processing improves conventional method because of reflection coefficient caused by single-channel seismic data regularization method
Or impedance transverse direction wild effect；2. carrying out constraint and Quality Control using well logging information preferably to protect so that inversion result is more acurrate
The boundary for having protected stratum, as the edge etc. of tomography, crack, pinching and sand body.In addition, it will be appreciated by those skilled in the art that
It is that step 101 necessarily executes before step 102~102, step 101 can be held in any time before step 104
Row.
Optionally, in step 103, obtaining the second reflection coefficient matrix may include:
It sets objective function E (m):
Wherein, d is the original post-stack seismic data body obtained after being overlapped to original multitrace seismogram, and G is son
Wave convolution matrix, m indicate the second reflection coefficient matrix to be asked, λiFor the i-th corresponding spatial regularization parameter, DiIt is i-th
The frequency wavenumber domain operator of corresponding representation space structure, | | * | |1Indicate 1 norm, | | * | |2Indicate 2 norms,
It can be considered as the residual energy of actual seismic data and synthetic seismogram,It can be considered as being added just
Change the factor, then to play the role of stable algorithm；
It iteratively solves above formula to obtain so that objective function E (m) is less than the matrix m of the first preset threshold, matrix m is to exist
The the second reflection coefficient matrix obtained under current spatial regularization parameter sequence.
Optionally, at step 104, the first reflection coefficient matrix and the second reflection coefficient can be judged by following method
Whether the error between matrix meets preset expected:
It sets objective function min (m):
Wherein, m indicates that the second reflection coefficient matrix, S indicate that the first reflection coefficient matrix, C indicate covariance matrix, λiFor
I-th corresponding spatial regularization parameter, DiFor the frequency wavenumber domain operator of the i-th corresponding representation space structure, | | * | |1
Indicate 1 norm, | | * | |pIndicate p norm, | | C (m-S) | |pIt can be considered as the difference of calculated theoretical value and well observation data,It can be considered as the regularization factors being added, to play the role of stable algorithm；
If min (m) judges between the first reflection coefficient matrix and the second reflection coefficient matrix less than the second preset threshold
Error meet preset expected, conversely, then judging that the error between the first reflection coefficient matrix and the second reflection coefficient matrix is not inconsistent
Close preset expected.
Embodiment 2
The invention also discloses a kind of guarantors for improving seismic resolution to construct well control space deconvolution device.The present embodiment
In, the first reflection coefficient matrix determination unit, initialization unit, the second reflection coefficient matrix determination unit, error judgment unit,
Spatial regularization parameter adjustment unit, convolution superpositing unit, in which:
First reflection coefficient matrix determination unit, for obtaining the first reflection coefficient matrix based on log；
Initialization unit, for providing initial spatial regularization parameter sequence λi, i=1,2 ..., N, wherein N is original
The total road number for including in the multitrace seismogram of beginning；
Second reflection coefficient matrix determination unit, using the regularization along direction in space, for based on current space just
Then change argument sequence and space deconvolution is carried out to obtain the second reflection coefficient matrix to the multitrace seismogram, starts Shi Congchu
Beginningization unit reception space regularization parameter sequence, the error described in error judgment unit judges do not meet the feelings of preset expected
Under condition, from spatial regularization parameter adjustment unit reception space regularization parameter sequence；
Error judgment unit, respectively from the first reflection coefficient matrix determination unit and the second reflection coefficient matrix determination unit
The first reflection coefficient matrix and the second reflection coefficient matrix are received, for judging the first reflection coefficient matrix and the second reflection coefficient
Whether the error between matrix meets preset expected；
The case where spatial regularization parameter adjustment unit, the error described in error judgment unit judges does not meet preset expected
Under, for adjusting current spatial regularization parameter sequence；
Convolution superpositing unit, in the case that the error described in error judgment unit judges meets preset expected, for by the
The second currently available reflection coefficient matrix of two reflection coefficient matrix determination units as final reflection coefficient matrix, and by its
In convolution, and the convolution result of superposition multiple tracks are carried out per the wavelet in corresponding reflection coefficient sequence and the road together, to obtain
Required post-stack seismic data body.
Optionally, in the second reflection coefficient matrix determination unit, obtaining the second reflection coefficient matrix may include:
It sets objective function E (m):
Wherein, d is the original post-stack seismic data body obtained after being overlapped to original multitrace seismogram, and G is son
Wave convolution matrix；M indicates the second reflection coefficient matrix to be asked, λiFor the i-th corresponding spatial regularization parameter, DiIt is i-th
The frequency wavenumber domain operator of corresponding representation space structure, | | * | |1Indicate 1 norm, | | * | |2Indicate 2 norms；
It iteratively solves above formula to obtain so that objective function E (m) is less than the matrix m of the first preset threshold, matrix m is to exist
The the second reflection coefficient matrix obtained under current spatial regularization parameter sequence.
Optionally, in error judgment unit, the first reflection coefficient matrix and second can be judged by following method instead
Whether the error penetrated between coefficient matrix meets preset expected:
It sets objective function min (m):
Wherein, m indicates that the second reflection coefficient matrix, S indicate that the first reflection coefficient matrix, C indicate covariance matrix, λiFor
I-th corresponding spatial regularization parameter, DiFor the frequency wavenumber domain operator of the i-th corresponding representation space structure, | | * | |1
Indicate 1 norm, | | * | |pIndicate p norm；
If min (m) judges between the first reflection coefficient matrix and the second reflection coefficient matrix less than the second preset threshold
Error meet preset expected, conversely, then judging that the error between the first reflection coefficient matrix and the second reflection coefficient matrix is not inconsistent
Close preset expected.
Using example
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
Fig. 2 shows the realization schematic diagrames of a specific application example according to the present invention.
Fig. 3 (a), (b), (c) show certain original noisy seismic data cube figure, without regularization deconvolution post-stack seismic data
Body figure, the post-stack seismic data body figure obtained using the present invention.
Fig. 3 (a) shows total 681 earthquake record, and time sampling interval is 2ms, the derived reflection system of Amoco model
After several and 30Hz dominant frequency Ricker wavelet convolution, add 20% random noise as original test data.Fig. 3 (b) is shown
Without regularization term deconvolution post-stack seismic data body.Fig. 3 (c) be using deconvolution post-stack seismic data body of the invention, wherein
Initially given spatial regularization parameter is 0.1, the information for being constrained to true model 0-5Hz of well.As seen from the figure, no canonical
It is obvious unstable for changing the deconvolution post-stack seismic data body of item；And deconvolution post-stack seismic data body of the invention is stable
, and low-frequency information is preferably restored.In addition, compared to firsthand information, using obtained post-stack seismic data body of the invention
Longitudinal resolution be improved, and lateral information has better continuity.
Fig. 4 shows the well control multichannel deconvolution of the A well obtained in another application example using the present invention, B well, C well
As a result with well filtering reflective index contrast figure.
The gas field D earth's surface is desert and semidesert environment, and the frequency and signal-to-noise ratio of seismic data are relatively low, point of seismic data
Resolution is also very low, at present the dominant frequency of three dimensional seismic data be 25Hz~30Hz, effective high frequency also be 50Hz, 10 meters or so or small
It is extremely difficult in the identification of 10 meters of sand bodies, therefore it is essential to improve seismic data resolution processing.Using the specific mistake of the present invention
Journey is as follows:
(1) selecting initial space regularization parameter is 20, and multiple tracks original seismic data is carried out space deconvolution, obtains ground
Shake reflectivity series；
(2) after being filtered to well curve, the reflectivity series of well are extracted；
(3) error for obtaining the reflectivity series of well and the progress operation of fractal series is as constraint；
(4) reflectivity series of well and fractal series error are larger at this time；Then, spatial regularization is modified
Parameter is 25, is repeated the above process；
(5) after iteration, (spatial regularization parameter is 30 at this time), output reflection coefficient when error is less than given threshold
Series, and it is subjected to convolution with the stabilization wavelet that original seismic data is extracted and is superimposed, obtain final high-resolution seismic exploration
Data volume.
Fig. 5 (a) is an original seismic profile in certain work area, and Fig. 5 (b) is to apply obtained seismic profile of the invention (most
Eventually 30) determining spatial regularization parameter is.It can be found that well control multichannel deconvolution method maintains the transverse direction of seismic response
Stability, continuity, and effectively increase the resolution ratio of seismic data.Treated, and seismic profile reflectance signature is abundant, meets
The feature that the sand body of geological knowledge mutually becomes faster bores well logging filtering postwave impedance contrast with real, coincide after deconvolution with real result of boring,
And it can identify that most thin reservoir is 8 meters of thick reservoirs.
The present invention can be system, method and/or computer program product.Computer program product may include computer
Readable storage medium storing program for executing, containing for making processor realize the computer-readable program instructions of various aspects of the invention.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.The selection of term used herein, purport
In the principle, practical application or improvement to the technology in market for best explaining each embodiment, or make the art
Other those of ordinary skill can understand each embodiment disclosed herein.
Claims (6)
1. a kind of guarantor for improving seismic resolution constructs well control space the Method of Deconvolution, this method comprises:
Step 101, the first reflection coefficient matrix is obtained based on log；
Step 102, initial spatial regularization parameter sequence λ is providedi, i=1,2 ..., N, N are original multitrace seismogram
In include total road number, enter step 103；
Step 103, using the regularization along direction in space, based on current spatial regularization parameter sequence to described original more
Road earthquake record carries out space deconvolution, obtains the second reflection coefficient matrix, enters step 104；
Step 104, if the error between obtained the first reflection coefficient matrix and the second reflection coefficient matrix meets preset expected,
106 are then entered step, otherwise, enters step 105；
Step 105, current spatial regularization parameter sequence is adjusted, step 103 is then return to；
It step 106, will be wherein per right together using the second currently available reflection coefficient matrix as final reflection coefficient matrix
The wavelet of the reflection coefficient sequence and the road answered carries out convolution, and is superimposed the convolution result of multiple tracks, with obtaining required poststack
Shake data volume.
2. according to the method described in claim 1, wherein, in step 103, obtaining the second reflection coefficient matrix includes:
It sets objective function E (m):
Wherein, d is the original post-stack seismic data body obtained after being overlapped to original multitrace seismogram, and G is wavelet pleat
Product matrix；M indicates the second reflection coefficient matrix to be asked, λiFor the i-th corresponding spatial regularization parameter, DiIt is corresponding for i-th
Representation space structure frequency wavenumber domain operator, | | * | |1Indicate 1 norm, | | * | |2Indicate 2 norms；
Iterative solution above formula obtains so that objective function E (m) is less than the matrix m of the first preset threshold, and matrix m is current
Spatial regularization parameter sequence under obtained the second reflection coefficient matrix.
3. according to the method described in claim 1, wherein, at step 104, judging the first reflection coefficient square by following method
Whether the error between battle array and the second reflection coefficient matrix meets preset expected:
It sets objective function min (m):
Wherein, m indicates that the second reflection coefficient matrix, S indicate that the first reflection coefficient matrix, C indicate covariance matrix, λiIt is i-th
Corresponding spatial regularization parameter, DiFor the frequency wavenumber domain operator of the i-th corresponding representation space structure, | | * | |1Indicate 1
Norm, | | * | |pIndicate p norm；
If min (m) judges the mistake between the first reflection coefficient matrix and the second reflection coefficient matrix less than the second preset threshold
Difference meets preset expected, conversely, the error for then judging between the first reflection coefficient matrix and the second reflection coefficient matrix do not meet it is pre-
Periodically hope.
4. a kind of guarantor for improving seismic resolution constructs well control space deconvolution device, which includes the first reflection coefficient matrix
Determination unit, initialization unit, the second reflection coefficient matrix determination unit, error judgment unit, spatial regularization parameter adjustment
Unit, convolution superpositing unit, in which:
First reflection coefficient matrix determination unit, for obtaining the first reflection coefficient matrix based on log；
Initialization unit, for providing initial spatial regularization parameter sequence λi, i=1,2 ..., N, wherein N is original more
The total road number for including in road earthquake record；
Second reflection coefficient matrix determination unit, using the regularization along direction in space, for based on current spatial regularization
Argument sequence carries out space deconvolution to the original multitrace seismogram to obtain the second reflection coefficient matrix, when beginning from
Initialization unit reception space regularization parameter sequence, the case where error judgment unit judges error does not meet preset expected
Under, from spatial regularization parameter adjustment unit reception space regularization parameter sequence；
Error judgment unit is received from the first reflection coefficient matrix determination unit and the second reflection coefficient matrix determination unit respectively
First reflection coefficient matrix and the second reflection coefficient matrix, for judging the first reflection coefficient matrix and the second reflection coefficient matrix
Between error whether meet preset expected；
Spatial regularization parameter adjustment unit, in the case that the error described in error judgment unit judges does not meet preset expected,
For adjusting current spatial regularization parameter sequence；
Convolution superpositing unit, in the case that the error described in error judgment unit judges meets preset expected, for anti-by second
The second currently available reflection coefficient matrix of coefficient matrix determination unit is penetrated as final reflection coefficient matrix, and will be wherein every
The wavelet in corresponding reflection coefficient sequence and the road carries out convolution, and the convolution result of superposition multiple tracks together, needed for obtaining
Post-stack seismic data body.
5. device according to claim 4, wherein in the second reflection coefficient matrix determination unit, obtain the second reflection
Coefficient matrix includes:
It sets objective function E (m):
Wherein, d is the original post-stack seismic data body obtained after being overlapped to original multitrace seismogram, and G is wavelet pleat
Product matrix；M indicates the second reflection coefficient matrix to be asked, λiFor the i-th corresponding spatial regularization parameter, DiIt is corresponding for i-th
Representation space structure frequency wavenumber domain operator, | | * | |1Indicate 1 norm, | | * | |2Indicate 2 norms；
Iterative solution above formula obtains so that objective function E (m) is less than the matrix m of the first preset threshold, and matrix m is current
Spatial regularization parameter sequence under obtained the second reflection coefficient matrix.
6. device according to claim 4, wherein in error judgment unit, judge the first reflection by following method
Whether the error between coefficient matrix and the second reflection coefficient matrix meets preset expected:
It sets objective function min (m):
Wherein, m indicates that the second reflection coefficient matrix, S indicate that the first reflection coefficient matrix, C indicate covariance matrix, λiIt is i-th
Corresponding spatial regularization parameter, DiFor the frequency wavenumber domain operator of the i-th corresponding representation space structure, | | * | |1Indicate 1
Norm, | | * | |pIndicate p norm；
If min (m) judges the mistake between the first reflection coefficient matrix and the second reflection coefficient matrix less than the second preset threshold
Difference meets preset expected, conversely, the error for then judging between the first reflection coefficient matrix and the second reflection coefficient matrix do not meet it is pre-
Periodically hope.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610624565.6A CN107678065B (en) | 2016-08-02 | 2016-08-02 | The guarantor for improving seismic resolution constructs well control space the Method of Deconvolution and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610624565.6A CN107678065B (en) | 2016-08-02 | 2016-08-02 | The guarantor for improving seismic resolution constructs well control space the Method of Deconvolution and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107678065A CN107678065A (en) | 2018-02-09 |
| CN107678065B true CN107678065B (en) | 2019-07-30 |
Family
ID=61133419
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610624565.6A Active CN107678065B (en) | 2016-08-02 | 2016-08-02 | The guarantor for improving seismic resolution constructs well control space the Method of Deconvolution and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107678065B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113534239A (en) * | 2020-04-18 | 2021-10-22 | 中国石油化工股份有限公司 | Waveform compression method based on signal separation |
| CN114910955B (en) * | 2022-05-10 | 2023-06-23 | 电子科技大学 | Data-driven seismic deconvolution method based on error constraint and sparse representation |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5173880A (en) * | 1989-12-26 | 1992-12-22 | Exxon Production Research Company | Method of generating seismic wavelets using seismic range equation |
| CN102937720A (en) * | 2011-08-15 | 2013-02-20 | 中国石油化工股份有限公司 | Method for improving seismic data resolution ratio for well control |
| CN103389513A (en) * | 2013-08-02 | 2013-11-13 | 中国石油集团川庆钻探工程有限公司 | Method for increasing seismic data resolution by applying sonic logging information constraint inversion |
| CN103713317A (en) * | 2012-09-29 | 2014-04-09 | 北京诺克斯达石油科技有限公司 | Time-varying wavelet based deconvolution method for frequency division processing |
| CN103954992A (en) * | 2014-04-03 | 2014-07-30 | 中国石油天然气股份有限公司 | Deconvolution method and device |
| CN105445800A (en) * | 2015-11-13 | 2016-03-30 | 中国石油化工股份有限公司 | Thick sand body top differentiation lithologic reservoir identification method |
-
2016
- 2016-08-02 CN CN201610624565.6A patent/CN107678065B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5173880A (en) * | 1989-12-26 | 1992-12-22 | Exxon Production Research Company | Method of generating seismic wavelets using seismic range equation |
| CN102937720A (en) * | 2011-08-15 | 2013-02-20 | 中国石油化工股份有限公司 | Method for improving seismic data resolution ratio for well control |
| CN103713317A (en) * | 2012-09-29 | 2014-04-09 | 北京诺克斯达石油科技有限公司 | Time-varying wavelet based deconvolution method for frequency division processing |
| CN103389513A (en) * | 2013-08-02 | 2013-11-13 | 中国石油集团川庆钻探工程有限公司 | Method for increasing seismic data resolution by applying sonic logging information constraint inversion |
| CN103954992A (en) * | 2014-04-03 | 2014-07-30 | 中国石油天然气股份有限公司 | Deconvolution method and device |
| CN105445800A (en) * | 2015-11-13 | 2016-03-30 | 中国石油化工股份有限公司 | Thick sand body top differentiation lithologic reservoir identification method |
Non-Patent Citations (1)
| Title |
|---|
| 分形反褶积方法研究及应用;王大兴 等;《石油物探》;19961225;第35卷(第4期);第1-10页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107678065A (en) | 2018-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xue et al. | Amplitude-preserving iterative deblending of simultaneous source seismic data using high-order radon transform | |
| Lu et al. | Seismic spectral decomposition using deconvolutive short-time Fourier transform spectrogram | |
| EP3076205B1 (en) | Method for survey data processing compensating for visco-acoustic effects in tilted transverse isotropy reverse time migration | |
| US7616524B1 (en) | Wavelet based intercept attribute for seismic exploration | |
| CN104237945B (en) | A kind of seismic data self adaptation high resolution processing method | |
| GB2414802A (en) | Seismic event correlation and Vp-Vs estimation | |
| CN107179550B (en) | A kind of seismic signal zero phase deconvolution method of data-driven | |
| WO2015014762A2 (en) | Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media | |
| CN107102356A (en) | Seismic signal high resolution data processing methods based on CEEMD | |
| CN107255831A (en) | A kind of extracting method of prestack frequency dispersion attribute | |
| CN107356965B (en) | Reflection coefficient inverting method for predicting reservoir based on weighted superposition Noise Elimination strategy | |
| CN112882099B (en) | Earthquake frequency band widening method and device, medium and electronic equipment | |
| Li et al. | Seismic coherence for discontinuity interpretation | |
| CN107678065B (en) | The guarantor for improving seismic resolution constructs well control space the Method of Deconvolution and device | |
| CN109358365A (en) | A kind of thin layer recognition methods based on well shake information matches | |
| Lloyd | An investigation of the role of low frequencies in seismic impedance inversion | |
| CN106125138B (en) | A kind of method that prestack CRP trace gathers open up frequency processing | |
| Marfurt et al. | Comparison of 3-d edge detection seismic attributes to Vinton Dome Louisiana | |
| Ha | Seismic data conditioning, attribute analysis, and machine-learning facies classification: Applications to Texas panhandle, Australia, New Zealand, and Gulf of Mexico | |
| Malviæ et al. | Cokriging geostatistical mapping and importance of quality of seismic attribute (s) | |
| Lyu | Full Wave-equation Based Passive Seismic Imaging and Multispectral Seismic Geometric Attributes | |
| Sanchis | Signal processing techniques for the enhancement of marine seismic data | |
| Lv et al. | Integrated characterization of deep karsted carbonates in the Tahe Oilfield, Tarim Basin | |
| Donoso Linares | Sparse-Layer Inversion as a Tool to Delineate Thin Sandstones, Lower Magdalena Valley Region, Colombia: A Case Study | |
| De Abreu et al. | The use of spectral decomposition to quality control seismic data |
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 |