Volcano channel identification method and device, electronic equipment and storage medium

Abstract

The invention discloses a volcano channel identification method and device, electronic equipment and a storage medium, and the method comprises the steps: carrying out the ACF processing of seismic data, so as to enhance the structural features in the seismic data; performing bilinear interpolation calculation on the seismic data subjected to ACF processing to obtain effective information of the seismic data; and regional difference calculation is carried out based on the seismic data subjected to ACF processing and bilinear interpolation, so that weak seismic signals are enhanced, and volcano channel identification is realized. According to the method, the volcano channel can be quickly and effectively identified.

Description

Technical Field

[0001] The present invention belongs to the field of geophysical exploration, and more specifically, relates to a method and device for identifying volcanic conduits, an electronic device, and a storage medium. Background Art

[0002] Volcanic rock eruption facies are well-developed in basins and are one of the main gas-producing strata of volcanic rocks. Currently, there is little research on volcanic conduits. Volcanic conduits play a role in connecting the upper and lower parts, and their response characteristics on seismic profiles are not obvious, which restricts the effective characterization of volcanic conduits. This has led to insufficient research on volcanic structures and understanding of favorable reservoirs. Through conventional geometric attributes in the past, such as high-precision coherence and other single attributes, there are large ambiguities and the recognition results are not clear. And through manual interpretation, not only is the workload large, but there are also certain interference factors, making the recognition results uncertain. Summary of the Invention

[0003] The object of the present invention is to propose a method and device for identifying volcanic conduits, an electronic device, and a storage medium to achieve rapid and effective identification of volcanic conduits.

[0004] To achieve the above object, in a first aspect, the present invention proposes a method for identifying volcanic conduits, including:

[0005] Performing ACF processing on seismic data to enhance the structural features in the seismic data;

[0006] Performing bilinear interpolation calculation on the seismic data after ACF processing to obtain effective information of the seismic data;

[0007] Based on the seismic data after ACF processing and bilinear interpolation, performing regional difference calculation to enhance weak seismic signals to achieve the identification of volcanic conduits.

[0008] Optionally, it further includes:

[0009] Judging whether the identified volcanic conduit is clear. If not, performing normalization processing on the seismic data after completing the regional difference calculation to obtain a clearer volcanic conduit.

[0010] Optionally, the calculation formula for performing ACF processing on the seismic data is:

[0011]

[0012] where I is the seismic data, t is the number of iterations, x and y are the horizontal and vertical coordinates of the pixel, is the divergence of the seismic data I in four directions, cN t cS xy cE xy cW xy cS xy cWxy represent the thermal conductivities in four directions, and ( is the smoothing correlation coefficient.

[0013] Optionally, the formula for performing bilinear interpolation calculation on the ACF-processed seismic data is:

[0014]

[0015]

[0016]

[0017] where I(R1) is the linear interpolation result of the inserted point R1 = (x, y 1 , y 1 ) between the pixel point (x 2 , y 1 ) and the pixel point (x 1 ), y 1 , y 2 ) and the pixel point (x 2 , y 2 ), and I(R2) is the linear interpolation result of the inserted point R2 = (x, y 2 ) between the pixel point (x

[0018] Optionally, the formula for performing regional difference calculation based on the ACF-processed and bilinearly interpolated seismic data is:

[0019]

[0020] where

[0021] where Var(x, y) is the local variance, (x + i, y + j) represents the pixel of the coordinate point centered on (x, y), and W s represents the number of pixels within the window, and M(x) represents the mean within the local window.

[0022] In a second aspect, the present invention provides an electronic device, which includes:

[0023] at least one processor; and,

[0024] a memory communicatively connected to the at least one processor; wherein,

[0025] the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute any one of the volcanic channel recognition methods in the first aspect.

[0026] In a third aspect, the present invention provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the volcanic channel recognition method according to any one of the first aspects.

[0027] In a fourth aspect, the present invention provides a volcanic channel recognition device, comprising:

[0028] an ACF processing module configured to perform ACF processing on seismic data to enhance the structural features in the seismic data;

[0029] a bilinear interpolation calculation module configured to perform bilinear interpolation calculation on the seismic data after ACF processing to obtain effective information of the seismic data;

[0030] a regional difference calculation module configured to perform regional difference calculation based on the seismic data after ACF processing and bilinear interpolation to enhance weak seismic signals for realizing the recognition of volcanic channels.

[0031] Optionally, the device further comprises a normalization processing module configured to determine whether the recognized volcanic channel is clear. If not, the normalization processing module performs normalization processing on the seismic data after the regional difference calculation to obtain a clearer volcanic channel.

[0032] Optionally, the calculation formula for performing ACF processing on the seismic data is as follows:

[0033]

[0034] where I is the seismic data, t is the number of iterations, x and y are the horizontal and vertical coordinates of the pixel, is the divergence of the seismic data I in four directions, cN t xy xy xy xy represent the thermal conductivities in four directions, and λ is the smoothing correlation coefficient.

[0035] The beneficial effects of the present invention are as follows:

[0036] The present invention first performs ACF processing on seismic data to enhance the structural features in the seismic data, then performs bilinear interpolation calculation on the seismic data after ACF processing to obtain effective information of the seismic data, and then performs regional difference calculation based on the seismic data after ACF processing and bilinear interpolation to enhance weak seismic signals, thereby realizing fast and accurate recognition of volcanic channels.

[0037] ​​​​The system of the present invention has other features and advantages, which will be apparent from the accompanying drawings incorporated herein and the subsequent detailed description, or will be described in detail in the accompanying drawings incorporated herein and the subsequent detailed description, and these drawings and detailed description together are used to explain the specific principles of the present invention. Description of the Drawings

[0038] By describing the exemplary embodiments of the present invention in more detail in conjunction with the accompanying drawings, the above and other objects, features, and advantages of the present invention will become more apparent. In the exemplary embodiments of the present invention, the same reference numerals generally represent the same components.

[0039] Figure 1 A step diagram showing a method for identifying a volcanic conduit according to the present invention is shown. Detailed Description of the Embodiments

[0040] As an important channel for the development of volcanic rock reservoirs, the seismic response characteristics of volcanic conduits are not obvious, and they cannot be identified through simple post-stack seismic attributes, resulting in the inability to further understand the reservoirs. At the same time, the workload of manual interpretation is large and there are subjective cognitive differences.

[0041] Therefore, the present invention proposes a method and device for identifying volcanic conduits, an electronic device, and a storage medium to achieve rapid and effective identification of volcanic conduits.

[0042] The present invention will be described in more detail below with reference to the accompanying drawings. Although the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to make the present invention more thorough and complete, and to fully convey the scope of the present invention to those skilled in the art.

[0043] Embodiment 1

[0044] As Figure 1 shown, this embodiment provides a method for identifying a volcanic conduit, including:

[0045] S1: Perform ACF processing on seismic data to enhance the structural features in the seismic data;

[0046] Specifically, the calculation formula for performing ACF processing on seismic data is:

[0047]

[0048] where I is the seismic data, t is the number of iterations, x and y are the horizontal and vertical coordinates of the pixel, is the divergence of the seismic data I in four directions, cN t ​xy , cS xy , cE xy , cW xy represent the thermal conductivities in four directions, I t = I x,y , then:

[0049]

[0050]

[0051]

[0052]

[0053] k and λ are smoothing correlation coefficients, where the larger the value of k, the smoother it is, and the less likely it is to retain edges; the same is true for λ, the larger the value, the smoother it is.

[0054] S2: Perform bilinear interpolation calculation on the seismic data processed by ACF to obtain the effective information of the seismic data;

[0055] Specifically, perform bilinear interpolation calculation on the data processed by ACF to obtain the effective information of the data, providing better data support for the subsequent identification of volcanic channels.

[0056] The formula for performing bilinear interpolation calculation on the seismic data processed by ACF is:

[0057]

[0058]

[0059]

[0060] where I(R1) is the linear interpolation result of the inserted point R1 = (x, y 1 , y 1 ) between the known pixel points (x 2 , y 1 ) and the known pixel points (x 1 ), y 1 , y 2 ) and the known pixel points (x 2 , y 2 ), that is, the pixel value of R1, I(R2) is the linear interpolation result of the inserted point R2 = (x, y 2 ) between the known pixel points (x 1

[0061]

[0061] 2 ) and the known pixel points (x 2 , y 2 ), that is, the pixel value of R2, and I(P) is the linear interpolation result of the inserted point P = (x, y) between R1 and R2 in the y direction, that is, the pixel value of the inserted point P to be found.

[0061] S3: Calculate the regional difference based on the seismic data processed by ACF and bilinear interpolation to enhance weak seismic signals for volcanic channel identification.

[0062] Specifically, the formula for calculating the regional difference based on the seismic data processed by ACF and bilinear interpolation is as follows:

[0063]

[0064] Among them,

[0065] Among them, Var(x, y) is the local variance, (x + i, y + j) represents the pixel of the coordinate point centered on (x, y), and W s represents the number of pixels within the window, and M(x) represents the mean within the local window.

[0066] In this embodiment, preferably, after step S3, it further includes:

[0067] S4: Determine whether the identified volcanic channel is clear. If not, perform normalization processing on the seismic data after the regional difference calculation to obtain a clearer volcanic channel.

[0068] Specifically, observe the volcanic channel identified through the regional difference calculation. If the volcanic channel is not prominent, the data can be further normalized to obtain a clearer volcanic channel.

[0069] Embodiment 2

[0070] This embodiment is described by comparing the identification results of using conventional attributes to identify volcanic channels with those of the volcanic channel identification method in Embodiment 1:

[0071] Extract conventional seismic attributes from the seismic data and select the coherence attribute to identify volcanic channels.

[0072] Using the same seismic data, first perform ACF processing, perform bilinear interpolation on the processed data, and then calculate the regional variance of the interpolated data to obtain the identification result of the volcanic channel.

[0073] Compare the two identification results:

[0074] The volcanic channels identified through the coherence attribute have multiple solutions and are not clearly identified. However, the volcanic channel identification results obtained through the method of Embodiment 1 can well identify and depict the volcanic channels and have good consistency with geological understanding.

[0075] Embodiment 3

[0076] This embodiment provides a volcanic channel identification device, including:

[0077] An ACF processing module, configured to perform ACF processing on seismic data to enhance the structural features in the seismic data;

[0078] A bilinear interpolation calculation module, configured to perform bilinear interpolation calculation on the seismic data after ACF processing to obtain the effective information of the seismic data;

[0079] A regional difference calculation module, configured to perform regional difference calculation based on the seismic data after ACF processing and bilinear interpolation, so as to enhance weak seismic signals to achieve the identification of volcanic channels.

[0080] Preferably, the device further includes a normalization processing module, configured to determine whether the identified volcanic channel is clear. If not, perform normalization processing on the seismic data after completing the regional difference calculation to obtain a clearer volcanic channel.

[0081] Among them, the calculation formula for performing ACF processing on seismic data is:

[0082]

[0083] Among them, I is the seismic data, t is the number of iterations, x and y are the horizontal and vertical coordinates of the pixel, is the divergence of the seismic data I in four directions, cN t is the divergence in four directions, cN xy , cS xy , cE xy , cW xy represent the thermal conductivities in four directions, and λ is the smoothing correlation coefficient.

[0084] Embodiment 4

[0085] This embodiment provides an electronic device, and the electronic device includes:

[0086] At least one processor; and,

[0087] A memory communicatively connected to the at least one processor; wherein,

[0088] The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor, so that the at least one processor can execute the volcanic channel identification method described in Embodiment 1.

[0089] An electronic device according to an embodiment of the present disclosure includes a memory and a processor. The memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, and the computer program products may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory, etc. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc.

[0090] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In an embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.

[0091] Those skilled in the art should understand that, in order to solve the technical problem of how to obtain good user experience effects, this embodiment may also include well-known structures such as communication buses, interfaces, etc., and these well-known structures should also be included in the protection scope of the present disclosure.

[0092] For the detailed description of this embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, and details will not be repeated here.

[0093] Embodiment 5

[0094] This embodiment provides a non-transitory computer-readable storage medium. The non-transitory computer-readable storage medium stores computer instructions, and the computer instructions are used to cause a computer to execute the volcanic channel recognition method described in Embodiment 1.

[0095] A computer-readable storage medium according to an embodiment of the present disclosure has non-transitory computer-readable instructions stored thereon. When the non-transitory computer-readable instructions are run by a processor, all or part of the steps of the methods of the foregoing embodiments of the present disclosure are executed.

[0096] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (such as: CD-ROM and DVD), magneto-optical storage media (such as: MO), magnetic storage media (such as: magnetic tape or removable hard disk), media with built-in rewritable non-volatile memory (such as: memory card) and media with built-in ROM (such as: ROM cartridge).

[0097] The various embodiments of the present invention have been described above. The above description is exemplary and not exhaustive, and is not limited to the disclosed embodiments. Many modifications and variations are obvious to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for identifying volcanic channels, characterized in that, it includes: Performing ACF processing on seismic data to enhance the structural features in the seismic data; Performing bilinear interpolation calculation on the seismic data after ACF processing to obtain the effective information of the seismic data; Performing regional difference calculation based on the seismic data after ACF processing and bilinear interpolation to enhance weak seismic signals, so as to realize the identification of volcanic channels.

2. The method for identifying volcanic channels according to claim 1, characterized in that, it further includes: Judging whether the identified volcanic channel is clear. If not, performing normalization processing on the seismic data after completing the regional difference calculation to obtain a clearer volcanic channel.

3. The method for identifying volcanic channels according to claim 1, characterized in that, The calculation formula for performing ACF processing on the seismic data is: Where I is the seismic data, t is the number of iterations, and x and y are the horizontal and vertical coordinates of the pixel. is the seismic data I t The divergence in four directions, cN xy , cS xy , cE xy , cW xy represent the thermal conductivities in four directions, and λ is the smoothing correlation coefficient.

4. The method for identifying volcanic channels according to claim 3, characterized in that, The calculation formula for performing bilinear interpolation calculation on the seismic data after ACF processing is: Among them, I(R1) is the linear interpolation result of the point R1=(x,y 1 ,y 1 ) inserted between the pixel point (x 2 ,y 1 ) and the pixel point (x 1 ), I(R2) is the linear interpolation result of the point R2=(x,y 1 ,y 2 ) inserted between the pixel point (x 2 ,y 2 ) and the pixel point (x 2 ), and I(P) is the linear interpolation result of the point P=(x, y) inserted between R1 and R2 in the y direction.

5. The method for identifying volcanic channels according to claim 4, characterized in that, The calculation formula for performing regional difference calculation based on the seismic data after ACF processing and bilinear interpolation is: Among them, Among them, Var(x, y) is the local variance, (x + i, y + j) represents the pixel of the coordinate point centered on (x, y), and W s represents the number of pixels within the window, and M(x) represents the mean value within the local window.

6. An electronic device, characterized in that, the electronic device includes: At least one processor; and, A memory communicatively connected to the at least one processor; wherein, The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the method for identifying volcanic channels according to any one of claims 1-5.

7. A non-transitory computer-readable storage medium, characterized in that, This non-transitory computer-readable storage medium stores computer instructions, and these computer instructions are used to make a computer execute the method for identifying volcanic channels according to any one of claims 1-5.

8. A volcanic channel identification device, characterized in that, it includes: An ACF processing module for performing ACF processing on seismic data to enhance the structural features in the seismic data; A bilinear interpolation calculation module for performing bilinear interpolation calculation on the seismic data after ACF processing to obtain the effective information of the seismic data; A regional difference calculation module for performing regional difference calculation based on the seismic data after ACF processing and bilinear interpolation to enhance weak seismic signals, so as to realize the identification of volcanic channels.

9. The volcanic channel identification device according to claim 8, characterized in that, It further includes a normalization processing module for judging whether the identified volcanic channel is clear. If not, performing normalization processing on the seismic data after completing the regional difference calculation to obtain a clearer volcanic channel.

10. The volcanic channel identification device according to claim 8, characterized in that, The calculation formula for performing ACF processing on the seismic data is: Among them, I is seismic data, t is the number of iterations, and x and y are the horizontal and vertical coordinates of the pixels. is the seismic data I t the divergence in four directions, cN xy , cS xy , cE xy , cW xy represent the thermal conductivities in four directions, and ( is the smoothing correlation coefficient.

Images