CN112835100A - Spectrum decomposition method and device - Google Patents
Spectrum decomposition method and device Download PDFInfo
- Publication number
- CN112835100A CN112835100A CN202110081366.6A CN202110081366A CN112835100A CN 112835100 A CN112835100 A CN 112835100A CN 202110081366 A CN202110081366 A CN 202110081366A CN 112835100 A CN112835100 A CN 112835100A
- Authority
- CN
- China
- Prior art keywords
- seismic
- data
- decomposition
- spectrum decomposition
- result
- 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.)
- Pending
Links
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 92
- 238000001228 spectrum Methods 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 34
- 230000003595 spectral effect Effects 0.000 claims abstract description 18
- 238000004590 computer program Methods 0.000 claims description 16
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 3
- 230000008021 deposition Effects 0.000 description 3
- 238000005056 compaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/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
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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 and a device for frequency spectrum decomposition, wherein the method for frequency spectrum decomposition comprises the following steps: and acquiring earthquake dip angle data according to the earthquake data of the preset work area. And performing spectral decomposition calculation on the seismic data. And according to the result of the spectrum decomposition calculation, combining the seismic dip angle data to obtain a spectrum decomposition result. The method solves the problem that the frequency spectrum decomposition result is transversely unstable, and further ensures the operation precision of subsequent reservoir thickness prediction, reservoir identification and the like.
Description
Technical Field
The invention relates to the technical field of oil and gas exploration, in particular to a frequency spectrum decomposition method and a frequency spectrum decomposition device.
Background
The frequency spectrum decomposition technology is an interpretation method of a frequency domain, is an important component in seismic attribute analysis, and is used for researching a thin layer structure and judging a deposition environment by using the change relation of frequency along with time.
At present, the related algorithms of the frequency spectrum decomposition technology are analyzed channel by channel according to seismic channels, but the frequency spectrum decomposition result is influenced by noise in seismic data in the process of analyzing the seismic channels channel by channel, so that the problem of transverse instability of the frequency spectrum decomposition result is caused, and the operation precision of subsequent reservoir thickness prediction, reservoir identification and the like is influenced.
Disclosure of Invention
The embodiment of the invention provides a frequency spectrum decomposition method, which is used for solving the problem that a frequency spectrum decomposition result is transversely unstable so as to ensure the operation precision of subsequent reservoir thickness prediction, reservoir identification and the like, and comprises the following steps:
acquiring earthquake dip angle data according to earthquake data of a preset work area;
performing spectral decomposition calculation on the seismic data;
and obtaining a spectrum decomposition result by combining the seismic dip angle data according to the result of spectrum decomposition calculation.
Optionally, obtaining a spectrum decomposition result by combining the seismic dip data according to a result of spectrum decomposition calculation, including:
configuring a filter according to the result of the frequency spectrum decomposition calculation and the seismic dip angle data;
and acquiring a spectrum decomposition result according to the configured filter.
Optionally, the seismic dip data includes: seismic dip information and seismic dip information.
Optionally, the method for performing spectrum decomposition calculation on the seismic data is a matching pursuit algorithm.
The embodiment of the invention also provides a frequency spectrum decomposition device, which is used for solving the problem that the frequency spectrum decomposition result is transversely unstable so as to ensure the operation precision of subsequent reservoir thickness prediction, reservoir identification and the like, and comprises the following components:
the data acquisition module is used for acquiring earthquake dip angle data according to earthquake data of a preset work area;
the calculation module is used for performing spectrum decomposition calculation on the seismic data;
and the frequency spectrum decomposition module is used for acquiring a frequency spectrum decomposition result by combining the seismic dip angle data according to a frequency spectrum decomposition calculation result.
Optionally, the spectrum decomposition module is further configured to:
configuring a filter according to the result of the frequency spectrum decomposition calculation and the seismic dip angle data;
and acquiring a spectrum decomposition result according to the configured filter.
Optionally, the seismic dip data includes: seismic dip information and seismic dip information.
Optionally, the method for performing spectrum decomposition calculation on the seismic data is a matching pursuit algorithm.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the method when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium, in which a computer program for executing the above method is stored.
According to the frequency spectrum decomposition method provided by the embodiment of the invention, the seismic dip angle data are obtained according to the seismic data of the preset work area, the frequency spectrum decomposition calculation is carried out on the seismic data, and then the frequency spectrum decomposition result is obtained by combining the seismic dip angle data according to the frequency spectrum decomposition calculation result, so that the frequency spectrum decomposition is carried out under the constraint of the seismic dip angle data, the problem that the frequency spectrum decomposition result is transversely unstable is solved, and the operation precision of subsequent reservoir thickness prediction, reservoir identification and the like is further ensured.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
FIG. 1 is a flow chart of a spectral decomposition method according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a spectrum decomposition apparatus according to an embodiment of the present invention;
FIG. 3 is a schematic seismic section of a post-stack seismic volume of a pre-engineered work area in an embodiment of the invention;
FIG. 4 is a schematic time slice of a pre-determined work area according to an embodiment of the present disclosure;
FIG. 5 is a schematic cross-sectional view of a preset work area data obtained by a prior art spectral decomposition method;
FIG. 6 is a schematic diagram of a preset work area data plane obtained by a prior art spectral decomposition method;
FIG. 7 is a schematic cross-sectional view of data of a predetermined work area according to an embodiment of the present invention;
FIG. 8 is a schematic plan view of a data plane of a predetermined work area according to an embodiment of the present invention;
FIG. 9 is a diagram illustrating an example of a time window in an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
Fig. 1 is a flowchart of a spectrum decomposition method according to an embodiment of the present invention, as shown in fig. 1, the method includes:
In this embodiment, the seismic dip data includes: seismic dip information and seismic dip information.
Specifically, the calculation formula of the seismic dip and the seismic dip is as follows:
seismic dip s:
s＝(p2+q2)1/2
earthquake tendency φ:
wherein p is the inclination angle in the x direction in the rectangular coordinate system, and q is the inclination angle in the y direction in the rectangular coordinate system.
And 102, carrying out spectrum decomposition calculation on the seismic data.
In this embodiment, in order to ensure high accuracy of spectral decomposition calculation, a matching pursuit algorithm may be used to perform spectral decomposition calculation on the seismic data.
And 103, acquiring a frequency spectrum decomposition result by combining the seismic dip angle data according to the frequency spectrum decomposition calculation result.
In this embodiment, step 103 includes:
configuring a filter according to the result of the frequency spectrum decomposition calculation and the seismic dip angle data;
and acquiring a spectrum decomposition result according to the configured filter.
In particular, the filter may be defined by a window function, as shown in figure 9,
are the 5 analysis points located within the analysis window function time window,
is that
α is a predetermined filter parameter.
As can be seen from fig. 1, according to the spectrum decomposition method provided in the embodiment of the present invention, the seismic dip data is obtained according to the seismic data of the preset work area, the spectrum decomposition calculation is performed on the seismic data, and then the spectrum decomposition result is obtained according to the result of the spectrum decomposition calculation by combining the seismic dip data, so that the spectrum decomposition is performed under the constraint of the seismic dip data, the problem of lateral instability of the spectrum decomposition result is solved, and the operation accuracy of subsequent reservoir thickness prediction, reservoir identification, and the like is ensured.
The invention is described below with specific examples of pre-set work zones:
FIGS. 3 and 4 show seismic sections and time slices, respectively, of a post-stack seismic volume. The arrows in fig. 3 indicate strong reflection points, which are caused by differential compaction of alluvial deposition systems, and the arrows in fig. 4 indicate the main river channels that are characteristic.
After spectral decomposition based on the match-track method, we have obtained the seismic profile and time slice of fig. 5 and 6, respectively, peak frequency. Since spectral analysis is sensitive to deposition anomalies, we find that the river has strong reflections, which can be demarcated in fig. 5 as river boundaries, characterized by a lower peak frequency response compared to the surrounding area. Meanwhile, due to the existence of spike anomaly, we can also barely track the lateral structure, and certainly, the spike anomaly can also distinguish the main river channel system in the time slice of fig. 6. Therefore, in seismic interpretation, it is necessary to implement horizontal constraints in order to impart immediate geological significance in seismic interpretation.
Fig. 7 and 8 show that the construction-constrained spectral decomposition is enhanced using different user-defined filtering parameters based on the match-and-trace approach. The strong reflections caused by differential compaction of the channel are retained in fig. 7, and the spike anomalies are greatly reduced (the user-defined filter parameter α is 0.25), which makes it possible to identify the channel in the time slice of fig. 8. If we set the user-defined filter parameter α to 1.0, we initialize much of the geological significance in FIG. 8 with the lateral variation constraint leads, which corresponds to the seismic section in FIG. 3. The main river in fig. 8 becomes the visible river, providing the interpreter with another reliable evidence for seismic interpretation.
Based on the same inventive concept, the embodiment of the present invention further provides a spectrum decomposition apparatus, as described in the following embodiments. Because the principle of the spectrum decomposition device for solving the problem is similar to the spectrum decomposition method, the implementation of the spectrum decomposition device can refer to the implementation of the spectrum decomposition method, and repeated details are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.
Fig. 2 is a schematic structural diagram of a spectrum decomposition apparatus according to an embodiment of the present invention, as shown in fig. 2, the apparatus includes:
the data acquisition module 201 is configured to acquire seismic dip angle data according to seismic data of a preset work area.
A calculating module 202, configured to perform a spectral decomposition calculation on the seismic data.
And the spectrum decomposition module 203 is used for acquiring a spectrum decomposition result by combining the seismic dip data according to the result of the spectrum decomposition calculation.
In an embodiment of the present invention, the spectrum decomposition module 203 is further configured to:
configuring a filter according to the result of the frequency spectrum decomposition calculation and the seismic dip angle data;
and acquiring a spectrum decomposition result according to the configured filter.
In an embodiment of the invention, the seismic dip data comprises: seismic dip information and seismic dip information.
In the embodiment of the invention, the method for performing spectrum decomposition calculation on the seismic data is a matching pursuit algorithm.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the above method when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium, in which a computer program for executing the above method is stored.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A method of spectral decomposition, comprising:
acquiring earthquake dip angle data according to earthquake data of a preset work area;
performing spectral decomposition calculation on the seismic data;
and obtaining a spectrum decomposition result by combining the seismic dip angle data according to the result of spectrum decomposition calculation.
2. The method of claim 1, wherein obtaining spectral decomposition results from the results of the spectral decomposition calculations in conjunction with the seismic dip data comprises:
configuring a filter according to the result of the frequency spectrum decomposition calculation and the seismic dip angle data;
and acquiring a spectrum decomposition result according to the configured filter.
3. The method of claim 1, wherein the seismic dip data comprises: seismic dip information and seismic dip information.
4. The method of claim 1, wherein the spectral decomposition of the seismic data is performed using a matching pursuit algorithm.
5. A spectral decomposition apparatus, comprising:
the data acquisition module is used for acquiring earthquake dip angle data according to earthquake data of a preset work area;
the calculation module is used for performing spectrum decomposition calculation on the seismic data;
and the frequency spectrum decomposition module is used for acquiring a frequency spectrum decomposition result by combining the seismic dip angle data according to a frequency spectrum decomposition calculation result.
6. The apparatus of claim 5, wherein the spectral decomposition module is further to:
configuring a filter according to the result of the frequency spectrum decomposition calculation and the seismic dip angle data;
and acquiring a spectrum decomposition result according to the configured filter.
7. The apparatus of claim 5, wherein the seismic dip data comprises: seismic dip information and seismic dip information.
8. The apparatus of claim 5, wherein the spectral decomposition of the seismic data is performed using a matching pursuit algorithm.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 4 when executing the computer program.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110081366.6A CN112835100A (en) | 2021-01-21 | 2021-01-21 | Spectrum decomposition method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110081366.6A CN112835100A (en) | 2021-01-21 | 2021-01-21 | Spectrum decomposition method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112835100A true CN112835100A (en) | 2021-05-25 |
Family
ID=75929287
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110081366.6A Pending CN112835100A (en) | 2021-01-21 | 2021-01-21 | Spectrum decomposition method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112835100A (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120143510A1 (en) * | 2007-05-25 | 2012-06-07 | Aftab Alam | High resolution attributes for seismic data processing and interpretation |
| CN103048690A (en) * | 2012-12-11 | 2013-04-17 | 成都理工大学 | Stratum reflection pickup technology of rapid matched projection decomposition extracted based on optimal seismic wavelet |
| CN104656132A (en) * | 2013-11-22 | 2015-05-27 | 中国石油天然气集团公司 | Method for improving storage efficiency of multi-azimuthal stratigraphic dip information |
| CN105353408A (en) * | 2015-11-20 | 2016-02-24 | 电子科技大学 | Wigner higher-order spectrum seismic signal spectral decomposition method based on matching pursuit |
| CN111487675A (en) * | 2020-03-25 | 2020-08-04 | 王仰华 | Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum |
-
2021
- 2021-01-21 CN CN202110081366.6A patent/CN112835100A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120143510A1 (en) * | 2007-05-25 | 2012-06-07 | Aftab Alam | High resolution attributes for seismic data processing and interpretation |
| CN103048690A (en) * | 2012-12-11 | 2013-04-17 | 成都理工大学 | Stratum reflection pickup technology of rapid matched projection decomposition extracted based on optimal seismic wavelet |
| CN104656132A (en) * | 2013-11-22 | 2015-05-27 | 中国石油天然气集团公司 | Method for improving storage efficiency of multi-azimuthal stratigraphic dip information |
| CN105353408A (en) * | 2015-11-20 | 2016-02-24 | 电子科技大学 | Wigner higher-order spectrum seismic signal spectral decomposition method based on matching pursuit |
| CN111487675A (en) * | 2020-03-25 | 2020-08-04 | 王仰华 | Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum |
Non-Patent Citations (1)
| Title |
|---|
| 杨海涛 等: "基于Wigner2Ville的谱分解效果分析", 《勘探地球物理进展》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Basir et al. | Multi-attribute ant-tracking and neural network for fault detection: A case study of an Iranian oilfield | |
| CN103454678B (en) | Method and system for determining seismic slice isochronism | |
| Ekinci et al. | On the effectiveness of directional derivative based filters on gravity anomalies for source edge approximation: synthetic simulations and a case study from the Aegean graben system (western Anatolia, Turkey) | |
| EA020070B1 (en) | Method for indexing a subsurface volume for the purpose of inferring geologic information | |
| CA2887175A1 (en) | System and method for processing 4d seismic data | |
| AU2014345427A1 (en) | Method and device for processing seismic signals | |
| CN110068860B (en) | Time-lapse seismic data processing method, device and system | |
| WO2015145195A1 (en) | Determining location and depth of subsurface magnetic sources | |
| AU2015290248A1 (en) | System and method for rock property estimation of subsurface geologic volumes | |
| CN108254785B (en) | Reservoir determination method and device | |
| CA2693247C (en) | Seismic attributes for structural analysis | |
| CN112835100A (en) | Spectrum decomposition method and device | |
| CN109581498B (en) | Reservoir thickness distribution determination method, system, device and readable medium | |
| CN112649887A (en) | Method and device for quantitatively dividing sequence stratum based on drilling data | |
| CN107656309B (en) | Method and device for determining fault distance | |
| Al-Dossary et al. | Estimating randomness using seismic disorder | |
| CN113900141B (en) | Oil-gas distribution prediction method and device | |
| CN114895372A (en) | Underground goaf detection comprehensive interpretation method and system | |
| CN112444871B (en) | Quantitative analysis method and equipment for crack spacing based on scattered wave seismic response characteristics | |
| CN112305597B (en) | Reservoir prediction method and device | |
| Alfonzo et al. | Analysis and calibration of 4D seismic data prior to 4D seismic inversion and history matching-Norne field case | |
| CN112363243A (en) | Method, device, computer equipment and medium for predicting high-quality reservoir of shale oil | |
| RU2789759C1 (en) | Method for determining the boundaries of the fractured zone | |
| CN113311485A (en) | Seismic sedimentary feature enhanced filtering method and device | |
| CN113325474B (en) | Method for discriminating biological reef |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210525 |