阅读说明：本技术 频谱分解方法及装置 (Spectrum decomposition method and device ) 是由 林腾飞 王雪玲 田中元 罗贝维 刘杏芳 于 2021-01-21 设计创作，主要内容包括：本发明公开了一种频谱分解方法及装置,该频谱分解方法包括：根据预设工区的地震数据获取地震倾角数据。对地震数据进行频谱分解计算。根据频谱分解计算的结果,结合地震倾角数据获取频谱分解结果。本发明解决了频谱分解结果出现横向不稳定的问题,进而保证了后续储层厚度预测、储层识别等的作业精度。(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.)

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.

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:

step 101, acquiring earthquake inclination angle data according to earthquake data of a preset work area.

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＝(p²+q²)^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.