An oil and gas detection method and device based on instantaneous energy frequency value and electronic equipment

Abstract

The application provides an oil and gas detection method and device based on instantaneous energy frequency value, a computer readable storage medium and an electronic device. The method comprises the following steps: for each seismic signal in a three-dimensional seismic data volume of a target area, decomposing the seismic signal into a data set composed of a plurality of IMF components, extracting the instantaneous amplitude and instantaneous frequency of each IMF component, calculating the instantaneous energy frequency value of each IMF component based on the instantaneous amplitude and instantaneous frequency, and calculating the instantaneous energy frequency value of the seismic signal based on the instantaneous energy frequency value of each IMF component; analyzing the change characteristics of energy and frequency in the stratum of the target area based on the instantaneous energy frequency value of each seismic signal in the three-dimensional seismic data volume and the instantaneous energy frequency value of each IMF component of each seismic signal, and further analyzing the oil and gas bearing property of the stratum. The method provided by the application can more accurately identify the energy and frequency anomalies in the stratum and can also assist in oil and gas identification.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geophysical exploration, and in particular to an oil and gas detection method and device based on instantaneous energy frequency value, a computer readable storage medium and an electronic device. BACKGROUND

[0002] Seismic wave causes the attenuation of seismic wave energy in the process of propagating to the underground. The attenuation of seismic wave is classified into two categories. One is the attenuation related to the propagation characteristics of seismic wave, and the other is the intrinsic attenuation related to the propagation medium, such as the absorption attenuation caused after passing through the oil and gas bearing formation. The second kind of attenuation characteristics can be used for detecting the oil and gas bearing property. At present, the techniques for detecting the oil and gas bearing property in the field of geophysical exploration mainly include two categories of post-stack and pre-stack. The oil and gas bearing property detection techniques based on pre-stack seismic signals mainly include AVO technique and pre-stack seismic inversion technique. The reservoir or fluid can be identified by extracting various AVO attributes or inverting various fluid factors and making crossplot analysis. This kind of technique is greatly affected by pre-stack data and has a large amount of calculation, but the accuracy of oil and gas bearing property detection is relatively higher than that of post-stack. The oil and gas bearing property detection techniques based on post-stack seismic signals mainly use the seismic spectral decomposition technique to extract the absorption coefficient, attenuation gradient and various frequency attributes, or use the low frequency shadowing technique to detect oil and gas according to the characteristics of high frequency attenuation and low frequency increase. The complex seismic trace analysis technique based on post-stack seismic signals extracts the "sweet spot" attribute which is also widely used for identifying reservoir or fluid, and the "sweet spot" attribute is represented by the ratio of instantaneous amplitude to instantaneous frequency. However, the instantaneous frequency directly extracted from the original seismic signal by using Hilbert transform cannot accurately represent the frequency components contained in the signal, and E.N.Huang and others have fully analyzed the instantaneous frequency and considered that only the instantaneous frequency extracted from the intrinsic mode function is meaningful. In view of the problems existing in the extraction of instantaneous frequency and the characteristics of low frequency increase and high frequency attenuation existing in the oil and gas bearing formation, it is necessary to study a more reasonable fluid identification method with faster calculation. SUMMARY

[0003] In view of the above problems, the embodiments of the present application provide an oil and gas detection method and device based on instantaneous energy frequency value, a computer readable storage medium and an electronic device.

[0004] In a first aspect, the embodiments of the present application provide an oil and gas detection method based on instantaneous energy frequency value, comprising:

[0005] S100, acquiring a three-dimensional seismic data body of a target area;

[0006] S200, for each seismic signal in the three-dimensional seismic data body, performing the following steps S210 to S240:

[0007] S210, decompose the seismic signal into a data set composed of several IMF components;

[0008] S220, extract the instantaneous amplitude and the instantaneous frequency of each of the IMF components;

[0009] S230, calculate the instantaneous energy-frequency value of each of the IMF components based on the instantaneous amplitude and the instantaneous frequency;

[0010] S240, calculate the instantaneous energy-frequency value of the seismic signal based on the instantaneous energy-frequency value of each of the IMF components;

[0011] S300, analyze the variation characteristics of energy and frequency in the target area formation based on the instantaneous energy-frequency value of each seismic signal in the three-dimensional seismic data body and the instantaneous energy-frequency value of each IMF component of each seismic signal, and analyze the oil and gas bearing property of the formation according to the variation characteristics of energy and frequency.

[0012] According to the embodiment of the present application, in the step S210, the EMD is used to decompose the seismic signal into a data set composed of several IMF components.

[0013] According to the embodiment of the present application, in the step S220, the Hilbert transform is performed on each of the IMF components to extract the instantaneous amplitude and the instantaneous frequency of each of the IMF components.

[0014] According to the embodiment of the present application, in the step S220, the instantaneous energy-frequency value of the i-th IMF component is calculated based on the instantaneous amplitude and the instantaneous frequency according to the following formula: i

[0015]

[0016] In the formula, N is the time window length, InsAMP i is the instantaneous amplitude, InsFre i is the instantaneous frequency.

[0017] According to the embodiment of the present application, in the step S230, the instantaneous energy-frequency value of the seismic signal is calculated based on the instantaneous energy-frequency value of each of the IMF components according to the following formula:

[0018]

[0019] In the formula, n is the number of decomposed IMF components.

[0020] According to the embodiment of the present application, in the step S300, the instantaneous energy-frequency value of each IMF component is used for comparative analysis to carry out the oil and gas bearing property prediction.

[0021] According to the embodiment of the present application, the method further comprises the following steps:​

[0022] S400, identifying energy and frequency anomalies in the formation based on the instantaneous energy-frequency value of each seismic signal in the three-dimensional seismic data volume and the instantaneous energy-frequency value of each IMF component of each seismic signal.

[0023] In a second aspect, the embodiments of the present application also provide an oil and gas detection device based on instantaneous energy-frequency value, characterized in that it comprises:

[0024] An acquisition module is configured to acquire a three-dimensional seismic data volume of a target area.

[0025] A decomposition module is configured to decompose each seismic signal in the three-dimensional seismic data volume into a data set composed of a plurality of IMF components, extract the instantaneous amplitude and instantaneous frequency of each IMF component, calculate the instantaneous energy-frequency value of each IMF component based on the instantaneous amplitude and instantaneous frequency, and calculate the instantaneous energy-frequency value of the seismic signal based on the instantaneous energy-frequency value of each IMF component.

[0026] An analysis module is configured to analyze the change characteristics of energy and frequency in the formation of the target area based on the instantaneous energy-frequency value of each seismic signal in the three-dimensional seismic data volume and the instantaneous energy-frequency value of each IMF component of each seismic signal, and analyze the oil and gas bearing property of the formation according to the change characteristics of energy and frequency.

[0027] In a third aspect, the embodiments of the present application provide a computer readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the oil and gas detection method based on instantaneous energy-frequency value as described in the first aspect.

[0028] In a fourth aspect, the embodiments of the present application provide an electronic device comprising:

[0029] A processor;

[0030] A memory for storing instructions executable by the processor;

[0031] The processor is configured to execute the instructions to implement the oil and gas detection method based on instantaneous energy-frequency value as described in the first aspect.

[0032] Compared with the prior art, the above technical solutions of the present application have the following beneficial effects:

[0033] The embodiment of the present application provides an oil and gas detection method based on instantaneous energy frequency value, which is used for judging or detecting oil and gas characteristics of a reservoir and improving the fluid identification capability. BRIEF DESCRIPTION OF DRAWINGS

[0034] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or the prior art description. Obviously, the drawings in the following description only show some of the embodiments of the present application, and other drawings can be obtained by those of ordinary skill in the art without any creative effort based on these drawings.

[0035] Figure 1 A workflow diagram of the oil and gas detection method based on the instantaneous energy frequency value of the embodiment of the present application is shown.

[0036] Figure 2 A structural schematic diagram of an electronic device of the embodiment of the present application is shown. DETAILED DESCRIPTION

[0037] In order to make the purpose, technical solutions and advantages of the embodiments of the present application more clear, the following will combine the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without any creative effort are within the protection scope of the present application.

[0038] Embodiment one

[0039] After the reservoir contains oil and gas, the difference between the reservoir impedance and the surrounding rock impedance becomes larger, and the seismic reflection energy is stronger. After the seismic wave passes through the oil and gas containing formation, the high frequency energy is attenuated, and the main frequency is reduced. According to this feature, the present application combines the characteristics of the "sweet spot" attribute, empirical mode decomposition and instantaneous frequency to construct a new fluid identification attribute, namely the instantaneous energy frequency value. The instantaneous energy frequency value attribute extracted by the present application can more accurately identify the energy and frequency anomalies in the formation. Since the frequency band ranges of the IMF components obtained by the empirical mode decomposition are different, and the characteristics of the seismic wave passing through the oil and gas containing formation can cause high frequency attenuation and low frequency increase, the instantaneous energy frequency value of different IMF components can also be used for comparative analysis to carry out oil and gas prediction.

[0040] Based on the above concept, the present application provides a method for oil and gas detection based on instantaneous energy-frequency value, which is used for extracting energy-frequency value attribute of target area based on post-stack seismic signal of the target area, and then detecting oil and gas bearing property of the stratum according to the energy-frequency value attribute. Figure 1 As shown in the figure, the method specifically realizes steps including but not limited to the following steps:

[0041] S100, acquiring three-dimensional seismic data volume of target area;

[0042] S200, for each seismic signal in the three-dimensional seismic data volume, performing the following steps S210 to S240:

[0043] S210, decomposing the seismic signal into a data set composed of several IMF components;

[0044] S220, extracting instantaneous amplitude and instantaneous frequency of each IMF component;

[0045] S230, calculating instantaneous energy-frequency value of each IMF component based on the instantaneous amplitude and the instantaneous frequency;

[0046] S240, calculating instantaneous energy-frequency value of the seismic signal based on the instantaneous energy-frequency value of each IMF component;

[0047] S300, analyzing the change characteristics of energy and frequency in the stratum of the target area based on the instantaneous energy-frequency value of each seismic signal in the three-dimensional seismic data volume and the instantaneous energy-frequency value of each IMF component of each seismic signal, and analyzing the oil and gas bearing property of the stratum according to the change characteristics of energy and frequency.

[0048] The implementation process of each step is described in detail below.

[0049] S100, acquiring three-dimensional seismic data volume of target area;

[0050] S200, for each seismic signal in the three-dimensional seismic data volume, performing the following steps S210 to S240:

[0051] S210, decomposing the seismic signal into a data set composed of several IMF components;

[0052] For post-stack seismic signal Signal(t) of the target area, the post-stack seismic signal Signal(t) is decomposed by EMD, and the seismic signal Signal(t) is decomposed into a data set composed of several IMF components, i.e. Where IMF i (t) is the i-th IMF component, RES n (t) is the residual.

[0053] S220, extracting the instantaneous amplitude and the instantaneous frequency of each IMF component;

[0054] Specifically, the IMF i (t) is subjected to Hilbert transform to extract the instantaneous amplitude InsAMP i (t) of the IMF i (t) and the instantaneous frequency InsFre i (t) of the IMF

[0055] S230, calculating the instantaneous energy-frequency value of each IMF component based on the instantaneous amplitude and the instantaneous frequency;

[0056] The energy-frequency value IEF i (t) of the i-th IMF component is calculated according to formula (1):

[0057]

[0058] In the formula, N is the length of the time window.

[0059] S240, calculating the instantaneous energy-frequency value of the seismic signal based on the instantaneous energy-frequency value of each IMF component;

[0060] The instantaneous energy-frequency value IEF(t) of the seismic signal is calculated according to formula (2):

[0061]

[0062] In the formula, n is the number of decomposed IMF components.

[0063] S300: the instantaneous energy-frequency value data body IEF cube and the instantaneous energy-frequency value IEF cube (i) of each component of the entire three-dimensional seismic data body are obtained by repeating steps S210 to S240 for any seismic signal in the three-dimensional seismic data body. cube and the instantaneous energy-frequency value IEF cube (i) of each component of the entire three-dimensional seismic data body.

[0064] Embodiment Two

[0065] The implementation process and technical effects of the technical solution of the present application will be described below in combination with an embodiment.

[0066] In the present embodiment, an oil and gas detection method based on instantaneous energy-frequency value can be implemented according to the following steps:

[0067] Firstly, the IMF component of the seismic signal Signal(t) is extracted. Specifically, the following steps are included:

[0068] (11) Initialization, r0(t) = Signal(t), let i = 1;

[0069] (12) Screening the i-th IMF component, specifically including the following steps:

[0070] (12a) Initialization, h0(t) = r i-1 (t), j = 1;

[0071] (12b) According to the waveform of the signal h j-1 (t), the positions and amplitude values of the local maximum points and local minimum points of h j-1 (t) are found out, and a boundary processing method is used to extend the local maximum points and local minimum points to both sides by 1-2 extreme points respectively;

[0072] (12c) The local minimum points of h j-1 (t) are differentiated by using the difference method to obtain the lower envelope line of h j-1 (t); the local maximum points of h j-1 (t) are differentiated to obtain the upper envelope line of h j-1 (t);

[0073] (12d) The average value m j-1 (t) of the upper envelope line and the lower envelope line obtained in step (12c) is calculated;

[0074] (12e) Subtraction: h j (t) = h j-1 (t) - m j-1 (t);

[0075] (12f) If the termination condition is met, then IMF i (t) = h j (t) (this step searches for an IMF component, i.e. the i-th IMF component, if the termination condition is met, and the subsequent steps (13) and (14) search for the next IMF component until the condition of step (14) is met, and multiple IMF components can be searched from the signal), and then go to step (13); otherwise, j = j + 1, and go to (12b);

[0076] In the above steps, the boundary processing method can adopt mirror extension, neural network prediction or waveform matching method, etc. Different processing methods can achieve different technical effects by being adopted for different types of signals.

[0077] In the above step, the difference between the local minimum point and the local maximum point can be calculated by using cubic spline difference, linear combination of B-spline function, local filtering method, piecewise power function interpolation method or high-order spline fitting method.

[0078] In the above step, the termination condition can be that the Cauchy criterion is satisfied or the number of extreme points is equal to the number of zero points or at most one difference.

[0079] (13) Calculate the residual signal: r i (t) = r i-1 (t) - IMF i (t).

[0080] (14) Determine whether the number of extreme points of r i (t) is still more than two, if yes, i = i + 1, go to step (12), if no, the decomposition is finished, r i (t) is the residual component.

[0081] Then, the instantaneous amplitude InsAMP i (t) and the instantaneous frequency InsFre i (k) of each IMF component are extracted. Specifically, the following steps are included:

[0082] (21) Calculate the imaginary part IMF i (t) of IMF i (t) by using Hilbert transform. im .

[0083] (22) Calculate the instantaneous amplitude InsAMP i (t) of IMF i (t) by using the following formula:

[0084]

[0085] (23) Calculate the instantaneous frequency InsFre i (t) of IMF i (k) by using the following formula:

[0086]

[0087] The extraction of the imaginary part of the seismic signal by using Hilbert transform is a recognized algorithm.

[0088] Subsequently, the instantaneous energy frequency value IEF(i,j) of each IMF component is extracted.

[0089] The energy frequency value IEF(i,j) of the i-th IMF component is calculated according to formula (1). N is the length of the time window.

[0090]

[0091] Subsequently, the instantaneous energy frequency value IEF(t) of the seismic signal Signal(t) is calculated. The instantaneous energy frequency value IEF(t) of the seismic signal is obtained according to formula (2).

[0092]

[0093] In the formula, n is the number of decomposed IMF components.

[0094] Thus, by repeatedly performing the above steps on any seismic signal in the three-dimensional seismic data volume, the instantaneous energy frequency value data volume IEF cube and the instantaneous energy frequency value IEF cube (i) of each component of the entire three-dimensional seismic data volume can be extracted.

[0095] Then, based on the instantaneous energy frequency value data volume IEF cube and the instantaneous energy frequency value IEF cube (i) of each component of the entire three-dimensional seismic data volume, the variation characteristics of energy and frequency in the target area stratum are analyzed, and then the oil and gas bearing property of the stratum is analyzed according to the variation characteristics of energy and frequency.

[0096] In summary, the present application proposes an oil and gas detection method based on instantaneous energy frequency value. According to the characteristics that the difference between reservoir impedance and surrounding rock impedance becomes larger after the reservoir contains oil and gas, the seismic reflection energy is stronger; and after the seismic wave passes through the oil and gas bearing stratum, the high frequency energy attenuates, the main frequency decreases, etc., the present application combines the characteristics of "sweet spot" attribute, empirical mode decomposition and instantaneous frequency to construct a new fluid identification attribute-instantaneous energy frequency value. The instantaneous energy frequency value attribute extracted by the present application can more accurately identify the energy and frequency anomalies in the stratum. Since the frequency band ranges of the IMF components obtained by empirical mode decomposition are different, and the characteristics that high frequency attenuates and low frequency increases when the seismic wave passes through the oil and gas bearing stratum, the instantaneous energy frequency values of different IMF components can also be used for comparative analysis to carry out oil and gas bearing property prediction.

[0097] Example Three

[0098] The following is an embodiment of the device of the present application, which can be used to execute the method embodiment of the present application. For details not disclosed in the device embodiment of the present application, please refer to the method embodiment of the present application.

[0099] The present embodiment provides an oil and gas detection device based on instantaneous energy frequency value, characterized in that it comprises:

[0100] An acquisition module is configured to acquire a three-dimensional seismic data volume of a target area.

[0101] a decomposition module, configured to decompose each seismic signal in the three-dimensional seismic data volume into a data set composed of a plurality of IMF components, extract the instantaneous amplitude and the instantaneous frequency of each IMF component, calculate the instantaneous energy-frequency value of each IMF component based on the instantaneous amplitude and the instantaneous frequency, calculate the instantaneous energy-frequency value of the seismic signal based on the instantaneous energy-frequency values of the IMF components of the seismic signal,

[0102] an analysis module, configured to analyze the variation characteristics of energy and frequency in the target area formation based on the instantaneous energy-frequency value of each seismic signal in the three-dimensional seismic data volume and the instantaneous energy-frequency value of each IMF component of each seismic signal, and analyze the oil and gas bearing property of the formation according to the variation characteristics of energy and frequency.

[0103] Embodiment four

[0104] The embodiment provides a computer readable medium, which stores a computer program, and the program is executed by a processor to implement each step of the oil and gas detection method based on instantaneous energy-frequency value.

[0105] It should be noted that the present application can implement all or part of the processes of the above-mentioned embodiments, which can be completed by a computer program instructing related hardware. The computer program can be stored in a computer readable storage medium, and the computer program can implement the steps of the above-mentioned embodiments when executed by a processor. The computer program includes computer program code, which can be in the form of source code, object code, executable files or some intermediate forms. The computer readable medium can include any entity or device, recording medium, U disk, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signal, telecommunication signal and software distribution medium, etc. that can carry the computer program code. Of course, there are other ways of readable storage medium, such as quantum memory, graphene memory, etc. It should be noted that the content of the computer readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction, for example, in some jurisdictions, according to legislation and patent practice, the computer readable medium does not include electrical carrier signals and telecommunication signals.

[0106] Embodiment five

[0107] Figure 2 is a structural schematic diagram of an electronic device according to an embodiment of the present application. As shown in Figure 2As shown in the hardware layer, the electronic device includes a processor, and optionally further includes an internal bus, a network interface, and a memory. The memory can include a memory such as a random-access memory (RAM), and can further include a non-volatile memory such as at least one disk memory. Of course, the electronic device can further include other hardware required by a business.

[0108] The processor, the network interface, and the memory can be connected to each other through the internal bus, which can be an industry standard architecture (ISA) bus, a peripheral component interconnect (PCI) bus, or an extended industry standard architecture (EISA) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, only a line segment is shown in the figure, but it does not mean that there is only one bus or only one type of bus.

[0109] The memory is used to store programs. Specifically, the program can include program code including computer operation instructions. The memory can include a memory and a non-volatile memory, and provide instructions and data to the processor. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs. The processor executes the program stored in the memory to perform all the steps of the aforementioned oil and gas detection method based on the instantaneous energy frequency value.

[0110] The communication bus mentioned in the above device can be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus, etc. The communication bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, only one thick line is shown in the figure, but it does not mean that there is only one bus or only one type of bus. The communication interface is used for communication between the above electronic device and other devices.

[0111] The bus includes hardware, software, or both, for coupling together various components within the system. For example, the bus can include an Accelerated Graphics Port (AGP) or other graphics bus, an Enhanced Industry Standard Architecture (EISA) bus, a Front Side Bus (FSB), a HyperTransport (HT) interconnect, an Industry Standard Architecture (ISA) bus, an Infmiband interconnect, a Low Pin Count (LPC) bus, a memory bus, a Micro Channel Architecture (MCA) bus, a Peripheral Component Interconnect (PCI) bus, a PCI-Express (PCI-X) bus, a Serial Advanced Technology Attachment (SATA) bus, a Video Electronics Standards Association local (VLB) bus, or another suitable bus or a combination of two or more of these. Where appropriate, the bus can include one or more buses. Although the present embodiments describe and show a particular bus, the present embodiments contemplate any suitable bus or interconnect.

[0112] The memory can include a Random Access Memory (RAM) or can include non-volatile memory (NVM), e.g., at least one disk memory. Optionally, the memory can also be at least one storage device located remotely from the aforementioned processor.

[0113] The memory can include mass storage for data or instructions. As an example and not by way of limitation, the memory can include a Hard Disk Drive (HDD), a floppy disk drive, flash memory, a Compact Disc (CD) or other optical disk, a Universal Serial Bus (USB) drive or other memory chip or cartridge, or a combination of two or more of these. Where appropriate, the memory can include removable or non-removable (or fixed) media. In particular embodiments, the memory is non-volatile, solid-state memory. In particular embodiments, the memory includes Read-Only Memory (ROM). Where appropriate, this ROM can be mask-programmed ROM, programmable ROM (PROM), erasable PROM (EPROM), electrically erasable PROM (EEPROM), electrically alterable ROM (EAROM), or flash memory or a combination of two or more of these. In particular embodiments, the memory includes ROM. The memory can be static, dynamic or both, as appropriate.

[0114] The processor described above can be a general processor, including a central processing unit (CPU), a network processor (NP), etc.; can also be a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field-programmable gate array (FPGA) or other programmable logic device, a discrete gate or transistor logic device, a discrete hardware component.

[0115] It should be explained that, for the convenience and brevity of description, only the division of the above functional units and modules is exemplified, and in actual application, the above functions can be completed by different functional units and modules according to needs, that is, the internal structure of the apparatus is divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiments can be integrated in one processing unit, or each unit can exist physically, or two or more units can be integrated in one unit, and the integrated unit can be realized in the form of hardware or in the form of software functional unit. In addition, the specific names of the functional units and modules are only for easy distinction, and do not limit the protection scope of the present application. The specific working process of the units and modules in the system can refer to the corresponding process in the foregoing method embodiments, which will not be repeated here.

[0116] The apparatus, device, system, module or unit illustrated in the above embodiments can be specifically realized by a computer chip or entity, or by a product with certain functions. A typical implementation device is a computer. Specifically, the computer can be, for example, a personal computer, a laptop computer, a vehicle-mounted human-computer interaction device, a cellular phone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an electronic mail device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

[0117] Although the present application provides method operation steps as described in the embodiments or flowcharts, more or fewer operation steps can be included based on conventional or non-inventive means. The order of steps listed in the embodiments is only one of the many step execution orders, and does not represent the only execution order. In actual device or terminal product execution, the method order shown in the embodiments or the drawings can be executed in sequence or in parallel (for example, in a parallel processor or multi-thread processing environment, or even in a distributed data processing environment).

[0118] The present application is described in reference to the flowchart and / or block diagrams of the methods, apparatus (systems) and computer program products according to embodiments of the present application. It will be understood that each block of the flowchart and / or block diagrams, and combinations of blocks in the flowchart and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can 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 and / or block diagram block or blocks. Figure 1 Figure 1

[0119] These computer program instructions can 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 instructions which implement the function specified in the flowchart and / or block diagram block or blocks. Figure 1 Figure 1

[0120] These computer program instructions can 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 and / or block diagram block or blocks. Figure 1 Figure 1

[0121] It should be noted that, in the description, relative terms such as first and second, and the like, can be used solely to distinguish one from another without necessarily implying any actual relationship or order between or among the entities or operations so distinguished. Also, the terms "comprises", "comprising", or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but can include other elements not expressly listed or inherent to such process, method, article, or apparatus. An element proceeded by "comprises... a" does not, without more constraints, exclude the presence of additional identical elements in the process, method, article, or apparatus that comprises the element.

[0122] ​​​​​​Various embodiments are described in related manner in the specification, and the same or similar parts among various embodiments can be referred to each other. Each embodiment focuses on the difference from other embodiments. In particular, for the device, electronic device and readable storage medium embodiments, since they are basically similar to the method embodiments, the description is relatively simple, and the relevant parts can be referred to the part of the method embodiment.

[0123] The above only describes the preferred embodiments of the present application, and is not intended to limit the protection scope of the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for detecting oil and gas based on instantaneous energy frequency values, characterized in that, Includes the following steps: S100, acquire the three-dimensional seismic data volume of the target area; S200, for each seismic signal in the three-dimensional seismic data volume, perform the following steps S210 to S240: S210 uses EMD to decompose the seismic signal into a data set consisting of several IMF components; S220, extract the instantaneous amplitude and instantaneous frequency of each IMF component; S230, Based on the instantaneous amplitude and instantaneous frequency, calculate the instantaneous energy frequency value of each of the IMF components; S240, Calculate the instantaneous energy frequency value of the seismic signal based on the instantaneous energy frequency value of each IMF component; S300, based on the instantaneous energy-frequency value of each seismic signal and the instantaneous energy-frequency value of each IMF component of each seismic signal in the three-dimensional seismic data volume, analyze the characteristics of energy and frequency variation in the strata of the target area, and analyze the oil and gas content of the strata according to the characteristics of energy and frequency variation; Wherein, the instantaneous energy frequency value (IEF) of the i-th IMF component i for: In the formula, N is the time window length, InsAMP i For instantaneous amplitude, InsFre i Instantaneous frequency; The instantaneous energy frequency value of the seismic signal is: In the formula, n is the number of IMF components in the decomposition.

2. The oil and gas detection method based on instantaneous energy frequency value as described in claim 1, characterized in that, In step S220, a Hilbert transform is performed on each of the IMF components to extract the instantaneous amplitude and instantaneous frequency of each IMF component.

3. The oil and gas detection method based on instantaneous energy frequency value as described in claim 1, characterized in that, In step S300, the instantaneous energy-frequency values ​​of each IMF component are compared and analyzed to predict the oil and gas content.

4. The oil and gas detection method based on instantaneous energy frequency value as described in claim 1, characterized in that, It also includes the following steps: S400, based on the instantaneous energy-frequency value of each seismic signal in the three-dimensional seismic data volume and the instantaneous energy-frequency value of each IMF component of each seismic signal, identifies energy and frequency anomalies in the formation.

5. An oil and gas detection device based on instantaneous energy frequency value, characterized in that, include: The acquisition module is used to acquire the three-dimensional seismic data volume of the target area; The decomposition module is used to decompose each seismic signal in the three-dimensional seismic data volume into a data set composed of several IMF components using EMD, extract the instantaneous amplitude and instantaneous frequency of each IMF component, calculate the instantaneous energy-frequency value of each IMF component based on the instantaneous amplitude and instantaneous frequency, and calculate the instantaneous energy-frequency value of the seismic signal based on the instantaneous energy-frequency value of each IMF component. The analysis module is used to analyze the energy and frequency variation characteristics in the strata of the target area based on the instantaneous energy and frequency values ​​of each seismic signal and each IMF component of each seismic signal in the three-dimensional seismic data volume, and to analyze the oil and gas content of the strata based on the energy and frequency variation characteristics. Wherein, the instantaneous energy frequency value (IEF) of the i-th IMF component i for: In the formula, N is the time window length, InsAMP i For instantaneous amplitude, InsFre i Instantaneous frequency; The instantaneous energy frequency value of the seismic signal is: In the formula, n is the number of IMF components in the decomposition.

6. A computer-readable storage medium, characterized in that, It stores a computer program, which, when executed by a processor, implements an oil and gas detection method based on instantaneous energy frequency value as described in any one of claims 1 to 4.

7. An electronic device comprising: processor; Memory used to store the processor's executable instructions; The processor is configured to execute the instructions to implement an oil and gas detection method based on instantaneous energy frequency value as described in any one of claims 1 to 4.

Images