Empirical mode decomposition method of seismic section

Abstract

The invention relates to the field of oil gas seismic exploration and particularly relates to an empirical mode decomposition method of a seismic section. The method relates to two-dimensional instantaneous attribute analysis of the seismic section and is a computing method for decomposing a single-component two-dimensional intrinsic mode function from the seismic section. In a two-dimensional empirical mode decomposition process, an extreme value detection method of two point adjacent regions of a time direction and an equality constrained tight framework region rarefaction curved surface interpolation method are adopted to solve an extreme value detection problem of the seismic section and an envelope surface interpolation problem of large-scale irregular data. According to the empirical mode decomposition method disclosed by the invention, a computing result is obviously better than that of an existing two-dimensional empirical mode decomposition rapid method, and the computing time is roughly equal.

Description

technical field

[0001] The invention relates to the field of oil and gas seismic exploration, and more specifically, relates to an empirical mode decomposition method of seismic sections in the analysis of seismic instantaneous attributes. Background technique

[0002] The traditional two-dimensional seismic section complex analysis technique (A.E.Barnes, Theory of 2D complex seismic trace analysis, Geophysics, Vol.61, No.1, 1996, 264–272) is to directly perform two-dimensional Hilbert transformation on the seismic section (Bidimensional Hilbert Transform, BHT), to form an analytical signal, and then calculate the instantaneous seismic characteristic attributes such as instantaneous amplitude, instantaneous phase, and instantaneous frequency, which reflect the micro-local changes in the process of seismic wave propagation, and can help determine the lithology, lithofacies changes and The properties of fluids contained in rock pores. However, this complex analysis technique ...

Examples

Embodiment Construction

[0029] The present invention will be further described below in conjunction with the accompanying drawings, but the embodiments of the present invention are not limited thereto.

[0030] A clear definition has been given for one-dimensional IMF, but there is no definition for BIMF. For a one-dimensional IMF, four conditions are met, the number of local maximum points and minimum points differs by at most one; there is a vibration mode, that is, there is only one maximum or pole between two zero-crossing points Small value; at any time, the mean line of the upper and lower envelopes is a zero function or tends to a zero function; the IMFs obtained by decomposing the same signal should be orthogonal to each other. Since the first two conditions are difficult to guarantee in the decomposition process, as long as the BIMF meets the last two conditions of the one-dimensional IMF in the calculation process, it can be used as an effective approximation of the two-dimensional single-c...

technical field

[0001] The invention relates to the field of oil and gas seismic exploration, and more specifically, relates to an empirical mode decomposition method of seismic sections in the analysis of seismic instantaneous attributes. Background technique

[0002] The traditional two-dimensional seismic section complex analysis technique (A.E.Barnes, Theory of 2D complex seismic trace analysis, Geophysics, Vol.61, No.1, 1996, 264–272) is to directly perform two-dimensional Hilbert transformation on the seismic section (Bidimensional Hilbert Transform, BHT), to form an analytical signal, and then calculate the instantaneous seismic characteristic attributes such as instantaneous amplitude, instantaneous phase, and instantaneous frequency, which reflect the micro-local changes in the process of seismic wave propagation, and can help determine the lithology, lithofacies changes and The properties of fluids contained in rock pores. However, this complex analysis technique ...

Examples

Embodiment Construction

[0029] The present invention will be further described below in conjunction with the accompanying drawings, but the embodiments of the present invention are not limited thereto.

[0030] A clear definition has been given for one-dimensional IMF, but there is no definition for BIMF. For a one-dimensional IMF, four conditions are met, the number of local maximum points and minimum points differs by at most one; there is a vibration mode, that is, there is only one maximum or pole between two zero-crossing points Small value; at any time, the mean line of the upper and lower envelopes is a zero function or tends to a zero function; the IMFs obtained by decomposing the same signal should be orthogonal to each other. Since the first two conditions are difficult to guarantee in the decomposition process, as long as the BIMF meets the last two conditions of the one-dimensional IMF in the calculation process, it can be used as an effective approximation of the two-dimensional single-c...

Images