Seismic data processing method and device for eliminating window effect, electronic equipment and storage medium

Abstract

The invention discloses a seismic data processing method for eliminating a window effect, electronic equipment, a storage medium and a device. The method comprises the following steps: inputting to-be-processed seismic data; performing window division on the to-be-processed seismic data based on the first time window length, the first space window length and the overlapping length to obtain a plurality of first windows and corresponding first window seismic data; adjusting each first window based on the second time window length and the second space window length to obtain a plurality of second windows and corresponding second window seismic data; performing interpolation or denoising processing on the second window seismic data to obtain third window seismic data; extracting the third window seismic data based on the first window to obtain fourth window seismic data; and combining all the fourth window seismic data, and carrying out linear weighted transition on the overlapped region to obtain the seismic data without the boundary effect. The method can effectively eliminate the boundary effect of the window, and ensures the stability and rationality of the calculation result.

Description

Technical Field

[0001] This invention belongs to the field of geophysical exploration technology, and more specifically, relates to a seismic data processing method, electronic device, storage medium, and apparatus for eliminating window effects. 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-shifted seismic techniques are becoming increasingly important, thus raising the requirements for seismic data quality. Currently, seismic data processing methods are becoming increasingly complex, and the computational load is also gradually increasing. Therefore, a windowed data processing workflow has emerged. This technology is increasingly being widely applied to seismic data noise suppression and data interpolation, which is of great significance for improving the quality of seismic data.

[0003] However, the handling of overlapping areas in different window data merging stages can affect the final result. Conventional methods use ramps to transition between different windows, but this can sometimes still produce anomalies at the boundaries, leading to poor denoising performance.

[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 seismic data processing method, electronic device, storage medium, and apparatus for eliminating window effects. This invention addresses the problem that when processing seismic data in multiple windows, significant boundary effects occur at the window boundaries during subsequent data merging, leading to abrupt changes in signal energy. By eliminating these boundary effects in seismic data processing, the invention ensures the stability and rationality of the calculation results, thereby obtaining better interpolated or denoised data.

[0006] To achieve the above objectives, the present invention proposes a seismic data processing method, electronic device, storage medium, and apparatus for eliminating window effects.

[0007] According to a first aspect of the present invention, a seismic data processing method for eliminating window effects is proposed, comprising:

[0008] Input the seismic data to be processed;

[0009] The seismic data to be processed is divided into windows based on the first time window length, the first spatial window length, and the overlap length, resulting in multiple first windows and corresponding first window seismic data.

[0010] Adjust each of the first windows based on the second time window length and the second spatial window length to obtain multiple second windows and corresponding second window seismic data;

[0011] The seismic data in the second window are interpolated or denoised to obtain the seismic data in the third window.

[0012] Based on the first window, the seismic data in the third window is extracted to obtain the seismic data in the fourth window;

[0013] Based on the first window, all the seismic data from the fourth window are merged, and the overlapping areas are linearly weighted to obtain seismic data that eliminates boundary effects.

[0014] Optionally, the length of the first time window is shorter than the length of the second time window.

[0015] Optionally, the length of the first spatial window is less than the length of the second spatial window.

[0016] Optionally, the center point of the second window is the same as that of the first window.

[0017] Optionally, the linear weighting of the linear weighting over-weighting is:

[0018] In the left window, w1 = (Oz) / O;

[0019] In the right window, w2 = z / O;

[0020] Where w1 is the linear weighted weight of the left window, w2 is the linear weighted weight of the right window, O is the window overlap length, and z is the coordinate of the overlapping region.

[0021] According to a second aspect of the present invention, a seismic data processing apparatus for eliminating window effects is provided, comprising:

[0022] The input module is used to input the seismic data to be processed;

[0023] The partitioning module is used to partition the seismic data to be processed into windows based on the first time window length, the first spatial window length, and the overlap length, so as to obtain multiple first windows and corresponding first window seismic data.

[0024] An adjustment module is used to adjust each of the first windows based on the second time window length and the second spatial window length to obtain multiple second windows and corresponding second window seismic data.

[0025] The processing module is used to interpolate or denoise the seismic data in the second window to obtain the seismic data in the third window.

[0026] The extraction module is used to extract the seismic data from the third window based on the first window to obtain the seismic data from the fourth window;

[0027] The merging module is used to merge all the fourth window seismic data based on the first window, and to perform linear weighting transition on the overlapping areas to obtain seismic data that eliminates boundary effects.

[0028] Optionally, the length of the first time window is shorter than the length of the second time window.

[0029] Optionally, the length of the first spatial window is less than the length of the second spatial window.

[0030] According to a third aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0031] At least one processor; and,

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

[0033] The memory stores instructions executable by the at least one processor, which, when executed, enable the at least one processor to perform the seismic data processing method for eliminating window effects as described in any of the first aspects.

[0034] 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 seismic data processing method for eliminating window effects as described in any of the first aspects.

[0035] The beneficial effects of this invention are as follows: By expanding the window based on the window division, the data range of the window is further expanded, and more data is introduced to participate in the calculation, so that the energy at the boundary is close to that of the adjacent window. In the subsequent data merging process, the basic range of the window is still used for the slope transition, thereby effectively eliminating the boundary effect of the window, ensuring the stability and rationality of the calculation results, and thus obtaining better interpolation or denoising data, improving the stability of seismic data processing results.

[0036] 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

[0037] 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.

[0038] Figure 1 A flowchart illustrating the steps of a seismic data processing method for eliminating window effects according to the present invention is shown.

[0039] Figure 2 A schematic diagram illustrating window division using a conventional method according to Embodiment 2 of the present invention is shown.

[0040] Figure 3 A schematic diagram of an enlarged window according to Embodiment 2 of the present invention is shown.

[0041] Figure 4 A schematic diagram of a noise model for seismic data to be processed according to Embodiment 2 of the present invention is shown.

[0042] Figure 5 This illustrates a conventional method for processing according to Embodiment 2 of the present invention. Figure 4 A schematic diagram of the results of the seismic data to be processed.

[0043] Figure 6 This illustrates a seismic data processing method based on Embodiment 2 of the present invention, which eliminates the window effect. Figure 4 A schematic diagram of the results of the seismic data to be processed. Detailed Implementation

[0044] The invention will now be described in more detail with reference to the accompanying drawings. While preferred embodiments of the invention are shown in the drawings, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that the invention will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.

[0045] like Figure 1 As shown, a seismic data processing method for eliminating window effects according to the present invention includes:

[0046] Input the seismic data to be processed;

[0047] Based on the first time window length, the first spatial window length, and the overlap length, the seismic data to be processed is divided into windows to obtain multiple first windows and corresponding first window seismic data.

[0048] Adjust each first window based on the second time window length and the second spatial window length to obtain multiple second windows and corresponding second window seismic data;

[0049] The seismic data in the second window is interpolated or denoised to obtain the seismic data in the third window.

[0050] Based on the first window, the seismic data in the third window is extracted to obtain the seismic data in the fourth window;

[0051] The first window is used to merge all the fourth window seismic data, and the overlapping areas are linearly weighted to obtain seismic data that eliminates boundary effects.

[0052] Specifically, this invention first divides the input seismic data into windows based on a set window range. Then, based on the window division, the window is expanded to further increase the data range. The expanded window data is then interpolated or denoised to introduce more data into the calculation, making the energy at the boundary close to that of adjacent windows. Data from the base window is extracted, and the extracted window data is merged from different windows. During the data merging process, the base range of the window is still used for a slope transition, thereby effectively eliminating the boundary effect of the window and ensuring the stability and rationality of the calculation results. This results in better interpolated or denoised data and improves the stability of the seismic data processing results.

[0053] In one example, the length of the first time window is less than the length of the second time window.

[0054] In one example, the length of the first spatial window is less than the length of the second spatial window.

[0055] In one example, the second window has the same center point as the first window.

[0056] In one example, the linearly weighted overweighting is:

[0057] In the left window, w1 = (Oz) / O;

[0058] In the right window, w2 = z / O;

[0059] Where w1 is the linear weighted weight of the left window, w2 is the linear weighted weight of the right window, O is the window overlap length, and z is the coordinate of the overlapping region.

[0060] The present invention will be further described below with reference to the accompanying drawings and specific embodiments, but this is not intended to limit the invention. It should be noted that, unless otherwise specified, the embodiments and features described in the embodiments of the present invention can be combined with each other.

[0061] Example 1

[0062] This embodiment provides a seismic data processing method for eliminating window effects, including:

[0063] Input the seismic data to be processed;

[0064] Based on the first time window length, the first spatial window length, and the overlap length, the seismic data to be processed is divided into windows to obtain multiple first windows and corresponding first window seismic data.

[0065] Adjust each first window based on the second time window length and the second spatial window length to obtain multiple second windows and corresponding second window seismic data;

[0066] The seismic data in the second window is interpolated or denoised to obtain the seismic data in the third window.

[0067] Based on the first window, the seismic data in the third window is extracted to obtain the seismic data in the fourth window;

[0068] The first window is used to merge all the fourth window seismic data, and the overlapping areas are linearly weighted to obtain seismic data that eliminates boundary effects.

[0069] The length of the first time window is less than the length of the second time window.

[0070] The length of the first spatial window is less than the length of the second spatial window.

[0071] The second window has the same center point as the first window.

[0072] The linear weighting of overweighted linear weighting is:

[0073] In the left window, w1 = (Oz) / O;

[0074] In the right window, w2 = z / O;

[0075] Where w1 is the linear weighted weight of the left window, w2 is the linear weighted weight of the right window, O is the window overlap length, and z is the coordinate of the overlapping region.

[0076] Example 2

[0077] This embodiment provides a seismic data processing method for eliminating window effects, including:

[0078] This method is applicable to interpolation and denoising algorithms with arbitrary window divisions. Taking adaptive subtraction as an example, the process is as follows:

[0079] Two-dimensional window parameter settings: time window length is win_t, spatial window length is win_x, time overlap length is overlap_t, spatial overlap length is overlap_x, and the length of the expanded time window. a dd_ w in_t, the length of the expanded space window add_ w in_x.

[0080] Adaptive subtraction parameters: The operator is filt, and the stability factor is λ.

[0081] (1) Input the seismic data to be denoised d(t, x) and the model data m(t, x);

[0082] (2) Figure 2 As shown, the seismic data d(t, x) and m(t, x) are divided into windows based on time, spatial window length, and overlap length. There are N windows. For the seismic data in a certain window, we use d i (t, x) and m i (t, x) represents, where 0 ≤ i < N.

[0083] (3) Figure 3 As shown, for the i-th window, we expand the temporal-spatial window from the center of the window, forming a window of length add_ w in_t and add_ w The expanded window of in_x and the corresponding expanded window seismic data are used... and express;

[0084] (4) Using existing adaptive subtraction techniques, the operator filt is mainly calculated by solving the toeplitz matrix:

[0085]

[0086] The Levinson-Durbin algorithm can be used to solve this problem and calculate the removed noise.

[0087]

[0088] (5) Extraction The original win_t and win_x data ranges in the window are used to obtain noise-removed seismic data r. i (t,x);

[0089] (6) Merge the noise-removed seismic data from different windows and perform linear weighting transition on the overlapping areas;

[0090] The window overlap length is O. At coordinate z within the overlap range, the linear weights are: left window: w1 = (Oz) / O; right window: w2 = z / O.

[0091] (7) After merging all window data, the final result r(t,x) is obtained.

[0092] Figure 4 The noise model of the seismic data to be processed is shown. Figure 5This demonstrates processing using conventional methods. Figure 4 The results of the seismic data to be processed Figure 6 This embodiment illustrates the process of processing. Figure 4 The results of the seismic data to be processed show that the method described in this implementation can effectively eliminate the boundary effect of the window, ensure the stability and rationality of the calculation results, thereby obtaining better interpolated or denoised data and improving the stability of the seismic data processing results.

[0093] Example 3

[0094] This embodiment provides a seismic data processing apparatus for eliminating window effects, including:

[0095] The input module is used to input the seismic data to be processed;

[0096] The partitioning module is used to partition the seismic data to be processed into windows based on the first time window length, the first spatial window length, and the overlap length, so as to obtain multiple first windows and corresponding first window seismic data.

[0097] The adjustment module is used to adjust each first window based on the second time window length and the second spatial window length to obtain multiple second windows and corresponding second window seismic data.

[0098] The processing module is used to interpolate or denoise the seismic data in the second window to obtain the seismic data in the third window.

[0099] The extraction module is used to extract seismic data from the third window based on the first window, and obtain seismic data from the fourth window.

[0100] The merging module is used to merge all the fourth window seismic data based on the first window, and to perform linear weighting transition on the overlapping areas to obtain seismic data that eliminates boundary effects.

[0101] The length of the first time window is less than the length of the second time window.

[0102] The length of the first spatial window is less than the length of the second spatial window.

[0103] The second window has the same center point as the first window.

[0104] The linear weighting of overweighted linear weighting is:

[0105] In the left window, w1 = (Oz) / O;

[0106] In the right window, w2 = z / O;

[0107] Where w1 is the linear weighted weight of the left window, w2 is the linear weighted weight of the right window, O is the window overlap length, and z is the coordinate of the overlapping region.

[0108] Example 4

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

[0110] At least one processor; and,

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

[0112] The memory stores instructions that can be executed by the at least one processor, which enables the at least one processor to perform the seismic data processing method for eliminating window effects in Embodiment 1.

[0113] 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.

[0114] 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.

[0115] 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.

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

[0117] Example 5

[0118] This disclosure provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the seismic data processing method for eliminating window effects as described in Embodiment 1.

[0119] 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.

[0120] 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).

[0121] 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 seismic data processing method for eliminating window effects, characterized in that, include: Input the seismic data to be processed; The seismic data to be processed is divided into windows based on the first time window length, the first spatial window length, and the overlap length, resulting in multiple first windows and corresponding first window seismic data. Adjust each of the first windows based on the second time window length and the second spatial window length to obtain multiple second windows and corresponding second window seismic data; The seismic data in the second window are interpolated or denoised to obtain the seismic data in the third window. Based on the first window, the seismic data in the third window is extracted to obtain the seismic data in the fourth window; Based on the first window, all the seismic data from the fourth window are merged, and the overlapping areas are linearly weighted to obtain seismic data that eliminates boundary effects.

2. The seismic data processing method for eliminating window effects according to claim 1, characterized in that, The length of the first time window is less than the length of the second time window.

3. The seismic data processing method for eliminating window effects according to claim 1, characterized in that, The length of the first spatial window is less than the length of the second spatial window.

4. The seismic data processing method for eliminating window effects according to claim 1, characterized in that, The second window has the same center point as the first window.

5. The seismic data processing method for eliminating window effects according to claim 1, characterized in that, The linear weighting of the over-weighted linear weighting is: In the left window, w1 = (Oz) / O; In the right window, w2 = z / O; Where w1 is the linear weighted weight of the left window, w2 is the linear weighted weight of the right window, O is the window overlap length, and z is the coordinate of the overlapping region.

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 seismic data processing method for eliminating window effects as described in 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 seismic data processing method for eliminating window effects as described in any one of claims 1-5.

8. A seismic data processing device for eliminating window effects, characterized in that, include: The input module is used to input the seismic data to be processed; The partitioning module is used to partition the seismic data to be processed into windows based on the first time window length, the first spatial window length, and the overlap length, so as to obtain multiple first windows and corresponding first window seismic data. An adjustment module is used to adjust each of the first windows based on the second time window length and the second spatial window length to obtain multiple second windows and corresponding second window seismic data. The processing module is used to interpolate or denoise the seismic data in the second window to obtain the seismic data in the third window. The extraction module is used to extract the seismic data from the third window based on the first window to obtain the seismic data from the fourth window; The merging module is used to merge all the seismic data from the fourth window based on the first window, and to perform linear weighting transition on the overlapping areas to obtain seismic data that eliminates boundary effects.

9. The seismic data processing apparatus for eliminating window effects according to claim 8, characterized in that, The length of the first time window is less than the length of the second time window.

10. The seismic data processing apparatus for eliminating window effects according to claim 8, characterized in that, The length of the first spatial window is less than the length of the second spatial window.

Images