CN112558161A - Well constraint earthquake low-frequency recovery method based on compressed sensing - Google Patents
Well constraint earthquake low-frequency recovery method based on compressed sensing Download PDFInfo
- Publication number
- CN112558161A CN112558161A CN202011546748.3A CN202011546748A CN112558161A CN 112558161 A CN112558161 A CN 112558161A CN 202011546748 A CN202011546748 A CN 202011546748A CN 112558161 A CN112558161 A CN 112558161A
- Authority
- CN
- China
- Prior art keywords
- well
- frequency
- low
- reflection coefficient
- compressed sensing
- 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
Images
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
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Abstract
The invention discloses a well-constrained earthquake low-frequency recovery method based on compressed sensing, which relates to the technical field of oil geophysical exploration and comprises the following steps: a1, well data processing: estimating reflection coefficients from well data
Description
Technical Field
The invention relates to the technical field of geophysical exploration of petroleum, in particular to a well constraint earthquake low-frequency recovery method based on compressed sensing.
Background
The low-frequency information can reflect the basic trend of a stratum and improve the speed analysis precision and the imaging precision of a deep structure, but the low-frequency information of conventionally acquired seismic data can be seriously polluted by noise, the sparse characteristic of signals utilized by a compressive sensing theory maps unknown complete seismic data into known low-dimensional observation data through a measurement matrix, and the complete seismic data is recovered by solving an underdetermined equation set.
Disclosure of Invention
Technical problem to be solved
Aiming at the defects of the prior art, the invention provides a well-constrained earthquake low-frequency recovery method based on compressed sensing, and solves the problem that the prior method is greatly influenced by noise and has poor low-frequency information recovery effect under the condition of heavy noise interference.
(II) technical scheme
In order to achieve the above purposes, the technical scheme adopted by the invention is as follows: a well constraint earthquake low-frequency recovery method based on compressed sensing comprises the following steps:
a2, extracting global static seismic wavelets and constructing a measurement matrix equation WF;
a3, initialization parameters: giving an initial reflection coefficient r0, the maximum iteration number Lin of an internal loop, the maximum iteration number Lout of an external loop and a residual error epsilon arbitrarily, wherein the initial iteration number k is 1;
a4、r0r0, outer loop: judge G-WFrk||2>εand k≤Lout,r0=rk-1And internal circulation: j is 0to Lin-1, and the reflection coefficient r is r at the time of terminationk,
a5, obtaining a spectrum R after frequency extension: r ═ Fr;
a6, obtaining a recovered spectrum G: g ═ L (r) + H (G), L, H are operators for extracting low-frequency and high-frequency spectral values, respectively;
a7, seismic data after low frequency recovery: g ═ F-1(G)。
Preferably, said T isλAs a function of the threshold value, said (WF)TIs the conjugate transpose of WF, the WF is a measurement matrix composed of wavelet operator and Fourier transform operator, the measurement matrix is composed of wavelet operator and Fourier transform operatorIs RwThe conjugate transpose of RwExtracting position information of corresponding well for sampling operator, wherein alpha value is larger than (WF)TThe maximum characteristic value of WF, and the lambda is the sparsity of the adjusting reflection coefficient.
Preferably, the inner loop and the outer loop perform: k is equal to k +1 and k is equal to k +1,the reflection coefficient r ═ rkAnd G is Fourier transform of the noisy data.
(III) advantageous effects
The invention has the beneficial effects that:
the well-constrained earthquake low-frequency recovery method based on compressed sensing comprises the steps of carrying out sparse inversion on noise-containing data through a certain noise-containing earthquake record with the signal-to-noise ratio SNR of 3.14, as shown in figure 1 and figure 2, respectively adopting a classical compressed sensing low-frequency compensation method and a classical compressed sensing low-frequency compensation method to obtain a reflection coefficient graph 3 and figure 4 and a frequency spectrum graph 5 and figure 6 of the earthquake record after frequency broadening, further showing that the precision of the reflection coefficient based on compressed sensing theory inversion is low, partial noise remains exist, the reflection coefficient inverted by the method has high precision and does not contain the noise residue, comparing the frequency spectrum graphs before and after frequency broadening with the reflection coefficient graphs 5 and figure 6, the low-frequency information compensation effect based on compressed sensing theory compensation is poor, particularly the arrow head of the 4-10Hz low-frequency information graph 4 is not compensated, and the method is not influenced by noise ratio reduction, so that the low-frequency information is recovered, the frequency band is widened to a certain extent, namely the inversion process of well logging data constraint compressed sensing is adopted in the method, and the influence of strong noise on the inversion process is overcome, so that the problem that the low-frequency information recovery effect is poor under the condition of heavy noise interference due to large noise influence in the conventional method is solved.
Drawings
FIG. 1 is a diagram of the low frequency compensation effect of the present invention when SNR is 3.14 compared with noisy seismic data and its spectrum a noisy data;
FIG. 2 is a diagram of the low frequency compensation effect of the present invention when SNR is 3.14 compared with noisy seismic data and its spectrum b noisy data;
fig. 3 is a schematic diagram of the comparison of the low frequency compensation effect of the present invention when SNR is 3.14 with the reflection coefficient based on compressed sensing estimation and the recovered spectrum c noisy data;
FIG. 4 is a diagram of a reflection coefficient spectrum based on compressed sensing estimation according to the present invention;
FIG. 5 is a graph illustrating the estimated reflection coefficient of the present invention;
FIG. 6 is a diagram illustrating the estimated reflection coefficient and the recovered spectrum f according to the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1 to 6, the present invention provides a technical solution: a well constraint earthquake low-frequency recovery method based on compressed sensing comprises the following steps:
a2, extracting global static seismic wavelets and constructing a measurement matrix equation WF;
a3, initialization parameters: giving an initial reflection coefficient r0, the maximum iteration number Lin of an internal loop, the maximum iteration number Lout of an external loop and a residual error epsilon arbitrarily, wherein the initial iteration number k is 1;
a4、r0r0, outer loop: judge G-WFrk||2>εand k≤Lout,r0=rk-1And internal circulation: j is 0to Lin-1, and the reflection coefficient r is r at the time of terminationk,
a5, obtaining a spectrum R after frequency extension: r ═ Fr;
a6, obtaining a recovered spectrum G: g ═ L (r) + H (G), L, H are operators for extracting low-frequency and high-frequency spectral values, respectively;
a7, seismic data after low frequency recovery: g ═ F-1(G)。
TλAs a function of the threshold value, said (WF)TIs the conjugate transpose of WF, the WF is a measurement matrix composed of wavelet operator and Fourier transform operator, the measurement matrix is composed of wavelet operator and Fourier transform operatorIs RwThe conjugate transpose of RwExtracting position information of corresponding well for sampling operator, wherein alpha value is larger than (WF)TThe maximum characteristic value of WF, and the lambda is the sparsity of the adjusting reflection coefficient.
Inner loop and outer loop execution: k is equal to k +1 and k is equal to k +1,the reflection coefficient r ═ rkAnd G is Fourier transform of the noisy data.
The method comprises the following operation steps:
S2, extracting the global static seismic wavelet and constructing a measurement matrix equation WF;
s3, initialization parameters: giving an initial reflection coefficient r0, the maximum iteration number Lin of an internal loop, the maximum iteration number Lout of an external loop and a residual error epsilon arbitrarily, wherein the initial iteration number k is 1;
S4、r0r0, outer loop: judge G-WFrk||2>εand k≤Lout,r0=rk-1And internal circulation: j is 0to Lin-1, and the reflection coefficient r is r at the time of terminationk,
S5, obtaining a spectrum R after frequency broadening: r ═ Fr;
s6, obtaining a restored frequency spectrum G: g ═ L (r) + H (G), L, H are operators for extracting low-frequency and high-frequency spectral values, respectively;
s7, seismic data after low-frequency recovery: g ═ F-1(G)。
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are only illustrative of the present invention and are not intended to limit the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (3)
1. A well constraint earthquake low-frequency recovery method based on compressed sensing comprises the following steps:
a2, extracting global static seismic wavelets and constructing a measurement matrix equation WF;
a3, initialization parameters: giving an initial reflection coefficient r0, the maximum iteration number Lin of an internal loop, the maximum iteration number Lout of an external loop and a residual error epsilon arbitrarily, wherein the initial iteration number k is 1;
a4、r0r0, outer loop: judge G-WFrk||2>εand k≤Lout,r0=rk-1And internal circulation: j is 0to Lin-1, and the reflection coefficient r is r at the time of terminationk,
a5, obtaining a spectrum R after frequency extension: r ═ Fr;
a6, obtaining a recovered spectrum G: g ═ L (r) + H (G), L, H are operators for extracting low-frequency and high-frequency spectral values, respectively;
a7, seismic data after low frequency recovery: g ═ F-1(G)。
2. The method for seismic low frequency recovery based on compressive sensing of well constraints as defined in claim 1, wherein: the T isλAs a function of the threshold value, said (WF)TIs the conjugate transpose of WF, the WF is a measurement matrix composed of wavelet operator and Fourier transform operator, the measurement matrix is composed of wavelet operator and Fourier transform operatorIs RwThe conjugate transpose of RwExtracting position information of corresponding well for sampling operator, wherein alpha value is larger than (WF)TThe maximum characteristic value of WF, and the lambda is the sparsity of the adjusting reflection coefficient.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011546748.3A CN112558161B (en) | 2020-12-24 | 2020-12-24 | Well constraint earthquake low-frequency recovery method based on compressed sensing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011546748.3A CN112558161B (en) | 2020-12-24 | 2020-12-24 | Well constraint earthquake low-frequency recovery method based on compressed sensing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112558161A true CN112558161A (en) | 2021-03-26 |
CN112558161B CN112558161B (en) | 2022-07-26 |
Family
ID=75032362
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011546748.3A Active CN112558161B (en) | 2020-12-24 | 2020-12-24 | Well constraint earthquake low-frequency recovery method based on compressed sensing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112558161B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102169190A (en) * | 2011-01-06 | 2011-08-31 | 中国科学院地质与地球物理研究所 | Well-constrained pre-stack elastic parameter inversing method for modulating supplemented subspace |
US20130236115A1 (en) * | 2012-03-09 | 2013-09-12 | U.S. Government As Represented By The Secretary Of The Army | Method and System for Recovery of Missing Spectral Information in Wideband Signal |
CN107015274A (en) * | 2017-04-12 | 2017-08-04 | 中国石油大学(华东) | One kind missing seismic exploration data recovery and rebuilding method |
CN109946739A (en) * | 2019-03-15 | 2019-06-28 | 成都理工大学 | A kind of seismic profile Enhancement Method based on compressive sensing theory |
CN110174702A (en) * | 2018-09-30 | 2019-08-27 | 中海油海南能源有限公司 | A kind of method and system that marine seismic data low frequency weak signal is restored |
CN110308483A (en) * | 2019-05-23 | 2019-10-08 | 中国石油天然气股份有限公司 | Reflection coefficient acquiring method and device based on multitask Bayes's compressed sensing |
WO2020146322A1 (en) * | 2019-01-08 | 2020-07-16 | Westerngeco Llc | Combining noise attenuation and wavefield reconstruction in seismic processing |
-
2020
- 2020-12-24 CN CN202011546748.3A patent/CN112558161B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102169190A (en) * | 2011-01-06 | 2011-08-31 | 中国科学院地质与地球物理研究所 | Well-constrained pre-stack elastic parameter inversing method for modulating supplemented subspace |
US20130236115A1 (en) * | 2012-03-09 | 2013-09-12 | U.S. Government As Represented By The Secretary Of The Army | Method and System for Recovery of Missing Spectral Information in Wideband Signal |
CN107015274A (en) * | 2017-04-12 | 2017-08-04 | 中国石油大学(华东) | One kind missing seismic exploration data recovery and rebuilding method |
CN110174702A (en) * | 2018-09-30 | 2019-08-27 | 中海油海南能源有限公司 | A kind of method and system that marine seismic data low frequency weak signal is restored |
WO2020146322A1 (en) * | 2019-01-08 | 2020-07-16 | Westerngeco Llc | Combining noise attenuation and wavefield reconstruction in seismic processing |
CN109946739A (en) * | 2019-03-15 | 2019-06-28 | 成都理工大学 | A kind of seismic profile Enhancement Method based on compressive sensing theory |
CN110308483A (en) * | 2019-05-23 | 2019-10-08 | 中国石油天然气股份有限公司 | Reflection coefficient acquiring method and device based on multitask Bayes's compressed sensing |
Non-Patent Citations (2)
Title |
---|
ZHANG JUN-HUA ET AL.: "Low-frequency data analysis and expansion", 《APPLIED GEOPHYSICS》 * |
孙苗苗等: "基于压缩感知的加权MCA地震数据重构方法", 《地球物理学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112558161B (en) | 2022-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Hankel low-rank approximation for seismic noise attenuation | |
Qiu et al. | Deep learning prior model for unsupervised seismic data random noise attenuation | |
CA2947678C (en) | Compressive sensing | |
CN110174702B (en) | Method and system for recovering low-frequency weak signals of marine seismic data | |
CN110031899B (en) | Weak signal extraction algorithm based on compressed sensing | |
CN107179550B (en) | A kind of seismic signal zero phase deconvolution method of data-driven | |
CN109143331B (en) | Seismic wavelet extraction method | |
CN110646841B (en) | Time-varying sparse deconvolution method and system | |
CN111399057B (en) | Seismic data noise suppression method based on non-convex sparse constraint | |
Wu et al. | Preliminary study on Dreamlet based compressive sensing data recovery | |
CN110687597B (en) | Wave impedance inversion method based on joint dictionary | |
CN110967735A (en) | Self-adaptive ghost wave suppression method and system | |
CN112558161B (en) | Well constraint earthquake low-frequency recovery method based on compressed sensing | |
CN110749923A (en) | Deconvolution method for improving resolution based on norm equation | |
Liu et al. | An improved Gaussian frequency domain sparse inversion method based on compressed sensing | |
Cao et al. | REDUNDANCY, FAST CURVELET TRANSFORM | |
LIU et al. | Iterative seismic random noise suppression method based on compressive sensing | |
Sun et al. | Seismic data denoising with correlation feature optimization via S-mean | |
CN110568489A (en) | Broadband inversion method of bulk medium | |
Wang et al. | A Joint Framework for Seismic Signal Denoising Using Total Generalized Variation and Shearlet Transform | |
Li et al. | Enhancement of the seismic data resolution through Q-compensated denoising based on dictionary learning | |
DONG et al. | Fast implementation technique of adaptive Kalman filtering deconvolution via dyadic wavelet transform | |
CN109085649B (en) | Seismic data denoising method based on wavelet transformation optimization | |
Chen et al. | Adaptive distributed acoustic sensing coupling noise removal method | |
CN117055107B (en) | Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20220629 Address after: 266555 Qian Tang Jiang Road, Huangdao District, Qingdao, Shandong 369 Applicant after: QINGDAO TECHNICAL College Applicant after: China University of Petroleum (East China) Address before: 266580 No. 66 Changjiang West Road, Huangdao District, Qingdao, Shandong. Applicant before: CHINA University OF PETROLEUM (EAST CHINA) |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |