Processing method and device for quickly smoothing seismic data amplitude, electronic equipment and storage medium

Abstract

The invention discloses a processing method and device for quickly smoothing seismic data amplitude, electronic equipment and a storage medium. The method comprises the following steps: inputting seismic data to be processed; calculating the amplitude coefficient of each point in each trace in the seismic data along the time direction; carrying out fast smooth calculation on the amplitude coefficient along the space direction; and calculating amplitude-smoothed seismic data based on the fast smoothing calculation result and the seismic data. According to the method, rapid smooth calculation is carried out on the amplitude coefficient in the space direction, data of all windows of the first window are read and sorted, positions of numerical values outside the windows of each subsequent window are found out through the dichotomy and then taken out, positions of newly added numerical values are found out through the dichotomy and then added, and finally, the amplitude coefficient is obtained. Repeated operations such as data reading and sorting on each window are avoided, the calculated amount of sorting is reduced, and the calculation efficiency is improved.

Description

Technical Field

[0001] This invention belongs to the field of geophysical exploration technology, and more specifically, relates to a method, electronic device, storage medium and apparatus for rapidly smoothing seismic data amplitude. Background Technology

[0002] With the deepening development of oil and gas exploration, exploration and development are gradually shifting towards complex oil and gas reservoirs. High-resolution seismic exploration and time-shift seismic techniques are becoming increasingly important, thus raising the requirements for seismic data quality. Currently, there is a growing demand for amplitude consistency processing in seismic data. Amplitude inconsistencies can lead to artifacts in subsequent migration imaging. This technique is being widely applied to seismic data amplitude processing, which is of great significance for improving the quality of seismic data.

[0003] However, existing smoothing methods require finding surrounding points for each point, sorting these values, excluding the maximum and minimum values, and then calculating the average of the remaining data. This results in repeated calculations such as sorting and averaging for each point, leading to reduced computational efficiency and making it difficult to meet the timeliness requirements of production.

[0004] The information disclosed in the background section of this invention is intended only to enhance the understanding of the general background of this invention and should not be construed as an admission or in any way implying that such information constitutes prior art known to those skilled in the art. Summary of the Invention

[0005] The purpose of this invention is to propose a method, electronic device, storage medium, and apparatus for processing seismic data amplitude quickly and smoothly, thereby solving the problem of large computational load and long computation time caused by repetitive operations such as sorting and averaging required for each point in traditional methods, and improving computational efficiency.

[0006] To achieve the above objectives, the present invention proposes a method, electronic device, storage medium, and apparatus for rapidly smoothing seismic data amplitude.

[0007] According to a first aspect of the present invention, a method for rapidly smoothing seismic data amplitude processing is proposed, comprising:

[0008] Calculate the amplitude coefficient of each point in each trace of the seismic data along the time direction;

[0009] Calculate the amplitude coefficient of each sampling point in the earthquake data;

[0010] The amplitude coefficient is calculated quickly and smoothly along the spatial direction;

[0011] Based on the results of the fast smoothing calculation and the earthquake data, the amplitude-smoothed earthquake data is calculated.

[0012] Optionally, the calculation expression of the amplitude coefficient is as follows:

[0013]

[0014] where d(t j , x i ) is the seismic data at the j-th time point of the i-th seismic trace, x i is the i-th seismic trace, 0 ≤ i < n, n is the total number of seismic traces, t j is the j-th time point, 0 ≤ j < T, T is the total number of time points of the seismic data, and agc_t is the window radius.

[0015] Optionally, the fast smoothing calculation of the amplitude coefficient along the spatial direction includes:

[0016] For the first sampling point, arrange the amplitude coefficients amp(t j , [0, 1,..., r]) in ascending order and store them in an array. After excluding the maximum and minimum values of the amplitude coefficients in the array, calculate the average value of the remaining amplitude coefficients, where r is the window radius;

[0017] For subsequent sampling points, during the window movement, when there is an amplitude coefficient in the array that goes out of the window range, determine the position of the amplitude coefficient in the array by the binary search method and then delete the amplitude coefficient from the array; when there is a new amplitude coefficient entering the window range, use the binary search method to determine the position of the new amplitude coefficient in the array, insert the new amplitude coefficient into the array at the position to obtain a new array, and calculate the average value of the remaining amplitude coefficients after excluding the maximum and minimum values of the amplitude coefficients in the array.

[0018] Optionally, the expression for calculating the average value of the remaining amplitude coefficients is as follows:

[0019]

[0020] where r1 is the number of the maximum and minimum values of the amplitude coefficients removed from the array, N is the number of amplitude coefficients in the array, that is, the number of sampling points, and tmp is the array.

[0021] Optionally, the calculation expression of the seismic data after amplitude smoothing is as follows:

[0022]

[0023] According to the second aspect of the present invention, a processing device for fast smoothing the amplitude of seismic data is proposed, including:

[0024] An input module for inputting the seismic data to be processed;

[0025] A first calculation module for calculating the amplitude coefficient of each point of each trace in the seismic data along the time direction;

[0026] A smoothing processing module for performing fast smoothing calculation on the amplitude coefficient along the spatial direction;

[0027] A second calculation module for calculating the seismic data with smoothed amplitude based on the fast smoothing calculation result and the seismic data.

[0028] Optionally, the calculation expression of the amplitude coefficient is:

[0029]

[0030] where d(t j , x i ) is the seismic data of the j-th time point of the i-th seismic trace, x i is the i-th seismic trace, 0 ≤ i < n, n is the total number of seismic traces, t j is the j-th time point, 0 ≤ j < T, T is the total number of time points of the seismic data, and agc_t is the window radius.

[0031] Optionally, the performing fast smoothing calculation on the amplitude coefficient along the spatial direction includes:

[0032] For the first sampling point, arrange the amplitude coefficients amp(t j , [0, 1,..., r]) in ascending order and store them in an array. After excluding the maximum and minimum amplitude coefficients in the array, calculate the average value of the remaining amplitude coefficients, where r is the window radius;

[0033] For subsequent sampling points, during the window movement, when there is an amplitude coefficient in the array that goes out of the window range, determine the position of the amplitude coefficient in the array by the bisection method and then delete the amplitude coefficient from the array; when there is a new amplitude coefficient entering the window range, determine the position of the new amplitude coefficient in the array by the bisection method, insert the new amplitude coefficient into the array at the position to obtain a new array, and after excluding the maximum and minimum amplitude coefficients in the array, calculate the average value of the remaining amplitude coefficients.

[0034] According to the third aspect of the present invention, an electronic device is provided, and the electronic device includes:

[0035] At least one processor; and,

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

[0037] The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the fast smoothing seismic data amplitude processing method described in any of the first aspects.

[0038] According to a fourth aspect of the invention, a non-transitory computer-readable storage medium is provided, which stores computer instructions for causing a computer to perform the fast smoothing seismic data amplitude processing method described in any of the first aspects.

[0039] The beneficial effects of this invention are as follows: This invention calculates the amplitude coefficient of each sampling point of the seismic data, and then performs a fast smoothing calculation on the amplitude coefficient along the spatial direction. Based on the fast smoothing calculation result and the seismic data, the smoothed seismic data is calculated. This solves the problem of large computational load and long calculation time in the smoothing calculation process of conventional algorithms, such as frequent sorting and averaging. By using a recursive approach, the position for removing or inserting data is found by binary search each time the window slides, thereby deleting or adding data, reducing the computational load of sorting, and improving computational efficiency.

[0040] The system of the present invention has other features and advantages that will be apparent from or will be set forth in detail in the accompanying drawings and following detailed description, which together serve to explain the particular principles of the invention. Attached Figure Description

[0041] The above and other objects, features and advantages of the present invention will become more apparent from the accompanying drawings, in which like reference numerals generally denote like parts.

[0042] Figure 1 A flowchart illustrating the steps of a method for rapidly smoothing seismic data amplitude according to the present invention is shown.

[0043] Figure 2 A schematic diagram of the seismic data to be processed according to Embodiment 2 of the present invention is shown.

[0044] Figure 3 A schematic diagram of seismic data after rapid smoothing processing according to Embodiment 2 of the present invention is shown. Detailed Implementation

[0045] 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 so that the present invention will be more thorough and complete, and will fully convey the scope of the present invention to those skilled in the art.

[0046] As Figure 1 shown, a method for processing seismic data with fast and smooth amplitude according to the present invention includes:

[0047] Inputting the seismic data to be processed;

[0048] Calculating the amplitude coefficient of each point of each trace in the seismic data along the time direction;

[0049] Performing fast and smooth calculation on the amplitude coefficient along the spatial direction;

[0050] Calculating the seismic data with smoothed amplitude based on the fast and smooth calculation result and the seismic data.

[0051] Specifically, the present invention first sets the automatic gain compensation parameter: window radius agc_t, sets the fast and smooth parameter: smooth radius r, and the number r1 of the maximum and minimum values to be removed; then inputs the seismic data d(t,x) to be processed, and there are n traces in total for the seismic data; for trace x i (0 <= i < n), calculating the amplitude coefficient of each sampling point of this trace:

[0052]

[0053] Performing fast and smooth calculation on the obtained amp along the spatial direction. At time slice t j below, 0 <= j < T, specifically as follows:

[0054] For the first sampling point, store the values of amp(t j ,[0,1,...,r]) in the array tmp. The current number of elements in the tmp array is N, and calculate

[0055] For subsequent sampling points, first, during the window movement, if there is data outside the window range, find the position of this value through the binary search method and delete it from the tmp array. If there is new data entering the window, use the binary search method to find the position and insert it into the tmp array, and calculate to obtain

[0056] For the obtained Calculate That is, the result after smoothing the seismic data.

[0057] In this invention, all data of all windows are read from the first window and sorted. For each subsequent window, the values outside the window are found and retrieved through the binary search method, and the newly added values are also found and added through the binary search method, avoiding repetitive operations such as data reading and sorting for each window, solving the problems of large computational amount and long calculation time in steps such as frequent sorting and calculating the average value in the conventional algorithm during the smoothing calculation process. Through a recursive method, the position for removing or inserting data is found through the binary search method each time the sliding window moves, thereby deleting or adding data, reducing the computational amount of sorting, and improving the calculation efficiency.

[0058] In one example, the calculation expression of the amplitude coefficient is:

[0059]

[0060] where d(t j , x i ) is the seismic data at the j-th time point of the i-th seismic trace, x i is the i-th seismic trace, 0 ≤ i < n, n is the total number of seismic data traces, t j is the j-th time point, 0 ≤ j < T, T is the total number of time points of the seismic data, and agc_t is the window radius.

[0061] In one example, the fast smoothing calculation of the amplitude coefficient along the spatial direction includes:

[0062] For the first sampling point, the amplitude coefficients amp(t j , [0, 1,..., r]) are sorted in ascending order and stored in an array. After excluding the maximum and minimum values of the amplitude coefficients in the array, the average value of the remaining amplitude coefficients is calculated, where r is the window radius;

[0063] For subsequent sampling points, during the process of window movement, when there are amplitude coefficients in the array that are out of the window range, the position of the amplitude coefficient in the array is determined through the binary search method and then the amplitude coefficient is deleted from the array; when there are new amplitude coefficients entering the window range, the position of the new amplitude coefficient in the array is determined through the binary search method, and the new amplitude coefficient is inserted into the array according to the position to obtain a new array. After excluding the maximum and minimum values of the amplitude coefficients in the array, the average value of the remaining amplitude coefficients is calculated.

[0064] In one example, the expression for calculating the average value of the remaining amplitude coefficients is:

[0065]

[0066] where r1 is the number of the maximum and minimum values of the amplitude coefficients removed from the array, N is the number of amplitude coefficients in the array, that is, the number of sampling points, and tmp is the array.

[0067] In one example, the calculation expression of the amplitude-smoothed seismic data is:

[0068]

[0069] The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments, but it is not a limitation of the present invention. It should be noted that, without conflict, the embodiments and features in the embodiments of the present invention can be combined with each other.

[0070] Embodiment 1

[0071] This embodiment provides a method for processing seismic data to quickly smooth the amplitude, including:

[0072] Input the seismic data to be processed;

[0073] Calculate the amplitude coefficient of each point in each trace of the seismic data along the time direction;

[0074] Perform a fast smoothing calculation on the amplitude coefficient along the spatial direction;

[0075] Calculate the amplitude-smoothed seismic data based on the fast smoothing calculation result and the seismic data.

[0076] The calculation expression of the amplitude coefficient is:

[0077]

[0078] where d(t j , x i ) is the seismic data of the i-th seismic trace at the j-th time point, x i is the i-th seismic trace, 0 ≤ i < n, n is the total number of seismic traces, t j is the j-th time point, 0 ≤ j < T, T is the total number of time points of the seismic data, and agc_t is the window radius.

[0079] Performing a fast smoothing calculation on the amplitude coefficient along the spatial direction includes:

[0080] For the first sampling point, arrange the amplitude coefficients amp(t j , [0, 1,..., r]) of all sampling points within the window in ascending order and store them in an array, and calculate the average value of the remaining amplitude coefficients after excluding the maximum and minimum values of the amplitude coefficients in the array, where r is the window radius;

[0081] For subsequent sampling points, during the window movement, when there are amplitude coefficients in the array that fall outside the window range, the position of the amplitude coefficient in the array is determined by the binary search method and then the amplitude coefficient is deleted from the array; when there are new amplitude coefficients entering the window range, the position of the new amplitude coefficient in the array is determined by the binary search method, and the new amplitude coefficient is inserted into the array according to the position to obtain a new array. After excluding the maximum and minimum values of the amplitude coefficients in the array, the average value of the remaining amplitude coefficients is calculated.

[0082] The expression for calculating the average value of the remaining amplitude coefficients is:

[0083]

[0084] where r1 is the number of the maximum and minimum values of the amplitude coefficients removed from the array, N is the number of amplitude coefficients in the array, that is, the number of sampling points, and tmp is the array.

[0085] The calculation expression for the seismic data after amplitude smoothing is:

[0086]

[0087] Example 2

[0088] This embodiment provides a method for processing seismic data to quickly smooth the amplitude, including:

[0089] Set the automatic gain compensation parameter: window radius agc_t,

[0090] Set the fast smoothing parameters: smoothing radius r, the number of the maximum and minimum values removed r1;

[0091] Input the seismic data d(t,x) to be processed, and there are n traces in total for the seismic data; for trace i (0 <= i < n), calculate the amplitude coefficient of each sampling point of this trace:

[0092]

[0093] Perform fast smoothing calculation on the obtained amp in the spatial direction. At time slice t j under the condition that 0 ≤ j < T, the specific steps are as follows:

[0094] For the first sampling point, store the values of amp(t j ,[0,1,...,r]] in the array tmp. The current number of elements in the tmp array is N, and calculate

[0095] For subsequent sampling points, first, during the window movement, if there is data outside the window range, use the binary search method to find the position of the value and delete it from the tmp array. If new data enters the window, use the binary search method to find the position and insert it into the tmp array, and then calculate to obtain

[0096] For the obtained Calculate That is, the result of the smoothed seismic data.

[0097] Figure 2 Shows the seismic data to be processed, Figure 3 Shows the use of the processing method of this embodiment for Figure 2 The processing result of the seismic data therein.

[0098] Embodiment III

[0099] This embodiment provides a processing device for quickly smoothing the amplitude of seismic data, including:

[0100] An input module for inputting the seismic data to be processed;

[0101] A first calculation module for calculating the amplitude coefficient of each point in each trace of the seismic data along the time direction;

[0102] A smoothing processing module for quickly smoothing and calculating the amplitude coefficient along the spatial direction;

[0103] A second calculation module for calculating the seismic data with smoothed amplitude based on the quick smoothing calculation result and the seismic data.

[0104] The calculation expression of the amplitude coefficient is:

[0105]

[0106] Where, d(t j ,x i ) is the seismic data of the i-th seismic trace at the j-th time point, x i is the i-th seismic trace, 0 ≤ i < n, n is the total number of seismic traces, t j is the j-th time point, 0 ≤ j < T, T is the total number of time points of the seismic data, and agc_t is the window radius.

[0107] Quickly smoothing and calculating the amplitude coefficient along the spatial direction includes:

[0108] For the first sampling point, arrange the amplitude coefficients amp(t j ,[0,1,...,r]) in ascending order and store them in an array, and calculate the average value of the remaining amplitude coefficients after excluding the maximum and minimum values of the amplitude coefficients in the array, where r is the window radius;

[0109] For subsequent sampling points, during the window movement, if an amplitude coefficient in the array leaves the window range, the amplitude coefficient is removed from the array after determining its position in the array using a binary search method. When a new amplitude coefficient enters the window range, the position of the new amplitude coefficient in the array is determined using a binary search method, and the new amplitude coefficient is inserted into the array according to its position to obtain a new array. After excluding the maximum and minimum values ​​of the amplitude coefficients in the array, the average value of the remaining amplitude coefficients is calculated.

[0110] Example 4

[0111] This disclosure also provides an electronic device, which includes:

[0112] At least one processor; and,

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

[0114] The memory stores instructions that can be executed by the at least one processor, which enables the at least one processor to perform the fast smoothing seismic data amplitude processing method in Embodiment 1.

[0115] An electronic device according to embodiments 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, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0116] The processor may be a central processing unit (CPU) or other form of processing unit with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of this disclosure, the processor is used to execute computer-readable instructions stored in the memory.

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

[0118] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0119] Example 5

[0120] This disclosure provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the fast smoothing seismic data amplitude processing method in Embodiment 1.

[0121] A computer-readable storage medium according to embodiments of the present disclosure stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the methods described in the foregoing embodiments of the present disclosure are performed.

[0122] The aforementioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or portable hard drive), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).

[0123] The various embodiments of the present invention have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for rapidly smoothing the amplitude of seismic data, characterized in that, include: Input the seismic data to be processed; Calculate the amplitude coefficient of each point in each trace of the seismic data along the time direction; The amplitude coefficient is calculated quickly and smoothly along the spatial direction; Based on the results of the fast smoothing calculation and the earthquake data, the amplitude-smoothed earthquake data is calculated.

2. The method for processing seismic data amplitude for rapid smoothing according to claim 1, characterized in that, The expression for calculating the amplitude coefficient is as follows: where, d(t j , x i ) is the seismic data at the j-th time point of the i-th seismic trace at (t j , x i ). x i is the i-th seismic trace, 0 ≤ i < n, where n is the total number of seismic traces, and t j is the j-th time slice, 0 ≤ t j < T, where T is the total time of the seismic data, and agc_t is the window radius.

3. The method for processing seismic data amplitude for rapid smoothing according to claim 2, characterized in that, The rapid and smooth calculation of the amplitude coefficient along the spatial direction includes: For the first sampling point, the amplitude coefficient amp(t) j The amplitude coefficients ([0,1,...,r]) are arranged in order of size and stored in an array. The average value of the remaining amplitude coefficients is calculated after excluding the maximum and minimum values ​​of the amplitude coefficients in the array, where r is the window radius. For subsequent sampling points, during the window movement, if an amplitude coefficient in the array leaves the window range, the amplitude coefficient is removed from the array after determining its position in the array using a bisection method. When a new amplitude coefficient enters the window range, the position of the new amplitude coefficient in the array is determined using a bisection method, and the new amplitude coefficient is inserted into the array at the determined position to obtain a new array. After excluding the maximum and minimum values ​​of the amplitude coefficients in the array, the average value of the remaining amplitude coefficients is calculated.

4. The method for processing seismic data amplitude for rapid smoothing according to claim 3, characterized in that, The expression for calculating the average value of the remaining amplitude coefficients is as follows: Where r1 is the number of maximum and minimum amplitude coefficients removed from the array, N is the number of amplitude coefficients in the array, i.e. the number of sampling points, and tmp is the array.

5. The method for processing seismic data amplitude for rapid smoothing according to claim 1, characterized in that, The calculation expression for the amplitude-smoothed seismic data is as follows:

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 that can be executed by the at least one processor to enable the at least one processor to perform the fast smoothing seismic data amplitude processing method according to any one of claims 1-5.

7. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium stores computer instructions for causing a computer to perform the fast smoothing seismic data amplitude processing method as described in any one of claims 1-5.

8. A processing device for rapidly smoothing seismic data amplitude, characterized in that, include: The input module is used to input the seismic data to be processed; The first calculation module is used to calculate the amplitude coefficient of each point in each trace of the seismic data along the time direction; A smoothing module is used to smooth the amplitude coefficient along the spatial direction; The second calculation module is used to calculate amplitude-smoothed seismic data based on the fast smoothing calculation results and the seismic data.

9. The processing apparatus for rapidly smoothing seismic data amplitude according to claim 8, characterized in that, The expression for calculating the amplitude coefficient is as follows: where, d(t j , x i ) is the seismic data at the j-th time point of the i-th seismic trace, x i is the i-th seismic trace, 0 ≤ i < n, n is the total number of seismic traces, t j is the j-th time point, 0 ≤ j < T, T is the total number of time points of the seismic data, and agc_t is the window radius.

10. The processing apparatus for rapidly smoothing seismic data amplitude according to claim 9, characterized in that, The smoothing process of the amplitude coefficient along the spatial direction includes: For the first sampling point, the amplitude coefficient amp(t) j The amplitude coefficients ([0,1,...,r]) are arranged in order of size and stored in an array. The average value of the remaining amplitude coefficients is calculated after excluding the maximum and minimum values ​​of the amplitude coefficients in the array, where r is the window radius. For subsequent sampling points, during the window movement, if an amplitude coefficient in the array leaves the window range, the amplitude coefficient is removed from the array after determining its position in the array using a bisection method. When a new amplitude coefficient enters the window range, the position of the new amplitude coefficient in the array is determined using a bisection method, and the new amplitude coefficient is inserted into the array at the determined position to obtain a new array. After excluding the maximum and minimum values ​​of the amplitude coefficients in the array, the average value of the remaining amplitude coefficients is calculated.

Images