CN117055107B - Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms - Google Patents
Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms Download PDFInfo
- Publication number
- CN117055107B CN117055107B CN202310786024.3A CN202310786024A CN117055107B CN 117055107 B CN117055107 B CN 117055107B CN 202310786024 A CN202310786024 A CN 202310786024A CN 117055107 B CN117055107 B CN 117055107B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- framelet
- seismic
- interpolation
- pseudo
- 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
- 238000000034 method Methods 0.000 title claims abstract description 99
- 230000009466 transformation Effects 0.000 title claims abstract description 37
- 230000003993 interaction Effects 0.000 title claims abstract description 16
- 230000009977 dual effect Effects 0.000 claims abstract description 25
- 238000012545 processing Methods 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000008602 contraction Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 238000002939 conjugate gradient method Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims 2
- 230000000694 effects Effects 0.000 abstract description 13
- 230000002411 adverse Effects 0.000 abstract description 7
- 238000013459 approach Methods 0.000 description 4
- 238000000844 transformation Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method based on the Framelet transformation and l p A pseudo-norm interaction seismic interpolation method relates to the field of oil and gas exploration and seismic data processing, and solves the problem that details of seismic data are ignored in the existing sparse transform-based method, and further adverse effects are caused on processing and interpretation of high-resolution data. The invention includes constructing a Framelet transform function F; performing the Framelet transformation on the seismic data, and performing the l on the seismic data after the Framelet transformation p Pseudo-norm constraint, constructing an objective function of seismic interpolation; determining and updating update formulas of the seismic data, the Lagrangian multiplier and the dual variables; and determining a final seismic data interpolation result according to the relation among the seismic data, the Lagrangian multiplier and the dual variables. The method is used for seismic data interpolation.
Description
Technical Field
The invention relates to the field of oil and gas exploration and seismic data processing, in particular to a method based on the Framelet transformation and the seismic data processingA pseudo-norm interaction seismic interpolation method.
Background
In seismic exploration, data acquisition is the most basic step. The quality of the acquired data affects the understanding of the subsurface structure. However, due to the limitation of the acquisition conditions, the problem of missing traces cannot be ignored, especially in areas with complex geological conditions, such as the Qiang pond basin. Incomplete seismic data will impair subsequent steps such as amplitude and offset analysis and interpretation of results. Therefore, seismic data reconstruction is necessary.
The seismic data interpolation method based on the sparse transformation is rapid in development. This approach does not require a priori information required for wave equation based approaches, and is relatively simple to calculate compared to predictive filtering based approaches. This approach assumes that the seismic data is sparse in a certain region, and missing traces break the sparsity. They reconstruct the dataThe problem is considered an optimization problem with sparse constraints. Constraint is generally associated withThe norms are related. Common transformations include Fourier transforms, radon transforms, curvelet transforms, wavelet transforms, seilet transforms, streamlet transforms, and shearlet transforms. These transformations provide a good tool for analyzing sparsity of seismic data. While these transformations have been widely used in practical processing, they tend to ignore details of the seismic data. Such loss of detail may adversely affect the processing and interpretation of the high resolution data.
Disclosure of Invention
The invention aims at: the invention provides a method based on the Framelet transformation sumThe seismic interpolation method based on pseudo-norm interaction solves the problem that details of seismic data are ignored in the existing sparse transform-based method, and further adverse effects are generated on processing and interpretation of high-resolution data, and the interpolation effect is improved.
The technical scheme adopted by the invention for achieving the purpose is as follows:
step 1: building a Framelet transform function;
Step 2: performing the Framelet transformation on the seismic data, and performing the Framelet transformation on the seismic dataPseudo-norm constraint, constructing an objective function of seismic interpolation;
step 3: determining seismic dataLagrangian multiplier term->And dual variables>Updating formulas of the three and +.>Lagrangian multiplier term->And dual variables>Updating;
processing the seismic data update sub-problem by using the cost Ma Yinli to finish the process ofIs updated according to the update of (a);
updating sub-problem of Lagrangian multiplier sub-item to complete pairing according to generalized contraction methodIs updated according to the update of (a);
the update sub-problem of the pair-pair variable is paired according to the alternate direction multiplier methodIs updated according to the update of (a);
step 4: and determining a final seismic data interpolation result according to the relation among the seismic data, the Lagrangian multiplier and the dual variables.
Step 1 comprises the following steps:
step 1.1: order theRepresenting a scale function->And->Representing two wavelet functions, < >>Representing scale factors, scale functionsThe specific definitions of the numbers and wavelet functions are as follows equation (1) and equation (2):
(1)
(2)
wherein the method comprises the steps of,/>And->Low-pass and high-pass filters for scale functions and wavelet functions, < >>Meaning 0,1,2;
step 1.2: from equation (1) and equation (2), the equation for the Framelet transform can be expressed as:
is a signal>And->A low frequency part and a high frequency part of the scale function and the wavelet function, respectively, < >>Representation->sub-Framelet transform +_>And->The specific expression of (2) is as follows:
(4)
(5)
step 2 comprises the following steps:
step 2.1: the relationship between the complete seismic data and the observations is represented by the following equation (6):
(6)
in the method, in the process of the invention,representing the complete seismic data, which is also the result expected by interpolation,/for example>Representing the observation result->Representing a sampling matrix, which may be a diagonal matrix;
step 2.2: by usingThe pseudo-norm is used as a main body of a reconstruction method, sparse constraint is carried out on the seismic data after the Framelet transformation, and an objective function of seismic interpolation is constructed, wherein the objective function is shown in the following formula (7):
(7)
wherein the method comprises the steps ofIs Lagrangian factor, +.>For the Framelet transform +_>,/>Representation->Norms (F/F)>Representation ofPseudo-norms;
step 2.3: solving equation (7) using the alternate direction multiplier method, equation (7) becomes whereinRepresentation ofThe Lagrangian multiplier term:
(8)
equation (8) is further expressed as:
(9)
is a dual variable, ++>Is->After that, the formula (9) is decomposed into three sub-problems as the following formulas (10) - (12), -, respectively>Is a dual parameter, its value is defined by itself;
(10)
(11)
(12)
step 3 comprises the following steps:
step 3.1: equation (10) is typicalThe norm problem, let equation (13) equal 0, can be applied to +.>Updating is performed as shown in the following formula (14):
(13)
(14)
wherein the method comprises the steps ofIs unit momentArray (S)>The method is a transpose of the Framelet transformation, so that the method is simplified, and is convenient to solve by using a conjugate gradient method;
step 3.2: solving the formula (11) by adopting a generalized contraction methodIs updated by the following formula:
(15)
step 3.3: according to the alternate direction multiplier method, can obtainThe updated formula of (c) should be:
(16)
step 4 comprises the following steps:
step 4.1: judgingIf true, let seismic data +.>Lagrangian multiplier term->Dual variable->The method comprises the steps of carrying out a first treatment on the surface of the If not, output ++>As a final interpolation result.
The invention has the following beneficial effects because the technical scheme is adopted:
the invention adopts the Framelet transformation to groundTransforming the seismic data by adoptingThe pseudo-norm is used as a main body of a reconstruction method, sparse constraint is carried out on the seismic data after the Framelet transformation, and a method based on the Framelet transformation and ++is provided by combining an alternate direction multiplier method>A pseudo-norm interaction seismic interpolation method. The method solves the problem that details of the seismic data are ignored in the existing sparse transform-based method, and further adverse effects are generated on processing and interpretation of high-resolution data, and achieves the effect of improving the signal-to-noise ratio of seismic data reconstruction.
Drawings
FIG. 1 shows a frame-based transform sum according to the present inventionA flow chart of a pseudo-norm interaction seismic interpolation method;
FIG. 2 is a seismic section of raw seismic data of the invention
FIG. 3 is a seismic section of the invention with a random loss of 40% of the seismic data;
FIG. 4 is an interpolation of a seismic section of the invention with a random loss of 40% of the seismic data;
fig. 5 is a residual of the interpolation result of the present invention.
Detailed Description
Hereinafter, embodiments of the present invention will be described in detail. While the invention will be described and illustrated in conjunction with certain specific embodiments, it will be understood that it is not intended to limit the invention to these embodiments alone. On the contrary, the invention is intended to cover modifications and equivalent arrangements included within the scope of the appended claims.
In addition, numerous specific details are set forth in the following description in order to provide a better illustration of the invention. It will be understood by those skilled in the art that the present invention may be practiced without these specific details.
The invention is described in detail below in connection with fig. 1-5.
The invention solves the technical problems that: the method solves the problem that details of the seismic data are ignored in the existing sparse transform-based method, and further adverse effects are generated on processing and interpretation of high-resolution data, and achieves the effect of improving the signal-to-noise ratio of seismic data reconstruction.
Technical means: based on Framelet transformation sumA pseudo-norm interaction seismic interpolation method comprises the following steps of
Step 1: building a Framelet transform function;
Step 2: performing the Framelet transformation on the seismic data, and performing the Framelet transformation on the seismic dataPseudo-norm constraint, constructing an objective function of seismic interpolation;
step 3: determining seismic dataLagrangian multiplier term->And dual variables>Updating formulas of the three and +.>Lagrangian multiplier term->And dual variables>Updating;
step 4: and determining a final seismic data interpolation result according to the relation among the seismic data, the Lagrangian multiplier and the dual variables.
Step 1 comprises the following steps:
step 1.1: order theRepresenting a scale function->And->Representing two wavelet functions, < >>The specific definitions of scale factors, scale functions and wavelet functions are represented by the following formulas (1) and (2):
(1)
(2)
wherein the method comprises the steps of,/>And->Low-pass and high-pass filters for scale functions and wavelet functions, < >>Meaning 0,1,2;
step 1.2: from equation (1) and equation (2), the equation for the Framelet transform can be expressed as:
is a signal>And->A low frequency part and a high frequency part of the scale function and the wavelet function, respectively, < >>Representation->sub-Framelet transform +_>And->The specific expression of (2) is as follows:
(4)
(5)
step 2 comprises the following steps:
step 2.1: the relationship between the complete seismic data and the observations is represented by the following equation (6):
(6)
in the method, in the process of the invention,representing the complete seismic data, which is also the result expected by interpolation,/for example>Representing the observation result->Representing a sampling matrix, which may be a diagonal matrix;
step 2.2: by usingThe pseudo-norm is used as a main body of a reconstruction method, sparse constraint is carried out on the seismic data after the Framelet transformation, and an objective function of seismic interpolation is constructed, wherein the objective function is shown in the following formula (7):
(7)
wherein the method comprises the steps ofIs Lagrangian factor, +.>For the Framelet transform +_>,/>Representation->Norms (F/F)>Representation ofPseudo-norms;
step 2.3: solving equation (7) using the alternate direction multiplier method, equation (7) becomes whereinRepresentation ofThe Lagrangian multiplier term:
(8)
equation (8) is further expressed as:
(9)
is a dual variable, ++>Is->After that, the formula (9) is decomposed into three sub-problems as the following formulas (10) - (12), -, respectively>Is a dual parameter, its value is defined by itself;
(10)
(11)
(12)
step 3 comprises the following steps:
step 3.1: equation (10) is typicalThe norm problem, let equation (13) equal 0, can be applied to +.>Updating is performed as shown in the following formula (14):
(13)
(14)
wherein the method comprises the steps ofIs a unitary matrix->The method is a transpose of the Framelet transformation, so that the method is simplified, and is convenient to solve by using a conjugate gradient method;
step 3.2: solving the formula (11) by adopting a generalized contraction methodIs updated by the following formula:
(15)
step 3.3: according to the alternate direction multiplier method, can obtainThe updated formula of (c) should be:
(16)
step 4 comprises the following steps:
step 4.1: judgingIf true, let seismic data +.>Lagrangian multiplier term->Dual variable->The method comprises the steps of carrying out a first treatment on the surface of the If not, output ++>As a final interpolation result.
The technical effects are as follows: the invention adopts the Framelet transformation to transform the seismic data, adopts the following steps ofThe pseudo-norm is used as a main body of a reconstruction method, sparse constraint is carried out on the seismic data after the Framelet transformation, and a method based on the Framelet transformation and ++is provided by combining an alternate direction multiplier method>A pseudo-norm interaction seismic interpolation method. The method solves the problem that details of the seismic data are ignored in the existing sparse transform-based method, and further adverse effects are generated on processing and interpretation of high-resolution data, and achieves the effect of improving the signal-to-noise ratio of seismic data reconstruction.
Example 1
As shown in fig. 1-5, the seismic data interpolation is as follows:
step 1.1: order theRepresenting a scale function->And->Representing two wavelet functions, < >>The specific definitions of scale factors, scale functions and wavelet functions are represented by the following formulas (1) and (2):
(1)
(2)
wherein the method comprises the steps of,/>And->Low-pass and high-pass filters for scale functions and wavelet functions, < >>Meaning 0,1,2;
step 1.2: from equation (1) and equation (2), the equation for the Framelet transform can be expressed as:
is a signal>And->A low frequency part and a high frequency part of the scale function and the wavelet function, respectively, < >>Representation->sub-Framelet transform +_>And->The specific expression of (2) is as follows:
(4)
(5)
step 2.1: the relationship between the complete seismic data and the observations is represented by the following equation (6):
(6)
in the method, in the process of the invention,representing the complete seismic data, which is also the result expected by interpolation,/for example>Representing the observation result->Representing a sampling matrix, which may be a diagonal matrix;
step 2.2: by usingThe pseudo-norm is used as a main body of a reconstruction method, sparse constraint is carried out on the seismic data after the Framelet transformation, and an objective function of seismic interpolation is constructed, wherein the objective function is shown in the following formula (7):
(7)
wherein the method comprises the steps ofIs Lagrangian factor, +.>For the Framelet transform +_>,/>Representation->Norms (F/F)>Representation ofPseudo-norms;
step 2.3: solving equation (7) using the alternate direction multiplier method, equation (7) becomes whereinRepresentation ofThe Lagrangian multiplier term:
(8)
equation (8) is further expressed as:
(9)
is a dual variable, ++>Is->After that, the formula (9) is decomposed into three sub-problems as the following formulas (10) - (12), -, respectively>Is a dual parameter, its value is defined by itself;
(10)
(11)
(12)
step 3.1: equation (10) is typicalThe norm problem, let equation (13) equal 0, can be applied to +.>Updating is performed as shown in the following formula (14):
(13)
(14)
wherein the method comprises the steps ofIs a single sheetBit matrix,/->The method is a transpose of the Framelet transformation, so that the method is simplified, and is convenient to solve by using a conjugate gradient method;
step 3.2: solving the formula (11) by adopting a generalized contraction methodIs updated by the following formula:
(15)
step 3.3: according to the alternate direction multiplier method, can obtainThe updated formula of (c) should be:
(16)
step 4.1: judgingIf true, let seismic data +.>Lagrangian multiplier term->Dual variable->The method comprises the steps of carrying out a first treatment on the surface of the If not, output ++>As a final interpolation result.
And (3) effect analysis: as shown in fig. 2-5, the interpolated result preferably reconstructs missing seismic data, which is closer to the original seismic data, thus proving the correctness of the method. The invention utilizes the Framelet transformation andthe pseudo-norm interaction method performs interpolation processing on the seismic data. The invention adopts the Framlet transformation to transform the seismic data, adopts +.>The pseudo-norm is used as a main body of a reconstruction method, sparse constraint is carried out on the seismic data after the Framelet transformation, and a method based on the Framelet transformation and ++is provided by combining an alternate direction multiplier method>A pseudo-norm interaction seismic interpolation method. The method solves the problem that details of the seismic data are ignored in the existing sparse transform-based method, and further adverse effects are generated on processing and interpretation of high-resolution data, and achieves the effect of improving the signal-to-noise ratio of seismic data reconstruction.
Claims (5)
1. Based on Framelet transformation sumA method of pseudo-norm interaction seismic interpolation comprising the steps of:
step 1: building a Framelet transform function;
Step 2: performing the Framelet transformation on the seismic data, and performing the Framelet transformation on the seismic dataThe pseudo-norm constraint is adopted, an objective function of seismic interpolation is constructed, and the objective function is decomposed into three sub-problems, namely a seismic data updating sub-problem, a Lagrangian multiplier updating sub-problem and a dual variable updating sub-problem;
step 3: determining seismic dataLagrangian multiplier term->And dual variables>Updating formulas of the three and +.>Lagrangian multiplier term->And dual variables>Updating;
processing the seismic data update sub-problem by using the cost Ma Yinli to finish the process ofIs updated according to the update of (a);
updating sub-problem of Lagrangian multiplier sub-item to complete pairing according to generalized contraction methodIs updated according to the update of (a);
the update sub-problem of the dual variable is completed according to the alternate direction multiplication methodIs updated according to the update of (a);
step 4: and determining a final seismic data interpolation result according to the relation among the updated seismic data, the Lagrangian multiplier and the dual variables.
2. A Framelet transform sum based method according to claim 1Pseudo-norm interaction groundThe vibration interpolation method is characterized in that: the step 1 comprises the following steps:
step 1.1: order theRepresenting a scale function->And->Representing two wavelet functions, < >>The specific definitions of scale factors, scale functions and wavelet functions are represented by the following formulas (1) and (2):
(1)
(2)
wherein the method comprises the steps of,/>And->Low-pass and high-pass filters for scale functions and wavelet functions, < >>Meaning 0,1,2;
step 1.2: from equation (1) and equation (2), the equation for the Framelet transform can be expressed as:
(3)
is a signal>And->A low frequency part and a high frequency part of the scale function and the wavelet function, respectively, < >>Representation->sub-Framelet transform +_>And->The specific expression of (2) is as follows:
(4)
(5)。
3. a Framelet transform sum based method according to claim 1The seismic interpolation method of pseudo-norm interaction is characterized in that: the step 2 comprises the following steps:
step 2.1: the relationship between the complete seismic data and the observations is represented by the following equation (6):
(6)
in the method, in the process of the invention,representing the complete seismic data, which is also the result expected by interpolation,/for example>Representing the observation result->Representing a sampling matrix, the sampling matrix being a diagonal matrix;
step 2.2: by usingThe pseudo-norm is used as a main body of a reconstruction method, sparse constraint is carried out on the seismic data after the Framelet transformation, and an objective function of seismic interpolation is constructed, wherein the objective function is shown in the following formula (7):
(7)
wherein the method comprises the steps ofIs Lagrangian factor, +.>For the Framelet transform +_>,/>Representation->Norms (F/F)>Representation->Pseudo-norms;
step 2.3: solving equation (7) using the alternate direction multiplier method, equation (7) becomes whereinRepresentation->The Lagrangian multiplier term:
(8)
equation (8) is further expressed as:
(9)
is a dual variable, ++>Is->After that, the formula (9) is decomposed into three sub-problems as the following formulas (10) - (12), -, respectively>Is a dual parameter, its value is defined by itself;
(10)
(11)
(12)。
4. a Framelet transform sum based method according to claim 3The seismic interpolation method of pseudo-norm interaction is characterized in that: the step 3 comprises the following steps:
step 3.1: equation (10) is typicalThe norm problem, let equation (13) equal 0, can be applied to +.>Updating is performed as shown in the following formula (14):
(13)
(14)
the formula (13) expands the formula (10), and derives m from the expanded formula, and specifically explains the updating step of m;
wherein the method comprises the steps ofIs a unitary matrix->The method is a transpose of the Framelet transformation, so that the formula (10) is simplified, and the method is convenient to solve by using a conjugate gradient method;
step 3.2: solving the formula (11) by adopting a generalized contraction methodIs updated by the following formula:
(15)
step 3.3: according to the alternate direction multiplier method, can obtainThe updated formula of (c) should be:
(16)。
5. a Framelet transform sum based method according to claim 1The seismic interpolation method of pseudo-norm interaction is characterized in that: the step 4 comprises the following steps:
step 4.1: judgingIf true, let seismic data +.>Pulling upGelang Ri multiplier sub-term->Dual variable->The method comprises the steps of carrying out a first treatment on the surface of the If not, output ++>As a final interpolation result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310786024.3A CN117055107B (en) | 2023-06-29 | 2023-06-29 | Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310786024.3A CN117055107B (en) | 2023-06-29 | 2023-06-29 | Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117055107A CN117055107A (en) | 2023-11-14 |
CN117055107B true CN117055107B (en) | 2024-04-09 |
Family
ID=88661494
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310786024.3A Active CN117055107B (en) | 2023-06-29 | 2023-06-29 | Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117055107B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4594693A (en) * | 1983-11-04 | 1986-06-10 | Mobil Oil Corporation | Seismic trace interpolation using f-k filtering |
CN110261912A (en) * | 2019-07-23 | 2019-09-20 | 河北地质大学 | The interpolation and denoising method and system of a kind of seismic data |
CN112526599A (en) * | 2019-09-17 | 2021-03-19 | 中国石油化工股份有限公司 | Wavelet phase estimation method and system based on weighted L1 norm sparsity criterion |
CN115469364A (en) * | 2022-10-21 | 2022-12-13 | 成都理工大学 | Earthquake denoising method and system based on stable Framelet transformation |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2784551A3 (en) * | 2013-03-26 | 2015-10-28 | CGG Services SA | System and method for interpolating seismic data by matching pursuit in fourier transform |
-
2023
- 2023-06-29 CN CN202310786024.3A patent/CN117055107B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4594693A (en) * | 1983-11-04 | 1986-06-10 | Mobil Oil Corporation | Seismic trace interpolation using f-k filtering |
CN110261912A (en) * | 2019-07-23 | 2019-09-20 | 河北地质大学 | The interpolation and denoising method and system of a kind of seismic data |
CN112526599A (en) * | 2019-09-17 | 2021-03-19 | 中国石油化工股份有限公司 | Wavelet phase estimation method and system based on weighted L1 norm sparsity criterion |
CN115469364A (en) * | 2022-10-21 | 2022-12-13 | 成都理工大学 | Earthquake denoising method and system based on stable Framelet transformation |
Non-Patent Citations (1)
Title |
---|
常凯 ; 张海江 ; 林叶 ; .基于样条插值与曲波变换压缩感知的井下微地震监测数据重建.物探化探计算技术.2016,(第06期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN117055107A (en) | 2023-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gao et al. | Irregular seismic data reconstruction based on exponential threshold model of POCS method | |
AU2007336115B2 (en) | Regularisation of irregularly sampled seismic data | |
CN110174702B (en) | Method and system for recovering low-frequency weak signals of marine seismic data | |
CN105549076B (en) | A kind of seismic data processing technique based on alternating direction method and full Theory of Variational Principles | |
Huang et al. | Robust seismic image interpolation with mathematical morphological constraint | |
CN112904426A (en) | Decoupling elastic wave reverse time migration method, system and application | |
Gan* et al. | Seismic data reconstruction via fast projection onto convex sets in the seislet transform domain | |
CN113608259B (en) | Seismic thin layer detection method based on ICEEMDAN constraint generalized S transformation | |
Yang et al. | Seislet-based morphological component analysis using scale-dependent exponential shrinkage | |
CN110646841B (en) | Time-varying sparse deconvolution method and system | |
CN111399057B (en) | Seismic data noise suppression method based on non-convex sparse constraint | |
Song et al. | An efficient wavefield inversion for transversely isotropic media with a vertical axis of symmetry | |
CN110244353B (en) | Seismic data regularization method based on sparse norm optimization algorithm | |
CN117055107B (en) | Seismic interpolation method based on interaction of Framelet transformation and Lp pseudo-norms | |
Min et al. | Fast convex set projection with deep prior for seismic interpolation | |
Wang et al. | Robust nonstationary local slope estimation | |
US9091789B2 (en) | Seismic signal processing method with Gaussian slowness-period packets | |
CN115469364B (en) | Seismic denoising method and system based on stable Framelet transformation | |
Bai et al. | Iterative sparse deconvolution using seislet-domain constraint | |
AU2011248987A1 (en) | Seismic signal processing method with Gaussian slowness-period packets | |
CN112149614B (en) | Pre-stack well-seismic combined intelligent denoising method | |
Pan et al. | Automatic Seismic Lithology Interpretation via Multi-attribute Integrated Deep Learning | |
Hao et al. | Nonstationary seismic inversion: joint estimation for acoustic impedance, attenuation factor and source wavelet | |
CN108345034B (en) | Seismic data regularization method | |
Pan et al. | The Interplay of Framelet Transform and l p Quasi-Norm to Interpolate 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 |