Regular interference wave elimination method and device based on independent interference source and medium

Abstract

The invention belongs to the technical field of geophysical prospecting engineering, and discloses a regular interference wave elimination method based on an independent interference source, a device thereof and a medium, and the method comprises the steps: combining an SPS format file of an observation system of the independent interference source with linear dynamic correction, two-dimensional Fourier forward transformation, two-dimensional Fourier inverse transformation and reaction dynamic correction, and obtaining an SPS format file of the observation system of the independent interference source; and regular interference waves in the related seismic shot gather data are eliminated. Apparatus and media based on the method are capable of storing and executing computer programs based on the method. According to the method, regular interference waves can be quickly, accurately and efficiently positioned and eliminated, and the working efficiency is high. The method is suitable for regular interference wave elimination in seismic data.

Description

Technical Field

[0001] The present invention belongs to the technical field of geophysical prospecting engineering, and relates to a method for eliminating regular interference waves, specifically a method for eliminating regular interference waves based on independent interference sources, as well as its equipment and medium. Background Art

[0002] The regular interference waves generated by mechanical vibration usually appear as hyperbolic event axes with different starting times on seismic records. This type of interference wave has strong energy, and it takes a long time and is difficult to completely eliminate using traditional elimination methods. Since this type of interference wave has the same hyperbolic characteristics as the effective wave, it is very difficult to completely eliminate it or it will damage the effective wave during elimination.

[0003] Currently, for this type of regular interference wave generated by mechanical vibration, the main elimination method is to directly delimit a time window at the position where the regular interference wave exists, transform the data within the time window to the F-K domain, and circle the position of the regular interference wave by observing the position difference between the effective wave and the regular interference wave for elimination. Since the regular interference wave has the same hyperbolic event axis shape as the effective wave, when transformed to the F-K domain, the effective wave and the regular interference wave cannot be completely separated, and the effective wave will be damaged while eliminating the regular interference wave, resulting in a poor elimination effect. Moreover, the delimited time window range cannot be too large and can only be removed in segments, with very low elimination efficiency.

[0004] How to quickly, accurately, and efficiently locate and eliminate regular interference waves and improve work efficiency is the key to improving exploration quality and ability. Summary of the Invention

[0005] To solve the above deficiencies in the prior art, the present invention aims to provide a method for eliminating regular interference waves based on independent interference sources, so as to achieve the purpose of quickly, accurately, and efficiently locating and eliminating regular interference waves.

[0006] A method for eliminating regular interference waves based on independent interference sources includes the following steps carried out in sequence: S1. Compile an SPS format file of the observation system related to the independent interference source, update it to the SPS format file of the original observation system to obtain a new observation system SPS format file, and locate the shot gather data and offset data related to the independent interference source in the new observation system SPS format file; Among them, the SPS format file refers to the shell processing support format file, which is the SPS format of seismic exploration auxiliary data.

[0007] S2. Select a single shot data from the shot gather data, limit the offset range, and perform linear dynamic correction to obtain a single shot data with flattened event axes.

[0008] S3. Take the single-shot data after in-phase axis leveling, perform two-dimensional forward Fourier transform, observe the differences between the regular interference waves and the effective waves in the F-K domain, eliminate the regions where the regular interference waves are located, and obtain the data in the F-K domain with regular interference waves eliminated.

[0009] S4. Perform two-dimensional inverse Fourier transform on the data with regular interference waves eliminated, and perform normal moveout correction to obtain the single-shot data in the time domain with regular interference waves eliminated.

[0010] S5. Eliminate the regular interference waves from the remaining single-shot data in the shot gather data according to steps S2 - S4.

[0011] The independent interference source is an interference source that excites regular interference waves and the complete regular interference waves excited by it can be found in a certain single-shot data.

[0012] Among them, the original observation system refers to the seismic data obtained by superimposing the data of the regular interference waves excited by the independent interference source and the seismic waves excited by the conventional seismic source.

[0013] As a limitation of the present invention, the compilation is to regard the independent interference source as an independent excitation seismic source, and according to the position coordinates, elevation of the independent interference source, and the information of the number of geophones and offset involved, compile the SPS format file of the observation system of the independent interference source.

[0014] As a further limitation of the present invention, the positioning is to use the position and receiving array information of the independent interference source to locate the shot gather data and offset data related to the independent interference source.

[0015] As a further further limitation of the present invention, the linear normal moveout correction is to perform linear normal moveout correction according to the apparent velocity differences between the regular interference waves and the effective waves, given the apparent velocity of the regular interference waves, and level the in-phase axis of the regular interference waves; the leveling is to level the in-phase axis of each regular interference wave hyperbola into a straight line, and the slope of the straight line is 0.

[0016] As another limitation of the present invention, the two-dimensional Fourier transform is the F-K transform.

[0017] Another object of the present invention is to provide a computer device and a medium to implement a computer program for running or storing and executing the above method for eliminating regular interference waves based on an independent interference source.

[0018] The present invention also provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, it implements the method for eliminating regular interference waves based on an independent interference source described in any one of the above technical solutions.

[0019] The present invention also provides a computer-readable storage medium storing a computer program for executing the method for eliminating regular interference waves based on independent interference sources according to any one of the above technical solutions.

[0020] Due to the adoption of the above technical solutions, the beneficial effects achieved by the present invention compared with the prior art are as follows: (1) The present invention updates the SPS format file of the observation system of the independent interference source to the SPS format file of the original observation system, and quickly and accurately locates the shot gather data involving the independent interference waves according to the information in the SPS format file of the independent interference source. Compared with manually finding out the independent interference waves one by one, the present invention greatly improves the work efficiency and accuracy.

[0021] (2) By limiting the offset range involved in the single-shot data, the present invention performs linear dynamic correction on the regular interference waves, flattens the in-phase axis of the regular interference waves, and obtains regular interference waves that are straight lines in the time domain. After two-dimensional Fourier transform, they become a point in the F-K domain, enabling the efficient and accurate elimination of regular interference waves.

[0022] (3) The present invention combines the independent interference source with linear dynamic correction, two-dimensional Fourier forward transform to the F-K domain to eliminate regular interference waves, and two-dimensional Fourier inverse transform to the time domain for reverse dynamic correction, which can quickly, accurately and efficiently locate and eliminate regular interference waves, while improving the signal-to-noise ratio of seismic data and work efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] The present invention will be further described in detail below with reference to the drawings and specific embodiments.

[0024] Figure 1 It is the single-shot data containing regular interference waves in Embodiment 1 of the present invention; Figure 2 It is the single-shot data after linear dynamic correction of the regular interference waves in Embodiment 1 of the present invention; Figure 3 It is the display of the selected time window position and the data within the time window in the single-shot data of Embodiment 1 of the present invention transformed into the F-K domain, Figure 3 In (a), the selected time window is located at 0 ms - 2600 ms, the time window width is 30 channels, and the display of the data within the time window transformed into the F-K domain, Figure 3 In (b), the selected time window is located at 3100 ms - 4900 ms, the time window width is 30 channels, and the display of the data within the time window transformed into the F-K domain; Figure 4 It is the single-shot data of the data after eliminating regular interference waves in the F-K domain of Embodiment 1 of the present invention transformed back to the time domain by two-dimensional Fourier inverse transform; Figure 5This is the single-shot data in the time domain after eliminating the regular interference wave and performing reverse correction in Embodiment 1 of the present invention. Detailed implementation manners

[0025] The present invention will be further described in detail below through specific embodiments and accompanying drawings. It should be understood that the described embodiments are only used to explain the present invention and do not limit the present invention.

[0026] Embodiment 1 A method for eliminating regular interference waves based on independent interference sources This embodiment is a method for eliminating regular interference waves based on independent interference sources. The specific method includes the following steps carried out in sequence: S0. Data acquisition work During the actual acquisition process in the work area, the field operators record the positions of the independent interference sources. After the data is returned, the indoor processing personnel, according to the actual positions of the independent interference sources, find a single-shot data with complete interference waves at the corresponding positions in the shot gather data, as shown in Figure 1 .

[0027] From Figure 1 , it can be known that the conventional seismic source is located at x = 250 m, the independent interference source is located at x = 60 m, the trace interval is 1 m, the regular interference wave excited by the independent interference source involves 30 traces. Thus, the coordinate position of the independent interference source is obtained as (60, 0), the number of traces is 30, the elevation is 0, the number of geophone points is 30, and the maximum offset is 15 m.

[0028] S1. Regard the independent interference source as an independent excitation seismic source. According to the position coordinates, elevation, the number of geophone points involved, offset and other information of the independent interference source, compile the observation system SPS format file of the independent interference source, update it into the original observation system SPS format file, and use the coordinate position and receiving array information of the independent interference source to locate the shot gather data and offset data related to the independent interference source.

[0029] S2. Limit Figure 1 the offset range involved in the single-shot data to 0 - 15. According to the apparent velocity difference between the regular interference wave and the effective wave, the apparent velocity of the regular interference wave is given as 800 m / s, and linear dynamic correction is performed to obtain the regular interference wave with the in-phase axis of each regular interference wave being a straight line and the slope being 0, as shown in Figure 2 .

[0030] S3. Open two time windows on the regular interference wave, perform two-dimensional Fourier forward transform on the data within the time windows to the F - K domain, and use the position difference between the regular interference wave and the effective wave in the F - K domain to delineate the range of the regular interference wave for elimination, as shown in Figure 3 . From Figure 3 , it can be known that Figure 3In (a), the selected time window is located at 0 ms - 2600 ms, the time window width is 30 channels, and the display of the data within the time window transformed into the F-K domain. Figure 3 In (b), the selected time window is located at 3100 ms - 4900 ms, the time window width is 30 channels, and the display of the data within the time window transformed into the F-K domain. Each hyperbolic phase axis of the regular interference wave within the time window is a straight line with a slope of 0. Therefore, the regular interference wave within the time window becomes a point in the F-K domain.

[0031] S4. The data in the F-K domain after eliminating the regular interference wave is inverse two-dimensional Fourier transformed into the time domain as shown in Figure 4 , and moveout correction is performed to obtain the single-shot data with the regular interference wave eliminated as shown in Figure 5 .

[0032] S5. According to steps S2 - S4, the regular interference wave is eliminated from the remaining single-shot data in the shot gather data.

[0033] Through the method of this embodiment, the rapid and efficient elimination of the regular interference wave in seismic data is achieved. Through the SPS format file information of the observation system of independent interference sources, the relevant shot gather data can be obtained quickly and efficiently. Compared with manual search, the working efficiency of the present invention is significantly improved. By giving the apparent velocity of the regular interference wave, linear moveout correction is performed on the regular interference wave to flatten the phase axis of the regular interference wave, and then it is two-dimensional Fourier transformed forward into the F-K domain to eliminate the regular interference wave. Subsequently, it is inverse two-dimensional Fourier transformed into the time domain for moveout correction to obtain the seismic data with the regular interference wave eliminated. This method can quickly, accurately, and efficiently achieve the elimination of the regular interference wave, while improving the signal-to-noise ratio and working efficiency.

[0034] Embodiment 2 A computer device This embodiment provides a computer device, which includes: a memory, a processor, and a computer program stored on the memory and executable on the processor to implement the method for eliminating regular interference waves based on independent interference sources in Embodiment 1.

[0035] The memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, and the computer program products may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory, etc. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc.

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

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

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

[0039] Embodiment 3 A computer-readable storage medium This embodiment provides a computer-readable storage medium storing a computer program, and when the computer program is executed by a processor, it implements the rule interference wave cancellation method based on independent interference sources in Embodiment 1.

[0040] The computer-readable storage medium stores non-transitory computer-readable instructions. When the non-transitory computer-readable instructions are run by a processor, all or part of the steps of the methods in the foregoing embodiments are executed.

[0041] The above 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 removable hard disk), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).

[0042] It should be noted that the above are only preferred embodiments of the present invention and are not used to limit the present invention. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art can still modify the technical solutions described in the above embodiments or perform equivalent replacements for some of the technical features. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention shall be included in the scope of protection of the claims of the present invention.

Claims

1. A method for eliminating regular interference waves based on independent interference sources, characterized in that, it includes the following steps carried out in sequence: S1. Compile an SPS format file of the observation system of the relevant independent interference source, update it to the SPS format file of the original observation system to obtain a new observation system SPS format file, and locate the shot gather data and offset data related to the independent interference source in the new observation system SPS format file; S2. Select a single-shot data from the shot gather data, limit the offset range, and perform linear dynamic correction to obtain a single-shot data with flattened event axes; S3. Take the single-shot data with flattened event axes, perform two-dimensional forward Fourier transform, observe the difference between regular interference waves and effective waves in the F-K domain, eliminate the area where regular interference waves are located, and obtain data with regular interference waves eliminated in the F-K domain; S4. Perform two-dimensional inverse Fourier transform on the data with regular interference waves eliminated, and perform reverse dynamic correction to obtain a single-shot data with regular interference waves eliminated in the time domain; S5. Eliminate regular interference waves for the remaining single-shot data in the shot gather data according to steps S2 - S4; The independent interference source is an interference source that excites regular interference waves and the complete regular interference waves it excites can be found in a certain single-shot data.

2. The method for eliminating regular interference waves based on independent interference sources according to claim 1, characterized in that, The compilation is to regard the independent interference source as an independent excitation source, and compile the SPS file of the observation system of the independent interference source according to the position coordinates, elevation of the independent interference source, and the number of geophone points and offset information it involves.

3. The method for eliminating regular interference waves based on independent interference sources according to claim 2, characterized in that, The positioning is to locate the shot gather data and offset data related to the independent interference source by using the position and receiving array information of the independent interference source.

4. The method for eliminating regular interference waves based on independent interference sources according to claim 3, characterized in that, The linear dynamic correction is to perform linear dynamic correction according to the apparent velocity difference between regular interference waves and effective waves, given the apparent velocity of regular interference waves, and flatten the event axes of regular interference waves; The flattening is to flatten the event axes of each regular interference wave hyperbola into a straight line, and the slope of the straight line is 0.

5. The method for eliminating regular interference waves based on independent interference sources according to any one of claims 1 - 4, characterized in that, The two-dimensional Fourier transform is the F-K transform.

6. A computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method for eliminating regular interference waves based on independent interference sources according to any one of claims 1 to 5.

7. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program for executing the method for eliminating regular interference waves based on independent interference sources according to any one of claims 1 to 5.

Images