CN107678065A - 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
- CN107678065A CN107678065A CN201610624565.6A CN201610624565A CN107678065A CN 107678065 A CN107678065 A CN 107678065A CN 201610624565 A CN201610624565 A CN 201610624565A CN 107678065 A CN107678065 A CN 107678065A
- Authority
- CN
- China
- Prior art keywords
- reflectance factor
- factor matrix
- msub
- mrow
- matrix
- 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/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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses the guarantor for improving seismic resolution to construct well control space the Method of Deconvolution and device.This method is used along direction in space Regularization Technique on the basis of conventional deconvolution model, realizes multiple tracks deconvolution simultaneously；Well information has been incorporated into inversion result by initial model, used restraint and Quality Control, has made result more accurate, has avoided space wild effect, while improving longitudinal frame, 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 follow-up 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
Make well control space the Method of Deconvolution and a kind of guarantor for improving seismic resolution constructs well control space deconvolution device.
Background technology
Poststack seismic data can regard the convolution of wavelet and reflectance factor as, and which results in band while wavelet and data
Limit.This means 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 yardstick 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 not do the Method of Deconvolution of any hypothesis to the phase of seismic wavelet, such as
Homomorphic deconvolution, minimum entropy deconvolution etc..These common methods have its respective excellent, shortcoming.
Predictive deconvolution is applied in surface seismic exploration outside multiple suppression, is additionally operable to eliminate the sides such as marine singing effect
The implementation method and application oil instrument .2009 of the predictive deconvolution operator extraction of face (Hu Liu multitrace seismograms, 23 (5)：
55-57).It is the least-square filtering based on lienar prediction theory, i.e. Wiener filtering.Any signal all exists random dry in itself
Disturb, the purpose for doing signal transacting seeks to the determination composition that extraction is included in random signal.So 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 wave filter simultaneously
Out.It has two basic assumption conditions, first, seismic wavelet is minimum phase, second, 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, make in the form of the signal of convolution relation composition first is transformed to be multiplied, then the signal of multiplication relation composition is changed
For the form (multiple plug spectrum) of addition.It is that after being filtered computing in multiple plug spectral domain, original function is transformed into by inverse function
The broad sense in domain is, linear filtering system.But homomorphic deconvolution require wavelet based on low frequency, reflection coefficient sequence using high frequency as
It is main, therefore both can separate in intermediary heat spectrum (in fact both partially overlap, thus can not divide both completely
From), so as to obtain seismic wavelet.Wang Jun etc. (improvement of homomorphic deconvolution and application geophysical prospectings for oil .2003,38
Volume：Minimum entropy deconvolution MED, L mould and PARSIMONY 27-30) is used to be obtained for object function from superposed average logarithmic spectrum
Optimal wavelet and anti-wavelet are automatically determined among seismic wavelet and anti-wavelet, reduces human factor.Then seismic channel is mixed
Phase deconvolution, thus homomorphic deconvolution and minimum entropy deconvolution, L mould deconvolution are combined, mutually learnt from other's strong points to offset one's weaknesses.
Week it is superfine (the well control the Method of Deconvolution lithologic deposits .2010 under signal to noise spectrum constraint, 7：72-76) by drawing
Enter the high-frequency information of well-log information, conventional the Method of Deconvolution is improved, realize 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
Improve.It is that well information is added during deconvolution carries out Quality Control with the similarity of the present invention, still, when the letter of data
Make an uproar when including high wave number noise and error than relatively low or seismic data, this method can cause the reflectance factor of estimation or impedance horizontal
It is unstable.
Li Zhina etc. (number of patent application CN104914466A, 2015, September), which is proposed, a kind of improves seismic data resolution
Method.The technology is to be related to the pre-stack processing that more subwaves developmentally shake data.It is primarily based on SRME and focusing transform theory pushes away
Lead the multidimensional wavelet deconvolution model drawn using more subwaves；Then returned using the unstable state for introducing shaping regularization adaptive
Matched filtering method separates the high-resolution data recovered in focusing domain by more subwaves, realizes the high-resolution of initial data
Rate is changed.The technology, which is realized, is based on wave theory, from spatially using restraint, can improve while longitudinal frame is improved
The lateral resolution of data.It is with the similarity of the present invention during deconvolution from spatially using restraint, but two
Person based on principle and realization method and differ, what is finally wanted is a little that the patent does not add the information of well and is subject to
Constraint and Quality Control.
In a word, the Method of Deconvolution conventional at present application more adds a regularization term along time orientation, and passes through
Inverting mode obtains a stabilization and optimal solution.But these conventional methods are mostly using a lot of processing mode.When
The signal to noise ratio of data is relatively low or when seismic data includes high wave number noise and error, and single track regularization method often causes to estimate
Reflectance factor or impedance it is laterally unstable.In addition, conventional deconvolution processing have ignored the spatial relationship between Dao Yu roads, and
And the information for not adding well uses restraint and Quality Control.
The content of the invention
The present invention proposes a kind of guarantor for improving seismic resolution and constructs 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 include：
Step 101, the first reflectance factor matrix is obtained based on log；
Step 102, initial spatial regularization parameter sequence λ is providedi, i=1,2 ..., N, N be original multiple tracks earthquake
The total road number included in record, into step, 103；
Step 103, using the regularization along direction in space, based on multiple tracks described in current spatial regularization parameter sequence pair
Earthquake record carries out space deconvolution, the second reflectance factor matrix is obtained, into step 104；
Step 104, if the error between obtained the first reflectance factor matrix and the second reflectance factor matrix meet it is predetermined
It is expected, then into step 106, otherwise, into step 105；
Step 105, current spatial regularization parameter sequence is adjusted, is then return to step 103；
Step 106, using the second currently available reflectance factor matrix as final reflectance factor matrix, by each of which
The wavelet of reflection coefficient sequence corresponding to 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, disclose a kind of guarantor for improving seismic resolution and construct well control space deconvolution dress
Put, the device includes the first reflectance factor matrix determining unit, initialization unit, the second reflectance factor matrix determining unit, mistake
Poor judging unit, spatial regularization parameter adjustment unit, convolution superpositing unit, wherein：
First reflectance factor matrix determining unit, for obtaining the first reflectance factor matrix based on log；
Initialization unit, for providing initial spatial regularization parameter sequence λi, i=1,2 ..., N, wherein N are original
The total road number included in the multitrace seismogram of beginning；
Second reflectance factor matrix determining 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 the multitrace seismogram to obtain the second reflectance factor matrix, from first during beginning
Beginningization unit reception space regularization parameter sequence, the situation of preset expected is not met in error described in error judgment unit judges
Under, from spatial regularization parameter adjustment unit reception space regularization parameter sequence；
Error judgment unit, respectively from the first reflectance factor matrix determining unit and the second reflectance factor matrix determining unit
The first reflectance factor matrix and the second reflectance factor matrix are received, for judging the first reflectance factor matrix and the second reflectance factor
Whether the error between matrix meets preset expected；
Spatial regularization parameter adjustment unit, the situation of preset expected is not met in error described in error judgment unit judges
Under, for adjusting current spatial regularization parameter sequence；
Convolution superpositing unit, in the case where error described in error judgment unit judges meets preset expected, for by
The second currently available reflectance factor matrix of two reflectance factor matrix determining units as final reflectance factor matrix, and by its
In per the wavelet in reflection coefficient sequence and the road corresponding to together carry out convolution, and the convolution result of superposition multiple tracks, to obtain
Required post-stack seismic data body.
In the present invention, well control space deconvolution is constructed by protecting, caused by solving single-channel seismic data regularization method
Reflectance factor or the transversely unstable problem of wave impedance；The information of well is incorporated into inversion result by initial model simultaneously,
Use restraint and Quality Control, reach the purpose for improving seismic data resolution.The present invention can substantially increase earthquake deconvolution letter
Number processing accuracy, there is high industrial utility value and popularizing application prospect.
Brief description of the drawings
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 number
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 is shown realizes schematic diagram according to one concrete application example of 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 wells, B wells, the well control multichannel deconvolution result of C wells and the well filtering reflective obtained using the present invention
Index contrast figure.
Fig. 5 (a), (b) show certain original seismic profile and the deconvolution achievement Profile Correlation figure obtained using the present invention.
Embodiment
The preferred embodiment of the present invention is more fully described below with reference to accompanying drawings.Although the present invention is shown in accompanying 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, these embodiments are provided so that the present invention is more thorough and complete, and can be by the present invention's
Scope is intactly 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 include：
Step 101, the first reflectance factor matrix is obtained based on log, for example, can be bent in the well logging to collecting
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 be original multiple tracks earthquake
The total road number included in record, into step, 103；
Step 103, using the regularization along direction in space, based on multiple tracks described in current spatial regularization parameter sequence pair
Earthquake record carries out space deconvolution, the second reflectance factor matrix is obtained, into step 104；
Step 104, if the error between obtained the first reflectance factor matrix and the second reflectance factor matrix meet it is predetermined
It is expected, then into step 106, otherwise, into step 105；
Step 105, current spatial regularization parameter sequence is adjusted, is then return to step 103；
Step 106, using the second currently available reflectance factor matrix as final reflectance factor matrix, by each of which
The wavelet of reflection coefficient sequence corresponding to 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 beneficial effect：1. utilization space Regularization Technique realizes that multichannel seismic data is same
Shi Jinhang deconvolution seismic data processings, improve conventional method reflectance factor caused by single-channel seismic data regularization method
Or impedance transverse direction wild effect；2. enter row constraint and Quality Control using well logging information so that inversion result is more accurate, preferably protects
The border on stratum is protected, as edge of tomography, crack, pinching and sand body etc..In addition, it will be appreciated by those skilled in the art that
It is that step 101 necessarily performs before step 102~102, step 101 can be held in the random time before step 104
OK.
Alternatively, in step 103, obtaining the second reflectance factor matrix can include：
Sets target 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
Ripple convolution matrix, m represent the second reflectance factor matrix to be asked, λiFor spatial regularization parameter corresponding to i-th, DiFor i-th
The frequency wavenumber domain operator of corresponding representation space structure, | | * | |11 norm is represented, | | * | |22 norms are represented,
Actual seismic data and the residual energy of synthetic seismogram can be considered as,It can be considered as adding just
Then change the factor, to play a part of stable algorithm；
Iterative above formula obtains the matrix m for causing object function E (m) to be less than the first predetermined threshold value, and matrix m as exists
The the second reflectance factor matrix obtained under current spatial regularization parameter sequence.
Alternatively, at step 104, the first reflectance factor matrix and the second reflectance factor can be judged by following method
Whether the error between matrix meets preset expected：
Sets target function min (m)：
Wherein, m represents the second reflectance factor matrix, and S represents the first reflectance factor matrix, and C represents covariance matrix, λiFor
Spatial regularization parameter corresponding to i-th, DiFor the frequency wavenumber domain operator of representation space structure corresponding to i-th, | | * | |1
1 norm is represented, | | * | |pP norms are represented, | | C (m-S) | |pThe difference of theoretical value and well observation data calculated can be considered as,The regularization factors added can be considered as, to play a part of stable algorithm；
If min (m) is less than the second predetermined threshold value, judge between the first reflectance factor matrix and the second reflectance factor matrix
Error meet preset expected, conversely, then judging that the error between the first reflectance factor matrix and the second reflectance factor matrix is not inconsistent
Close preset expected.
Embodiment 2
The invention also discloses a kind of guarantor for improving seismic resolution to construct well control space deconvolution device.The present embodiment
In, the first reflectance factor matrix determining unit, initialization unit, the second reflectance factor matrix determining unit, error judgment unit,
Spatial regularization parameter adjustment unit, convolution superpositing unit, wherein：
First reflectance factor matrix determining unit, for obtaining the first reflectance factor matrix based on log；
Initialization unit, for providing initial spatial regularization parameter sequence λi, i=1,2 ..., N, wherein N are original
The total road number included in the multitrace seismogram of beginning；
Second reflectance factor matrix determining 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 the multitrace seismogram to obtain the second reflectance factor matrix, from first during beginning
Beginningization unit reception space regularization parameter sequence, the feelings of preset expected are not met in error described in error judgment unit judges
Under condition, from spatial regularization parameter adjustment unit reception space regularization parameter sequence；
Error judgment unit, respectively from the first reflectance factor matrix determining unit and the second reflectance factor matrix determining unit
The first reflectance factor matrix and the second reflectance factor matrix are received, for judging the first reflectance factor matrix and the second reflectance factor
Whether the error between matrix meets preset expected；
Spatial regularization parameter adjustment unit, the situation of preset expected is not met in error described in error judgment unit judges
Under, for adjusting current spatial regularization parameter sequence；
Convolution superpositing unit, in the case where error described in error judgment unit judges meets preset expected, for by
The second currently available reflectance factor matrix of two reflectance factor matrix determining units as final reflectance factor matrix, and by its
In per the wavelet in reflection coefficient sequence and the road corresponding to together carry out convolution, and the convolution result of superposition multiple tracks, to obtain
Required post-stack seismic data body.
Alternatively, in the second reflectance factor matrix determining unit, obtaining the second reflectance factor matrix can include：
Sets target 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
Ripple convolution matrix；M represents the second reflectance factor matrix to be asked, λiFor spatial regularization parameter corresponding to i-th, DiFor i-th
The frequency wavenumber domain operator of corresponding representation space structure, | | * | |11 norm is represented, | | * | |2Represent 2 norms；
Iterative above formula obtains the matrix m for causing object function E (m) to be less than the first predetermined threshold value, and matrix m as exists
The the second reflectance factor matrix obtained under current spatial regularization parameter sequence.
Alternatively, in error judgment unit, the first reflectance factor matrix and second anti-can be judged by following method
Whether the error penetrated between coefficient matrix meets preset expected：
Sets target function min (m)：
Wherein, m represents the second reflectance factor matrix, and S represents the first reflectance factor matrix, and C represents covariance matrix, λiFor
Spatial regularization parameter corresponding to i-th, DiFor the frequency wavenumber domain operator of representation space structure corresponding to i-th, | | * | |1
1 norm is represented, | | * | |pRepresent p norms；
If min (m) is less than the second predetermined threshold value, judge between the first reflectance factor matrix and the second reflectance factor matrix
Error meet preset expected, conversely, then judging that the error between the first reflectance factor matrix and the second reflectance factor matrix is not inconsistent
Close preset expected.
Using example
For ease of understanding the scheme of the embodiment of the present invention and its effect, a concrete application example given below.This area
It should be understood to the one skilled in the art that the example, only for the purposes of understanding the present invention, its any detail is not intended to be limited in any way
The system present invention.
Fig. 2 is shown realizes schematic diagram according to one concrete application example of 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 681 earthquake record altogether, and time sampling interval is 2ms, and the derived reflection of Amoco models is
Number is with after 30Hz dominant frequency Ricker wavelet convolutions, having added 20% random noise as original test data.Fig. 3 (b) is shown
Without regularization term deconvolution post-stack seismic data body.Fig. 3 (c) is the deconvolution post-stack seismic data body for applying the present 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
The deconvolution post-stack seismic data body for changing item is obvious unstable；And the deconvolution post-stack seismic data body of the present invention is stable
, and low-frequency information is preferably recovered.In addition, compared to firsthand information, the post-stack seismic data body obtained using the present invention
Longitudinal frame be improved, and horizontal information has more preferable continuity.
Fig. 4 shows A wells, B wells, the well control multichannel deconvolution of C wells obtained in another application example using the present invention
As a result with well filtering reflective index contrast figure.
D gas fields 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, and the dominant frequency of three dimensional seismic data is 25Hz~30Hz at present, and effective high frequency is also 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) it is 20 to select initial space regularization parameter, and multiple tracks original seismic data is carried out into space deconvolution, obtains ground
Shake reflectivity series；
(2) after being filtered processing to well curve, the reflectivity series of well are extracted；
(3) error for obtaining the reflectivity series of well and fractal series progress computing is as constraint；
(4) now the reflectivity series of well and fractal series error are larger；Then, spatial regularization is changed
Parameter is 25, repeats said process；
(5) after iteration, (now spatial regularization parameter is 30), output reflection coefficient when error is less than given threshold
Series, and it is subjected to convolution with the stable wavelet of original seismic data extraction 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
30) spatial regularization parameter determined eventually 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.Seismic profile reflectance signature after processing enriches, and meets
The feature that the sand body of geological knowledge mutually accelerates, well logging filtering postwave impedance contrast is bored with real, is 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 can 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.
It is described above various embodiments of the present invention, described above is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.In the case of without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes will be apparent from for the those of ordinary skill in art field.The selection of term used herein, purport
The principle of each embodiment, practical application or improvement to the technology in market are best being explained, or is making the art
Other those of ordinary skill are understood that 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 includes：
Step 101, the first reflectance factor matrix is obtained based on log；
Step 102, initial spatial regularization parameter sequence λ is providedi, i=1,2 ..., N, N be original multitrace seismogram
In total road number for including, into step 103；
Step 103, using the regularization along direction in space, based on original more described in current spatial regularization parameter sequence pair
Road earthquake record carries out space deconvolution, the second reflectance factor matrix is obtained, into step 104；
Step 104, if the error between obtained the first reflectance factor matrix and the second reflectance factor matrix meets preset expected,
Then enter step 106, otherwise, into step 105；
Step 105, current spatial regularization parameter sequence is adjusted, is then return to step 103；
Step 106, using the second currently available reflectance factor matrix as final reflectance factor matrix, by each of which road pair
The wavelet of the reflection coefficient sequence answered and the road carries out convolution, and is superimposed the convolution result of multiple tracks, with obtaining required poststack
Shake data volume.
2. according to the method for claim 1, wherein, in step 103, obtaining the second reflectance factor matrix includes：
Sets target function E (m)：
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<mo>|</mo>
<mi>d</mi>
<mo>-</mo>
<mi>G</mi>
<mi>m</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mi>m</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>,</mo>
</mrow>
Wherein, d is the original post-stack seismic data obtained after being overlapped to original multitrace seismogram
Body, G are wavelet convolution matrix；M represents the second reflectance factor matrix to be asked, λiJoin for spatial regularization corresponding to i-th
Number, DiFor the frequency wavenumber domain operator of representation space structure corresponding to i-th, | | * | |11 norm is represented, | | * | |2Represent 2 models
Number；
Iterative above formula obtains the matrix m for causing object function E (m) to be less than the first predetermined threshold value, and matrix m is current
Spatial regularization parameter sequence under obtained the second reflectance factor matrix.
3. according to the method for claim 1, wherein, at step 104, the first reflectance factor square is judged by following method
Whether the error between battle array and the second reflectance factor matrix meets preset expected：
Sets target function min (m)：
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msub>
<mo>|</mo>
<mi>p</mi>
</msub>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mi>m</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>,</mo>
</mrow>
Wherein, m represents the second reflectance factor matrix, and S represents the first reflectance factor matrix, and C represents covariance matrix, λiFor i-th
Corresponding spatial regularization parameter, DiFor the frequency wavenumber domain operator of representation space structure corresponding to i-th, | | * | |1Represent 1
Norm, | | * | |pRepresent p norms；
If min (m) is less than the second predetermined threshold value, the mistake between the first reflectance factor matrix and the second reflectance factor matrix is judged
Difference meets preset expected, conversely, the error for then judging between the first reflectance factor matrix and the second reflectance factor matrix do not meet it is pre-
Periodically hope.
4. a kind of guarantor for improving seismic resolution constructs well control space deconvolution device, the device includes the first reflectance factor matrix
Determining unit, initialization unit, the second reflectance factor matrix determining unit, error judgment unit, spatial regularization parameter adjustment
Unit, convolution superpositing unit, wherein：
First reflectance factor matrix determining unit, for obtaining the first reflectance factor matrix based on log；
Initialization unit, for providing initial spatial regularization parameter sequence λi, i=1,2 ..., N, wherein N are original more
The total road number included in road earthquake record；
Second reflectance factor matrix determining 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 reflectance factor matrix, during beginning from
Initialization unit reception space regularization parameter sequence, the feelings of preset expected are not met in error described in error judgment unit judges
Under condition, from spatial regularization parameter adjustment unit reception space regularization parameter sequence；
Error judgment unit, received respectively from the first reflectance factor matrix determining unit and the second reflectance factor matrix determining unit
First reflectance factor matrix and the second reflectance factor matrix, for judging the first reflectance factor matrix and the second reflectance factor matrix
Between error whether meet preset expected；
Spatial regularization parameter adjustment unit, in the case where 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 where error described in error judgment unit judges meets preset expected, for anti-by second
The second currently available reflectance factor matrix of coefficient matrix determining unit is penetrated as final reflectance factor matrix, and will be wherein every
The wavelet in reflection coefficient sequence corresponding to together and the road carries out convolution, and the convolution result of superposition multiple tracks, needed for obtaining
Post-stack seismic data body.
5. device according to claim 4, wherein, in the second reflectance factor matrix determining unit, obtain the second reflection
Coefficient matrix includes：
Sets target function E (m)：
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<mo>|</mo>
<mi>d</mi>
<mo>-</mo>
<mi>G</mi>
<mi>m</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mi>m</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>,</mo>
</mrow>
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 represents the second reflectance factor matrix to be asked, λiFor spatial regularization parameter corresponding to i-th, DiIt is corresponding for i-th
Representation space structure frequency wavenumber domain operator, | | * | |11 norm is represented, | | * | |2Represent 2 norms；
Iterative above formula obtains the matrix m for causing object function E (m) to be less than the first predetermined threshold value, and matrix m is current
Spatial regularization parameter sequence under obtained the second reflectance factor matrix.
6. device according to claim 4, wherein, in error judgment unit, the first reflection is judged by following method
Whether the error between coefficient matrix and the second reflectance factor matrix meets preset expected：
Sets target function min (m)：
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msub>
<mo>|</mo>
<mi>p</mi>
</msub>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mi>m</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>,</mo>
</mrow>
Wherein, m represents the second reflectance factor matrix, and S represents the first reflectance factor matrix, and C represents covariance matrix, λiFor i-th
Corresponding spatial regularization parameter, DiFor the frequency wavenumber domain operator of representation space structure corresponding to i-th, | | * | |1Represent 1
Norm, | | * | |pRepresent p norms；
If min (m) is less than the second predetermined threshold value, the mistake between the first reflectance factor matrix and the second reflectance factor matrix is judged
Difference meets preset expected, conversely, the error for then judging between the first reflectance factor matrix and the second reflectance factor 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 true CN107678065A (en) | 2018-02-09 |
| CN107678065B 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) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113534239A (en) * | 2020-04-18 | 2021-10-22 | 中国石油化工股份有限公司 | Waveform compression method based on signal separation |
| CN114910955A (en) * | 2022-05-10 | 2022-08-16 | 电子科技大学 | 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 |
|---|
| 王大兴 等: "分形反褶积方法研究及应用", 《石油物探》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113534239A (en) * | 2020-04-18 | 2021-10-22 | 中国石油化工股份有限公司 | Waveform compression method based on signal separation |
| CN114910955A (en) * | 2022-05-10 | 2022-08-16 | 电子科技大学 | Data-driven seismic deconvolution method based on error constraint and sparse representation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107678065B (en) | 2019-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3076205B1 (en) | Method for survey data processing compensating for visco-acoustic effects in tilted transverse isotropy reverse time migration | |
| US9470811B2 (en) | Creating a high resolution velocity model using seismic tomography and impedance inversion | |
| US6931324B2 (en) | Method for determining formation quality factor from seismic data | |
| CN104849756B (en) | A kind of seismic data resolution that improves strengthens the method for effective weak signal energy | |
| US7616524B1 (en) | Wavelet based intercept attribute for seismic exploration | |
| CN104237945B (en) | A kind of seismic data self adaptation high resolution processing method | |
| CN112882099B (en) | Earthquake frequency band widening method and device, medium and electronic equipment | |
| CN107179550B (en) | A kind of seismic signal zero phase deconvolution method of data-driven | |
| Khosro Anjom et al. | Full-waveform matching of VP and VS models from surface waves | |
| CN103364826A (en) | An earthquake blind source deconvolution method based on independent component analysis | |
| Li et al. | Seismic coherence for discontinuity interpretation | |
| Xue et al. | Q-factor estimation by compensation of amplitude spectra in synchrosqueezed wavelet domain | |
| CN107678065B (en) | The guarantor for improving seismic resolution constructs well control space the Method of Deconvolution and device | |
| CN113156500A (en) | Data-driven rapid construction constraint prestack seismic multi-channel inversion method | |
| Domenzain et al. | Joint full-waveform ground-penetrating radar and electrical resistivity inversion applied to field data acquired on the surface | |
| Mozayan et al. | Blocky inversion of multichannel elastic impedance for elastic parameters | |
| Lloyd | An investigation of the role of low frequencies in seismic impedance inversion | |
| Malviæ et al. | Cokriging geostatistical mapping and importance of quality of seismic attribute (s) | |
| Ha | Seismic data conditioning, attribute analysis, and machine-learning facies classification: Applications to Texas panhandle, Australia, New Zealand, and Gulf of Mexico | |
| Jiang et al. | Source-Independent Full-Waveform Inversion Based on Convolutional Wasserstein Distance Objective Function | |
| Sun et al. | Pure–P data inversion and its application in thin gas shale bed prediction | |
| Tian et al. | Upscaling method based on sedimentary cycle prior and its application to well‐to‐seismic tie | |
| De Abreu et al. | The use of spectral decomposition to quality control seismic data | |
| CN106125138A (en) | The method that frequency processes opened up by a kind of prestack CRP road collection | |
| Cheng et al. | Structural interval-Q estimation by using semi-blind time-variant sparse deconvolution |
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 |