Earthquake coherent body calculation method and device, electronic equipment and readable storage medium

Abstract

The invention provides a seismic coherent body calculation method and device, electronic equipment and a computer readable storage medium. The method comprises the following steps: sequentially extracting three-dimensional sampling data volumes from three-dimensional post-stack seismic data volumes of a to-be-processed region by taking a single data sampling point as a center according to a preset calculation window; determining a local event slope at the central position of each three-dimensional sampling data volume by using a direct search method; and according to the local event slope at the central position of each three-dimensional sampling data volume and a coherence calculation formula based on a variational principle, determining a coherence value at the central position of each three-dimensional sampling data volume, and obtaining a seismic coherence body result of the to-be-processed region. Based on the variational principle, the adaptive standard trace is introduced to calculate the seismic coherent body, so that the influence of the transverse change of the amplitude of seismic data on the calculation of the coherent body can be eliminated, and the resolution capability of discontinuous bodies such as faults and unconformity surfaces can be remarkably improved.

Description

Technical Field

[0001] The present invention relates to the field of geophysical technology, and in particular to a method, device, electronic device and computer-readable storage medium for calculating seismic coherence volumes. Background Art

[0002] In the field of exploration seismology, discontinuity characterization is a key issue in seismic data interpretation. Especially when accurately interpreting faults, reservoirs, and oil and gas migration processes, fast and high-resolution discontinuity characterization tools will significantly improve the efficiency and accuracy of seismic interpretation. Coherence volumes, a key seismic attribute, are currently one of the most important tools for discontinuity characterization. The results of coherence volume calculations can be used not only for fault interpretation in exploration areas but also for monitoring the quality of seismic data processing.

[0003] In the prior art, there are three main generations of coherence volume calculation methods. The first generation coherence volume algorithm (hereinafter referred to as the C1 algorithm) is a cross-correlation-based coherence volume algorithm. This method calculates the cross-correlation between a selected central trace and two adjacent traces, and uses the cross-correlation result to measure the coherence value. The second generation coherence volume algorithm (hereinafter referred to as the C2 algorithm) is a coherence volume algorithm based on multi-trace similarity. This algorithm selects a standard trace and calculates the sum of the least squares errors between each seismic trace within the selected window and the standard trace, and uses this error sum to measure the coherence value. The third generation coherence volume algorithm (hereinafter referred to as the C3 algorithm) is a coherence volume algorithm based on characteristic structure. This algorithm measures the coherence value of the corresponding seismic data by calculating the maximum eigenvalue of the covariance matrix of the seismic data within a specific window.

[0004] While the C1 algorithm performs well with high-quality data, it suffers from poor noise immunity when faced with seismic data with a low signal-to-noise ratio, often failing to accurately distinguish discontinuities. The C2 algorithm can produce erroneous results when faced with seismic data with dramatic lateral amplitude variations. These erroneous results can reduce resolution and even lead to misjudgments in seismic interpretation. While the C3 algorithm offers higher resolution, its computational efficiency is low. Typically, the C3 algorithm takes ten times as long to compute the same data set as the C2 algorithm, a shortcoming that hinders its practical application. Summary of the Invention

[0005] The present application provides a method, device, electronic device and computer-readable storage medium for calculating a seismic coherence volume. The calculation method is improved on the basis of the C2 algorithm and can improve the resolution while maintaining the corresponding calculation efficiency.

[0006] In a first aspect, the present application discloses a method for calculating a seismic coherence volume, comprising:

[0007] According to the preset calculation window, with a single data sampling point as the center, three-dimensional sampling data volumes are sequentially extracted from the three-dimensional post-stack seismic data volume of the area to be processed;

[0008] Determining the local event slope at the center position of each three-dimensional sampling data volume using a direct search method;

[0009] Determining the coherence value at the center of each three-dimensional sampling data volume according to the local event slope at the center of each three-dimensional sampling data volume and a coherence calculation formula based on the variational principle, to obtain a seismic coherence volume result of the area to be processed;

[0010] The coherence calculation formula based on the variational principle is:

[0011]

[0012] Where: C v represents the coherence value based on the variational principle, u represents the three-dimensional sampling data volume, t represents the sampling time at the center position of the three-dimensional sampling data volume, △t is the time sampling interval, represents the local event slope at the center of the 3D sampling data volume, M represents the number of time samples in the 3D sampling data volume, i represents the index of the time sampling point, N represents the total number of channels in the 3D sampling data volume, j represents the index of the channel number, x and y represent the spatial coordinates of the seismic channel relative to the center of the 3D sampling data volume, Represents the average of all trace data along the spatial direction of the three-dimensional sampling data volume, a(j) represents the adaptive coefficient based on the variational principle,

[0013] In a possible implementation, the process of determining the local event slope at the center position of each three-dimensional sampling data volume using a direct search method includes:

[0014] Determining a plurality of slope sampling points according to a preset value range of the local event slope in the first direction and the second direction and a sampling interval;

[0015] At the center position of each of the three-dimensional sampling data volumes, a coherence value corresponding to each slope sampling point is calculated using a coherence calculation formula based on multi-channel similarity;

[0016] A slope sampling point corresponding to a minimum coherence value is selected at the center position of each of the three-dimensional sampling data volumes as the local event slope at the center position of each of the three-dimensional sampling data volumes.

[0017] In a possible implementation, the calculation window is a sampling range consisting of a preset number of seismic traces and a preset length of time sampling points around a single data sampling point.

[0018] In a possible implementation, the number of seismic traces in the calculation window is 3, 5, or 9.

[0019] In a second aspect, the present application further discloses a seismic coherence volume calculation device, comprising:

[0020] The data volume selection module is used to extract three-dimensional sampling data volumes from the three-dimensional post-stack seismic data volume of the area to be processed in sequence based on a preset calculation window and with a single data sampling point as the center;

[0021] a slope calculation module, configured to determine the local event slope at the center position of each of the three-dimensional sampling data volumes by using a direct search method;

[0022] a coherence value calculation module, configured to determine the coherence value at the center of each three-dimensional sampling data volume according to the local event slope at the center of each three-dimensional sampling data volume and a coherence calculation formula based on the variational principle, thereby obtaining a seismic coherence volume result of the area to be processed;

[0023] The coherence calculation formula based on the variational principle is:

[0024]

[0025] Where: C v represents the coherence value based on the variational principle, u represents the three-dimensional sampling data volume, t represents the sampling time at the center position of the three-dimensional sampling data volume, △t is the time sampling interval, represents the local event slope at the center of the 3D sampling data volume, M represents the number of time samples in the 3D sampling data volume, i represents the index of the time sampling point, N represents the total number of channels in the 3D sampling data volume, j represents the index of the channel number, x and y represent the spatial coordinates of the seismic channel relative to the center of the 3D sampling data volume, Represents the average of all trace data along the spatial direction of the three-dimensional sampling data volume, a(j) represents the adaptive coefficient based on the variational principle,

[0026] In a possible implementation, the slope calculation module includes:

[0027] A slope sampling point determination unit, configured to determine a plurality of slope sampling points according to a preset value range of the local event slope in the first direction and the second direction and a sampling interval;

[0028] a slope sampling point coherence value calculation unit, configured to calculate the coherence value corresponding to each slope sampling point at the center position of each three-dimensional sampling data volume using a coherence calculation formula based on multi-channel similarity;

[0029] The slope determining unit is configured to select the local event slope of the slope sampling point corresponding to the minimum coherence value at the center position of each of the three-dimensional sampling data volumes as the local event slope at the center position of each of the three-dimensional sampling data volumes.

[0030] In a third aspect, the present application further provides an electronic device, comprising: a memory, a processor;

[0031] The memory stores computer-executable instructions;

[0032] The processor executes the computer-executable instructions stored in the memory, so that the processor executes the above first aspect and / or various possible implementations of the first aspect.

[0033] In a fourth aspect, the present application also provides a computer-readable storage medium, in which computer-executable instructions are stored. When the computer-executable instructions are executed by a processor, they are used to implement the above first aspect and / or various possible implementations of the first aspect.

[0034] The seismic coherence volume calculation method, apparatus, device, and computer-readable storage medium provided in the embodiments of the present application sequentially extract three-dimensional sampling data volumes from a three-dimensional post-stack seismic data volume of a to-be-processed area based on a preset calculation window and with a single data sampling point as the center; determine the local event slope at the center position of each three-dimensional sampling data volume using a direct search method; and determine the coherence value at the center position of each three-dimensional sampling data volume based on the local event slope at the center position of each three-dimensional sampling data volume and a coherence calculation formula based on the variational principle to obtain a seismic coherence volume result for the to-be-processed area. The coherence calculation formula is: Where: C v represents the coherence value based on the variational principle, u represents the three-dimensional sampling data volume, t represents the sampling time at the center position of the three-dimensional sampling data volume, △t is the time sampling interval, represents the local event slope at the center of the 3D sampling data volume, M represents the number of time samples in the 3D sampling data volume, i represents the index of the time sampling point, N represents the total number of channels in the 3D sampling data volume, j represents the index of the channel number, x and y represent the spatial coordinates of the seismic channel relative to the center of the 3D sampling data volume, Represents the average of all trace data along the spatial direction of the three-dimensional sampling data volume, a(j) represents the adaptive coefficient based on the variational principle, Since this application uses the coherence calculation formula of the adaptive standard trace to calculate the coherence value of each three-dimensional sampling data volume, the influence of the lateral variation of the seismic data amplitude on the coherence value calculation is reduced. This method not only maintains the high efficiency of the C2 algorithm, but also has the resolution of the C3 algorithm, significantly improving the ability to resolve discontinuities such as faults and unconformities. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0036] Figure 1 A schematic flow chart of a method for calculating a seismic coherence volume provided in one embodiment of the present application;

[0037] Figure 2 A schematic diagram of various calculation windows provided in one embodiment of the present application;

[0038] Figure 3 3D post-stack seismic data volume A is calculated using the C2 algorithm, the C3 algorithm, and the embodiment of the present application for a horizontal slice comparison diagram at a certain time position;

[0039] Figure 4 A comparison diagram of vertical slices of a coherence volume at a certain spatial position calculated using the C2 algorithm, the C3 algorithm, and the embodiment of the present application for a 3D post-stack seismic data volume A;

[0040] Figure 5 3D post-stack seismic data volume B is calculated using the C2 algorithm, the C3 algorithm, and the embodiment of the present application for a horizontal slice comparison diagram at a certain time position;

[0041] Figure 6 A comparison diagram of vertical slices of a coherence volume at a certain spatial position calculated using the C2 algorithm, the C3 algorithm, and the embodiment of the present application for a 3D post-stack seismic data volume B;

[0042] The above drawings illustrate specific embodiments of the present application, which will be described in more detail below. These drawings and the textual description are not intended to limit the scope of the present application in any way, but rather to illustrate the concepts of the present application to those skilled in the art by reference to specific embodiments. DETAILED DESCRIPTION

[0043] Exemplary embodiments will be described in detail herein, with examples illustrated in the accompanying drawings. In the following description, when referring to the drawings, identical numerals in different figures represent identical or similar elements, unless otherwise indicated. The embodiments described in the following exemplary embodiments are not intended to represent all embodiments consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with certain aspects of the present application, as detailed in the appended claims.

[0044] Among existing coherence volume calculation methods, the C2 algorithm can produce erroneous results when used with seismic data that exhibits dramatic lateral amplitude variations. These erroneous results can reduce resolution and even lead to misjudgments in seismic interpretation. While the C3 algorithm offers higher resolution, its computational efficiency is low. Typically, the C3 algorithm takes ten times longer to calculate the same data volume than the C2 algorithm, a drawback that hinders its practical application.

[0045] The present application provides a method, apparatus, device, and computer-readable storage medium for calculating a seismic coherence volume. Based on the variational principle, an adaptive standard trace is introduced to calculate the least squares error between the seismic data and the adaptive standard trace within a calculation window. This method can eliminate the influence of lateral variations in the amplitude of the seismic data on the coherence volume calculation, thereby improving the C2 algorithm. While maintaining the computational efficiency of the C2 algorithm, this new method significantly enhances the ability to resolve discontinuities such as faults and unconformities.

[0046] The following specific embodiments describe in detail the technical solution of the present application and how the technical solution of the present application solves the above-mentioned technical problems. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments. The embodiments of the present application will be described below in conjunction with the accompanying drawings.

[0047] Figure 1 The flowchart of the method for calculating the seismic coherence volume provided in one embodiment of the present application is as follows. Figure 1 As shown, the earthquake coherence volume calculation method provided in this embodiment may include the following steps:

[0048] S110: According to a preset calculation window, with a single data sampling point as the center, three-dimensional sampling data volumes are sequentially extracted from the three-dimensional post-stack seismic data volume of the area to be processed.

[0049] In the embodiment of the present application, the area to be processed is a target area for which a corresponding coherence volume is to be obtained. In the target area, data acquisition is performed using a seismic data acquisition device, and then the acquired seismic data is subjected to three-dimensional stacking and migration processing to obtain a three-dimensional post-stack seismic data volume.

[0050] In order to obtain the seismic coherence volume results based on the above-mentioned three-dimensional post-stack seismic data volume, it is necessary to divide it into a series of three-dimensional sampling data volumes based on a preset calculation window with a single data sampling point as the center, and calculate the coherence value at the center position of each three-dimensional sampling data volume. Finally, the seismic coherence volume results corresponding to the three-dimensional post-stack seismic data volume of the area to be processed are obtained.

[0051] In the embodiment of the present application, the pre-set calculation window is: a sampling range consisting of a preset number of seismic traces and a preset length of time sampling points around a single data sampling point.

[0052] For example, the number of seismic traces in the calculation window of a single data sampling point can be as follows: Figure 2 As shown. Among them, Figure 2 -(a), Figure 2 -(b), Figure 2 -(c) are calculation windows containing 3, 5 and 9 seismic traces respectively. The solid curves in the figure represent the selected seismic traces, and the dashed-solid lines represent the unselected seismic traces.

[0053] Specifically, Figure 2 -The calculation window shown in (a) consists of the seismic trace where the sampling point is located and the two surrounding seismic traces; Figure 2 -The calculation window shown in (b) consists of the seismic trace where the sampling point is located and the four seismic traces before, after, and on the left and right sides; Figure 2 -The calculation window shown in (c) consists of the seismic trace where the sampling point is located and the eight surrounding seismic traces.

[0054] It can be understood that when selecting a three-dimensional sampling data volume from a three-dimensional post-stack seismic data volume, only three-dimensional sampling data volumes that can satisfy the requirement of having sufficient data sampling points within the calculation window are selected. The three-dimensional sampling data volume includes a center point located at the center of the three-dimensional sampling data volume and all data sampling points within the calculation window centered on the center point. Some data sampling points in the three-dimensional post-stack seismic data volume that are not suitable as the center points of the three-dimensional sampling data volume will not be used as center points. All data sampling points in the three-dimensional post-stack seismic data volume that can be used as center points will be used as center points and the corresponding three-dimensional sampling data volume will be extracted.

[0055] In the time direction, the time sampling window in the calculation window is preferably in the range of 20ms to 100ms. Of course, the time sampling window can be adjusted according to actual needs and is not specifically limited here. The number of sampling points in the time sampling window is guaranteed to be an odd number.

[0056] It should be noted that the three-dimensional sampling data volume extracted with the data sampling point as the center in the three-dimensional post-stack seismic data volume is equivalent to a rectangular data space.

[0057] S120: Determine the local event slope at the center position of each three-dimensional sampling data volume using a direct search method.

[0058] In the embodiment of the present application, the process of determining the local event slope at the center position of each three-dimensional sampling data volume using the direct search method may specifically include steps S121 to S123.

[0059] S121: Determine a plurality of slope sampling points according to a preset value range of the local event slope in the first direction and the second direction and a sampling interval.

[0060] For example, suppose the slope of the local event axis in the first direction (X direction) and the second direction (Y direction) are both in the range of [-p max ,p max ], the sampling interval is Δp. For example, p max It is 0.1ms / m and Δp is 0.01ms / m.

[0061] S122: At the center position of each three-dimensional sampling data volume, a coherence calculation formula based on multi-channel similarity is used to calculate the coherence value corresponding to each slope sampling point.

[0062] In this embodiment, for each slope sampling point (p x ,p y ), the coherence calculation formula based on multi-channel similarity (i.e., C2 algorithm) used to calculate its corresponding coherence value is:

[0063]

[0064] Wherein, C2 represents the coherence value calculated by the C2 algorithm, u represents the three-dimensional sampling data volume obtained in step S110, t represents the sampling time at the center position of the three-dimensional sampling data volume, Δt is the time sampling interval, (p x ,p y ) represents the slope value along the X and Y directions, M represents the number of time samples in the three-dimensional sampling data volume, i represents the index of the time sampling point, N represents the total number of channels in the three-dimensional sampling data volume, j represents the index of the channel number, and x and y represent the spatial coordinates of the seismic channel relative to the center of the three-dimensional sampling data volume. The function a(x) is transformed into a variational problem. a (x,a(x)=∫ t (u(t,x)-u a (t,x)) 2 )dt.

[0065] According to the variational principle, Q a Satisfies the Euler-Lagrange equation:

[0066]

[0067] By solving the above equation, we can obtain the function a(x). The solution to this variational problem means that among all possible forms of a(x), the form of a(x) that satisfies the Euler-Lagrange equation will make the function Q reach a global minimum.

[0068] The above formula can be used to obtain the expression of function a(x):

[0069]

[0070] Based on the function a(x) obtained by solving the variational problem, this embodiment performs a discrete representation on the function a(x) to obtain the adaptive coefficient a(j), which is specifically:

[0071]

[0072] At this time, the adaptive coefficient of discrete expression is brought into the least squares error Q for discrete expression as follows:

[0073] Normalize the least squares error of the discrete representation. After normalization, the coherence calculation formula based on the variational principle can be obtained. The specific calculation formula is as follows:

[0074]

[0075] In this embodiment, after the local event axis slope is obtained through steps S121-S123, Then, it is brought into the normalized coherence calculation formula to obtain the coherence value based on the local event slope. The specific calculation formula is also the coherence calculation formula:

[0076]

[0077] Where: C v represents the coherence value based on the variational principle, u represents the three-dimensional sampling data volume taken out in step S110, t represents the sampling time at the center position of the three-dimensional sampling data volume, Δt is the time sampling interval, represents the local event slope at the center of the 3D sampling data volume obtained in step S120, M represents the number of time samples in the 3D sampling data volume, i represents the index of the time sampling point, N represents the total number of traces in the 3D sampling data volume, j represents the index of the trace number, x and y represent the spatial coordinates of the seismic trace relative to the center of the 3D sampling data volume, Represents the average of all trace data along the spatial direction of the three-dimensional sampling data volume, a(j) represents the adaptive coefficient based on the variational principle,

[0078] To this end, in a specific embodiment of the present application, the above-mentioned coherence calculation formula based on the variational principle can be used to calculate the least squares error between the adaptive standard trace and each seismic trace in the calculation window according to the local event axis slope at the center position of each three-dimensional sampling data volume, and determine the coherence value at the center position of each three-dimensional sampling data volume, thereby obtaining the seismic coherence volume result of the area to be processed.

[0079] The seismic coherence volume calculation method provided in the embodiment of the present application is based on the variational principle and improves the calculation formula of the C2 algorithm by introducing an adaptive standard trace. This method can not only handle the problem of lateral variation in the amplitude of seismic data, but also maintains the high efficiency of the C2 algorithm while also having the resolution of the C3 algorithm, significantly improving the ability to resolve discontinuities such as faults and unconformities.

[0080] The following compares the seismic coherence volume calculation results obtained for different three-dimensional post-stack seismic data volumes using the C2 algorithm, the C3 algorithm, and the coherence volume calculation method provided in the embodiments of the present application.

[0081] like Figure 3 As shown, Figure 3 -(a) is a horizontal slice of the 3D post-stack seismic data volume A at a certain time position; Figure 3 -(b) Calculation results using the C2 algorithm; Figure 3 -(c) Calculation results using the C3 algorithm; Figure 3 -(d) Calculation results using the coherence calculation formula based on the variational principle.

[0082] like Figure 4 As shown, Figure 4 -(a) A vertical slice of 3D post-stack seismic data volume A at a certain spatial position; Figure 4 -(b) Calculation results using the C2 algorithm; Figure 4 -(c) Calculation results using the C3 algorithm; Figure 4 -(d) Calculation results using the coherence calculation formula based on the variational principle.

[0083] like Figure 5 As shown, Figure 5 -(a) A horizontal slice of 3D post-stack seismic data volume B at a certain time position; Figure 5 -(b) Calculation results using the C2 algorithm; Figure 5 -(c) Calculation results using the C3 algorithm; Figure 5 -(d) Calculation results using the coherence calculation formula based on the variational principle.

[0084] like Figure 6 As shown, Figure 6 -(a) A vertical slice of the 3D post-stack seismic data volume B at a certain spatial position; Figure 6 -(b) Calculation results using the C2 algorithm; Figure 6 -(c) Calculation results using the C3 algorithm; Figure 6 -(d) Calculation results using the coherence calculation formula based on the variational principle.

[0085] FIG7 is a schematic diagram of the structure of a seismic coherence volume calculation device provided in one embodiment of the present application. As shown in FIG7 , the seismic coherence volume calculation device 70 provided in this embodiment includes:

[0086] The data volume selection module 701 is used to extract three-dimensional sampling data volumes from the three-dimensional post-stack seismic data volume of the area to be processed in sequence based on a preset calculation window and with a single data sampling point as the center;

[0087] The slope calculation module 702 is used to determine the local event slope at the center position of each three-dimensional sampling data volume using a direct search method;

[0088] a coherence value calculation module 703, configured to determine the coherence value at the center of each three-dimensional sampling data volume based on the local event slope at the center of each three-dimensional sampling data volume and a coherence calculation formula based on the variational principle, thereby obtaining a seismic coherence volume result of the area to be processed;

[0089] The coherence calculation formula based on the variational principle is:

[0090]

[0091] Where: C v represents the coherence value based on the variational principle, u represents the three-dimensional sampling data volume, t represents the sampling time at the center position of the three-dimensional sampling data volume, △t is the time sampling interval, represents the local event slope at the center of the 3D sampling data volume, M represents the number of time samples in the 3D sampling data volume, i represents the index of the time sampling point, N represents the total number of channels in the 3D sampling data volume, j represents the index of the channel number, x and y represent the spatial coordinates of the seismic channel relative to the center of the 3D sampling data volume, Represents the average of all trace data along the spatial direction of the three-dimensional sampling data volume, a(j) represents the adaptive coefficient based on the variational principle,

[0092] In a possible implementation, the slope calculation module 702 may include a slope sampling point determination unit, a slope sampling point coherence value calculation unit, and a slope determination unit.

[0093] A slope sampling point determination unit, configured to determine a plurality of slope sampling points of the local event slope according to a preset value range of the local event slope in the first direction and the second direction and a sampling interval;

[0094] A slope sampling point coherence value calculation unit is used to calculate the coherence value corresponding to each slope sampling point at the center position of each three-dimensional sampling data volume using a coherence calculation formula based on multi-channel similarity;

[0095] The slope determination unit is used to select the local event slope of the slope sampling point corresponding to the minimum coherence value at the center position of each three-dimensional sampling data volume as the local event slope at the center position of each three-dimensional sampling data volume.

[0096] In a possible implementation, the calculation window is a sampling range consisting of a preset number of seismic traces and a preset length of time sampling points around a single data sampling point.

[0097] In a possible implementation, the number of seismic traces in the calculation window is 3, 5, or 9.

[0098] The seismic coherence volume calculation device provided in this embodiment can execute the method provided in the above method embodiment. Its implementation principle and technical effects are similar and will not be described in detail in this embodiment.

[0099] Figure 8 is a schematic diagram of the structure of a seismic coherence volume calculation device provided in one embodiment of the present application. As shown in Figure 8, the seismic coherence volume calculation device 80 provided in this embodiment includes: at least one processor 801 and a memory 802. Optionally, the device 80 also includes a communication component 803. The processor 801, the memory 802, and the communication component 803 are connected via a bus 804.

[0100] During the specific implementation process, at least one processor 801 executes the computer-executable instructions stored in the memory 802, so that the at least one processor 801 performs the above method.

[0101] The specific implementation process of the processor 801 can be found in the above method embodiment. Its implementation principle and technical effects are similar and will not be repeated here in this embodiment.

[0102] In the above embodiments, it should be understood that the processor may be a central processing unit (CPU), other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASICs), etc. A general-purpose processor may be a microprocessor or any conventional processor. The steps of the method disclosed in the present invention may be directly executed by a hardware processor or by a combination of hardware and software modules in the processor.

[0103] The memory may include a high-speed memory (Random Access Memory, RAM), and may also include a non-volatile memory (NVM), such as at least one disk memory.

[0104] The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus. Buses can be classified into address buses, data buses, and control buses. For ease of illustration, the buses in the drawings of this application are not limited to just one bus or just one type of bus.

[0105] An embodiment of the present application also provides a computer program product, including a computer program, which implements the above-mentioned seismic coherence volume calculation method when executed by a processor.

[0106] The present application also provides a computer-readable storage medium, in which computer-executable instructions are stored. When a processor executes the computer-executable instructions, the above-mentioned seismic coherence volume calculation method is implemented.

[0107] The above-mentioned readable storage medium can be implemented by any type of volatile or non-volatile memory device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic disk or optical disk. The readable storage medium can be any available medium that can be accessed by a general-purpose or special-purpose computer.

[0108] An exemplary readable storage medium is coupled to a processor so that the processor can read information from the readable storage medium and write information to the readable storage medium. Of course, the readable storage medium can also be an integral part of the processor. The processor and the readable storage medium can be located in an application specific integrated circuit (ASIC). Of course, the processor and the readable storage medium can also exist in the device as discrete components.

[0109] The division of units is merely a logical functional division; actual implementations may employ alternative divisions, such as combining or integrating multiple units or components into another system, or omitting or disabling certain features. Furthermore, any direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection between devices or units, either through an interface, electrical, mechanical, or other means.

[0110] Units described as separate components may or may not be physically separate, and components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of these units may be selected to achieve the purpose of this embodiment according to actual needs.

[0111] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0112] If the function is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, and other media that can store program code.

[0113] Those skilled in the art will appreciate that all or part of the steps in the above-described method embodiments can be implemented using hardware associated with program instructions. The aforementioned program can be stored in a computer-readable storage medium. When executed, the program performs the steps of the above-described method embodiments. The aforementioned storage medium includes various media capable of storing program code, such as ROM, RAM, magnetic disks, or optical disks.

[0114] Finally, it should be noted that those skilled in the art will readily identify other embodiments of the present invention after considering the specification and practicing the invention disclosed herein. The present invention is intended to cover any variations, uses, or adaptations of the present invention that follow the general principles of the present invention and include common knowledge or customary techniques in the art not disclosed herein. The present invention is not limited to the precise structure described above and illustrated in the accompanying drawings, and various modifications and variations may be made without departing from the scope thereof. The scope of the present invention is limited solely by the appended claims.

Claims

1. A method for calculating a seismic coherence volume, characterized in that: include: According to the preset calculation window, with a single data sampling point as the center, three-dimensional sampling data volumes are sequentially extracted from the three-dimensional post-stack seismic data volume of the area to be processed; Determining the local event slope at the center position of each three-dimensional sampling data volume using a direct search method; Determining the coherence value at the center of each three-dimensional sampling data volume according to the local event slope at the center of each three-dimensional sampling data volume and a coherence calculation formula based on the variational principle, to obtain a seismic coherence volume result of the area to be processed; The coherence calculation formula based on the variational principle is: Where: C v represents the coherence value based on the variational principle, u represents the three-dimensional sampling data volume, t represents the sampling time at the center position of the three-dimensional sampling data volume, △t is the time sampling interval, represents the local event slope at the center of the 3D sampling data volume, M represents the number of time samples in the 3D sampling data volume, i represents the index of the time sampling point, N represents the total number of channels in the 3D sampling data volume, j represents the index of the channel number, x and y represent the spatial coordinates of the seismic channel relative to the center of the 3D sampling data volume, Represents the average of all trace data along the spatial direction of the three-dimensional sampling data volume, a(j) represents the adaptive coefficient based on the variational principle, 2. The method for calculating the seismic coherence volume according to claim 1, wherein: The process of determining the local event slope at the center position of each three-dimensional sampling data volume by using a direct search method includes: Determining a plurality of slope sampling points according to a preset value range of the local event slope in the first direction and the second direction and a sampling interval; At the center position of each of the three-dimensional sampling data volumes, a coherence value corresponding to each slope sampling point is calculated using a coherence calculation formula based on multi-channel similarity; The local event slope of the slope sampling point corresponding to the minimum coherence value at the center position of each three-dimensional sampling data volume is selected as the local event slope at the center position of each three-dimensional sampling data volume.

3. The method for calculating seismic coherence volume according to claim 1, wherein: The calculation window is a sampling range consisting of a preset number of seismic traces and a preset length of time sampling points around a single data sampling point.

4. The method for calculating the seismic coherence volume according to claim 3, wherein: The number of seismic traces in the calculation window is 3, 5 or 9.

5. A seismic coherence volume calculation device, characterized in that: include: The data volume selection module is used to extract three-dimensional sampling data volumes from the three-dimensional post-stack seismic data volume of the area to be processed in sequence based on a preset calculation window and with a single data sampling point as the center; a slope calculation module, configured to determine the local event slope at the center position of each of the three-dimensional sampling data volumes by using a direct search method; a coherence value calculation module, configured to determine the coherence value at the center of each three-dimensional sampling data volume according to the local event slope at the center of each three-dimensional sampling data volume and a coherence calculation formula based on the variational principle, thereby obtaining a seismic coherence volume result of the area to be processed; The coherence calculation formula based on the variational principle is: Where: C v represents the coherence value based on the variational principle, u represents the three-dimensional sampling data volume, t represents the sampling time at the center position of the three-dimensional sampling data volume, △t is the time sampling interval, represents the local event slope at the center of the 3D sampling data volume, M represents the number of time samples in the 3D sampling data volume, i represents the index of the time sampling point, N represents the total number of channels in the 3D sampling data volume, j represents the index of the channel number, x and y represent the spatial coordinates of the seismic channel relative to the center of the 3D sampling data volume, Represents the average of all trace data along the spatial direction of the three-dimensional sampling data volume, a(j) represents the adaptive coefficient based on the variational principle, 6. The earthquake coherence volume calculation device according to claim 5, characterized in that: The slope calculation module includes: A slope sampling point determination unit, configured to determine a plurality of slope sampling points according to a preset value range of the local event slope in the first direction and the second direction and a sampling interval; a slope sampling point coherence value calculation unit, configured to calculate the coherence value corresponding to each slope sampling point at the center position of each three-dimensional sampling data volume using a coherence calculation formula based on multi-channel similarity; The slope determining unit is configured to select the local event slope of the slope sampling point corresponding to the minimum coherence value at the center position of each of the three-dimensional sampling data volumes as the local event slope at the center position of each of the three-dimensional sampling data volumes.

7. An electronic device, characterized in that: include: Memory, processor; The memory stores computer-executable instructions; The processor executes the computer-executable instructions stored in the memory, so that the processor executes the seismic coherence volume calculation method according to any one of claims 1 to 4.

8. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer-executable instructions, which, when executed by a processor, are used to implement the seismic coherence volume calculation method according to any one of claims 1 to 4.

Images