Method, device, medium and equipment for suppressing harmonic interference of vibroseis

Abstract

The application discloses a method, device, medium and equipment for suppressing harmonic interference of a controllable source, and the method comprises the following steps: obtaining the amplitudes of second-order harmonic, third-order harmonic and fourth-order harmonic signals based on original data of the controllable source; obtaining a preset signal waveform of each order of harmonic based on the scanning start and end frequencies of a fundamental wave and the amplitudes of each order of harmonic signal; obtaining order harmonic related data based on the preset signal waveform of each order of harmonic and the original data of the controllable source for the preset signal waveform of each order of harmonic; and subtracting the order harmonic related data from the original data of the controllable source to obtain controllable source data after the order harmonic is removed. The method for suppressing harmonic interference of the controllable source based on the correlation characteristics can better suppress the harmonic signals of the controllable source, the signal-to-noise ratio of the obtained controllable source data is obviously improved, the harmonic interference suppression effect is obvious, and the controllable source data imaging effect is effectively improved.

Description

Technical Field

[0001] This invention belongs to the field of oil and gas geophysical exploration technology, specifically relating to a method, device, medium, and electronic equipment for suppressing controllable source harmonic interference based on relevant characteristics. Background Technology

[0002] In the process of seismic exploration and production, seismic waves can usually be generated using a controlled source, and seismic data can be collected using instruments such as detectors.

[0003] Controlled-source seismic (CRS) technology is highly favored in the exploration field due to its advantages such as safety, environmental friendliness, low energy consumption, high energy output, and flexible operation. It is currently one of the most important seismic acquisition technologies in seismic exploration, playing an increasingly important role in seismic exploration of complex surfaces in western China. However, due to the unique excitation method of CRS, harmonic interference is inevitably generated during excitation. This is a significant noise difference from well-shot excitation. The presence of harmonics degrades the quality of CRS data, reduces data resolution, and poses significant challenges to subsequent processing, imaging, and interpretation of CRS data. Therefore, it is necessary to suppress harmonic interference to maximize the effectiveness of CRS. Consequently, harmonic interference suppression has become a key research focus in CRS technology. Currently, the main method for harmonic suppression in CRS is predictive deconvolution, which predicts harmonic signals and then suppresses them during processing. However, due to the difficulty of accurately predicting harmonic interference, current methods can only suppress a portion of the harmonic interference, not completely eliminate it.

[0004] The remaining harmonic interference is generally suppressed by the high coverage number of the observation system.

[0005] Therefore, a method with a significant harmonic interference suppression effect is needed. Summary of the Invention

[0006] The purpose of this invention is to propose a method that has a significant effect on suppressing harmonic interference.

[0007] In a first aspect, the present invention provides a method for suppressing harmonic interference from a controllable seismic source based on correlation characteristics, comprising: obtaining the amplitudes of second-order, third-order, and fourth-order harmonic signals respectively based on the original data of the controllable seismic source; obtaining a preset signal waveform for each harmonic based on the scanning start and end frequencies of the fundamental wave and the amplitude of each harmonic signal; for each preset signal waveform, obtaining harmonic correlation data based on the preset signal waveform and the original data of the controllable seismic source; and subtracting the original data of the controllable seismic source from the harmonic correlation data to obtain controllable seismic source data after removing the harmonics.

[0008] Optionally, the amplitudes of the second-order, third-order, and fourth-order harmonic signals are obtained through the following steps: extracting the force signal data from the original data of the controllable source; transforming the force signal data to the time-frequency domain; determining the amplitude of each harmonic signal in the time-frequency domain data; and obtaining the amplitudes of the second-order, third-order, and fourth-order harmonic signals according to the order of the harmonic signals.

[0009] Optionally, the preset signal waveform for each harmonic can be obtained through the following steps: design the scanning start and end frequencies of each harmonic signal based on the scanning start and end frequencies of the fundamental wave; for each harmonic, obtain the preset signal waveform of the harmonic based on the scanning start and end frequencies and amplitude values ​​of the harmonic.

[0010] Optionally, the order harmonic correlation data can be obtained through the following steps: performing correlation processing on the preset signal waveform of the order harmonic and the original data of the controllable source; converting the correlation-processed data to the time-frequency domain; and extracting the order harmonic correlation data from the converted data.

[0011] Secondly, the present invention also provides an electronic device, the electronic device comprising: a memory storing executable instructions; and a processor that executes the executable instructions in the memory to implement the above-described method for suppressing controllable source harmonic interference based on relevant characteristics.

[0012] Thirdly, the present invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for suppressing controllable source harmonic interference based on relevant characteristics.

[0013] Fourthly, the present invention also provides an apparatus for suppressing harmonic interference from a controllable seismic source based on relevant characteristics, comprising: an amplitude acquisition module, which acquires the amplitudes of second-order, third-order, and fourth-order harmonic signals based on the original data of the controllable seismic source; a harmonic preset signal establishment module, which acquires preset signal waveforms for each harmonic based on the scanning start and end frequencies of the fundamental wave and the amplitude of each harmonic signal; and a harmonic removal module, which, for each preset signal waveform, acquires harmonic correlation data based on the preset signal waveform and the original data of the controllable seismic source, and subtracts the original data of the controllable seismic source from the harmonic correlation data to obtain controllable seismic source data after removing the harmonics.

[0014] Optionally, the amplitudes of the second-order, third-order, and fourth-order harmonic signals are obtained through the following steps: extracting the force signal data from the original data of the controllable source; transforming the force signal data to the time-frequency domain; determining the amplitude of each harmonic signal in the time-frequency domain data; and obtaining the amplitudes of the second-order, third-order, and fourth-order harmonic signals according to the order of the harmonic signals.

[0015] Optionally, the preset signal waveform for each harmonic can be obtained through the following steps: design the scanning start and end frequencies of each harmonic signal based on the scanning start and end frequencies of the fundamental wave; for each harmonic, obtain the preset signal waveform of the harmonic based on the scanning start and end frequencies and amplitude values ​​of the harmonic.

[0016] Optionally, the order harmonic correlation data can be obtained through the following steps: performing correlation processing on the preset signal waveform of the order harmonic and the original data of the controllable source; converting the correlation-processed data to the time-frequency domain; and extracting the order harmonic correlation data from the converted data.

[0017] The beneficial effects of the present invention are as follows: the method of suppressing harmonic interference of controllable seismic source based on correlation characteristics can effectively suppress harmonic signals of controllable seismic source, significantly improve the signal-to-noise ratio of the obtained controllable seismic source data, and effectively improve the imaging effect of controllable seismic source data.

[0018] The present invention has other features and advantages, which 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

[0019] The above and other objects, features and advantages of the present invention will become more apparent from the more detailed description of exemplary embodiments of the invention in conjunction with the accompanying drawings, wherein the same reference numerals generally represent the same components in the exemplary embodiments of the invention.

[0020] Figure 1 A flowchart of a method for suppressing controllable source harmonic interference based on correlation characteristics according to an embodiment of the present invention is shown.

[0021] Figure 2 The diagram shows a harmonic suppression process before harmonic suppression, according to an embodiment of the present invention, of a method for suppressing harmonic interference from a controllable seismic source based on correlation characteristics.

[0022] Figure 3 The diagram shows a harmonic suppression result of a method for suppressing controllable source harmonic interference based on correlation characteristics according to an embodiment of the present invention.

[0023] Figure 4 A structural block diagram of a device for suppressing controllable source harmonic interference based on correlation characteristics according to an embodiment of the present invention is shown.

[0024] Explanation of reference numerals in the attached figures

[0025] 102. Amplitude acquisition module; 104. Harmonic preset signal establishment module; 106. Harmonic removal module. Detailed Implementation

[0026] Preferred embodiments of the invention will now be described in more detail. While preferred embodiments of the invention are described below, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0027] This invention provides a method for suppressing harmonic interference from a controllable seismic source based on correlation characteristics, comprising: obtaining the amplitudes of second-order, third-order, and fourth-order harmonic signals based on the original data of the controllable seismic source; obtaining a preset signal waveform for each harmonic based on the scanning start and end frequencies of the fundamental wave and the amplitude of each harmonic signal; obtaining harmonic correlation data for each harmonic preset signal waveform based on the harmonic preset signal waveform and the original data of the controllable seismic source; and subtracting the harmonic correlation data from the original data of the controllable seismic source to obtain controllable seismic source data after removing the harmonics.

[0028] Specifically, by using the regular characteristics between harmonic interference and fundamental wave signal, preset signal waveforms for each order of harmonics are designed. Correlation techniques are used to separate each order of harmonics and the fundamental wave in the time domain to obtain harmonic correlation data. Then, a separation method is used to separate the harmonic correlation data from the original data of the controllable source to obtain the original data of the controllable source without harmonic influence.

[0029] According to the exemplary implementation, the method for suppressing harmonic interference of controllable seismic sources based on relevant characteristics can effectively suppress harmonic signals of controllable seismic sources, significantly improve the signal-to-noise ratio of the obtained controllable seismic source data, and effectively improve the imaging effect of controllable seismic source data.

[0030] As an optional approach, the amplitudes of the second, third, and fourth harmonic signals are obtained through the following steps: extracting the force signal data from the original data of the controllable source; transforming the force signal data to the time-frequency domain; determining the amplitude of each harmonic signal in the time-frequency domain data; and obtaining the amplitudes of the second, third, and fourth harmonic signals respectively according to the order of the harmonic signal.

[0031] Specifically, the force signal of the raw data of the controllable seismic source is extracted, the force signal is transformed into the time-frequency domain, and the amplitude of each harmonic signal is read in the time-frequency domain data. Based on the order and amplitude of each harmonic, the amplitudes of the second-order, third-order, and fourth-order harmonic signals are obtained.

[0032] As an optional approach, the preset signal waveform for each harmonic can be obtained through the following steps: design the scanning start and end frequencies of each harmonic signal based on the scanning start and end frequencies of the fundamental wave; for each harmonic, obtain the preset signal waveform of that harmonic based on the scanning start and end frequencies and amplitude values ​​of that harmonic.

[0033] Specifically, based on the scanning start and end frequencies of the fundamental frequency, the scanning start and end frequencies of the second, third, and fourth harmonics are designed respectively. For example, if the fundamental frequency is 10-80Hz, then the second harmonic frequency is 20-160Hz, the third harmonic is 30-240Hz, and the fourth harmonic is 40-320Hz. Harmonics above the fourth order have lower energy and can be ignored; only the fourth harmonic needs to be calculated. Based on the start and end frequencies and amplitude values ​​of each harmonic, preset signal waveforms for each harmonic are designed, and SEGY data is output.

[0034] As an optional approach, harmonic correlation data can be obtained through the following steps: performing correlation processing on the preset harmonic signal waveform and the original data of the controllable source; converting the correlation-processed data to the time-frequency domain; and extracting harmonic correlation data from the converted data.

[0035] Specifically, the designed preset signal waveforms for each harmonic order are correlated with the original data from the controllable seismic source. This correlates the corresponding harmonic signals in the original data to time zero, while other harmonics, including the fundamental wave, are shifted to either side of time zero. For example, if a second-order harmonic signal is correlated with the original data, the second-order harmonic signal in the original data is concentrated at time zero. The fundamental wave, being lower in order than the second-order harmonic, has its correlated data on the positive time axis, while the third and fourth-order harmonics, being higher in order than the second-order harmonic, have their correlated data on the negative time axis. Switching the correlated data to the time-frequency domain reveals that the second-order harmonic is separated from other data in the time domain. The correlated second-order harmonic data in the time domain is then extracted to obtain the second-order harmonic correlation data. This extracted second-order harmonic correlation data is then convolved with the designed preset second-order harmonic signal waveform, transforming it back to the waveform of the second-order harmonic in the original data. This method allows for the sequential correlation of each harmonic signal with the original data, followed by the extraction of each harmonic noise.

[0036] Secondly, the present invention also provides an electronic device, comprising: a memory storing executable instructions; and a processor executing the executable instructions in the memory to implement the above-described method for suppressing controllable source harmonic interference based on relevant characteristics.

[0037] Thirdly, the present invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for suppressing controllable source harmonic interference based on relevant characteristics.

[0038] Fourthly, the present invention also provides an apparatus for suppressing harmonic interference from a controllable seismic source based on relevant characteristics, comprising: an amplitude acquisition module, which acquires the amplitudes of second-order, third-order, and fourth-order harmonic signals based on the original data of the controllable seismic source; a harmonic preset signal establishment module, which acquires the preset signal waveform of each harmonic based on the scanning start and end frequencies of the fundamental wave and the amplitude of each harmonic signal; and a harmonic removal module, which, for each harmonic preset signal waveform, acquires harmonic correlation data based on the harmonic preset signal waveform and the original data of the controllable seismic source, and subtracts the harmonic correlation data from the original data of the controllable seismic source to obtain the controllable seismic source data after removing the harmonics.

[0039] Specifically, by using the regular characteristics between harmonic interference and fundamental wave signal, preset signal waveforms for each order of harmonics are designed. Correlation techniques are used to separate each order of harmonics and the fundamental wave in the time domain to obtain harmonic correlation data. Then, a separation method is used to separate the harmonic correlation data from the original data of the controllable source to obtain the original data of the controllable source without harmonic influence.

[0040] According to the exemplary implementation, the method for suppressing harmonic interference of controllable seismic sources based on relevant characteristics can effectively suppress harmonic signals of controllable seismic sources, significantly improve the signal-to-noise ratio of the obtained controllable seismic source data, and effectively improve the imaging effect of controllable seismic source data.

[0041] As an optional approach, the amplitudes of the second, third, and fourth harmonic signals are obtained through the following steps: extracting the force signal data from the original data of the controllable source; transforming the force signal data to the time-frequency domain; determining the amplitude of each harmonic signal in the time-frequency domain data; and obtaining the amplitudes of the second, third, and fourth harmonic signals respectively according to the order of the harmonic signal.

[0042] Specifically, the force signal of the raw data of the controllable seismic source is extracted, the force signal is transformed into the time-frequency domain, and the amplitude of each harmonic signal is read in the time-frequency domain data. Based on the order and amplitude of each harmonic, the amplitudes of the second-order, third-order, and fourth-order harmonic signals are obtained.

[0043] As an optional approach, the preset signal waveform for each harmonic can be obtained through the following steps: design the scanning start and end frequencies of each harmonic signal based on the scanning start and end frequencies of the fundamental wave; for each harmonic, obtain the preset signal waveform of that harmonic based on the scanning start and end frequencies and amplitude values ​​of that harmonic.

[0044] Specifically, based on the scanning start and end frequencies of the fundamental frequency, the scanning start and end frequencies of the second, third, and fourth harmonics are designed respectively. For example, if the fundamental frequency is 10-80Hz, then the second harmonic frequency is 20-160Hz, the third harmonic is 30-240Hz, and the fourth harmonic is 40-320Hz. Harmonics above the fourth order have lower energy and can be ignored; only the fourth harmonic needs to be calculated. Based on the start and end frequencies and amplitude values ​​of each harmonic, preset signal waveforms for each harmonic are designed, and SEGY data is output.

[0045] As an optional approach, harmonic correlation data can be obtained through the following steps: performing correlation processing on the preset harmonic signal waveform and the original data of the controllable source; converting the correlation-processed data to the time-frequency domain; and extracting harmonic correlation data from the converted data.

[0046] Specifically, the designed preset signal waveforms for each harmonic order are correlated with the original data from the controllable seismic source. This correlates the corresponding harmonic signals in the original data to time zero, while other harmonics, including the fundamental wave, are shifted to either side of time zero. For example, if a second-order harmonic signal is correlated with the original data, the second-order harmonic signal in the original data is concentrated at time zero. The fundamental wave, being lower in order than the second-order harmonic, has its correlated data on the positive time axis, while the third and fourth-order harmonics, being higher in order than the second-order harmonic, have their correlated data on the negative time axis. Switching the correlated data to the time-frequency domain reveals that the second-order harmonic is separated from other data in the time domain. The correlated second-order harmonic data in the time domain is then extracted to obtain the second-order harmonic correlation data. This extracted second-order harmonic correlation data is then convolved with the designed preset second-order harmonic signal waveform, transforming it back to the waveform of the second-order harmonic in the original data. This method allows for the sequential correlation of each harmonic signal with the original data, followed by the extraction of each harmonic noise.

[0047] Example 1

[0048] Figure 1 A flowchart of a method for suppressing controllable source harmonic interference based on correlation characteristics according to an embodiment of the present invention is shown. Figure 2 The diagram shows a harmonic suppression process before harmonic suppression, according to an embodiment of the present invention, of a method for suppressing harmonic interference from a controllable seismic source based on correlation characteristics. Figure 3 The diagram shows a harmonic suppression result of a method for suppressing controllable source harmonic interference based on correlation characteristics according to an embodiment of the present invention.

[0049] Combination Figure 1 , Figure 2 and Figure 3 As shown, the method for suppressing harmonic interference from a controllable seismic source based on correlation characteristics includes:

[0050] Step 1: Based on the raw data from the controllable seismic source, obtain the amplitudes of the second-order, third-order, and fourth-order harmonic signals respectively;

[0051] Step 2: Based on the scanning start and end frequencies of the fundamental wave and the amplitude of each harmonic signal, obtain the preset signal waveform for each harmonic.

[0052] Step 3: For each harmonic preset signal waveform, based on the harmonic preset signal waveform and the original data of the controllable source, obtain the harmonic related data. Subtract the harmonic related data from the original data of the controllable source to obtain the controllable source data after removing the harmonics.

[0053] The amplitudes of the second, third, and fourth harmonic signals are obtained through the following steps: extracting the force signal data from the original data of the controllable source; transforming the force signal data to the time-frequency domain; determining the amplitude of each harmonic signal in the time-frequency domain data; and obtaining the amplitudes of the second, third, and fourth harmonic signals according to the order of the harmonic signal.

[0054] The preset signal waveform for each harmonic is obtained through the following steps: the scan start and end frequencies of each harmonic signal are designed based on the scan start and end frequencies of the fundamental wave; for each harmonic, the preset signal waveform for that harmonic is obtained based on the scan start and end frequencies and amplitude values ​​of that harmonic.

[0055] The harmonic correlation data is obtained through the following steps: the harmonic preset signal waveform is correlated with the original data of the controllable source; the correlated data is converted to the time-frequency domain; and the harmonic correlation data is extracted from the converted data.

[0056] Figure 2 and Figure 3 The image shows the profile effect before and after suppressing controllable source harmonics using the method of this application. It can be clearly seen that before harmonic suppression, there are a large number of harmonic effects in the profile, which seriously affect the data quality and signal-to-noise ratio. After harmonic suppression, the profile signal-to-noise ratio is significantly improved, and the harmonic interference suppression effect is good.

[0057] Example 2

[0058] Figure 4 A structural block diagram of a device for suppressing controllable source harmonic interference based on correlation characteristics according to an embodiment of the present invention is shown.

[0059] like Figure 4 As shown, the device for suppressing harmonic interference from a controllable seismic source based on relevant characteristics includes:

[0060] The amplitude acquisition module 102 obtains the amplitudes of the second-order harmonic, third-order harmonic, and fourth-order harmonic signals based on the raw data from the controllable seismic source.

[0061] The harmonic preset signal establishment module 104 obtains the waveform of each harmonic preset signal based on the scanning start and end frequencies of the fundamental wave and the amplitude of each harmonic signal.

[0062] The harmonic removal module 106, for each order harmonic preset signal waveform, obtains order harmonic related data based on the order harmonic preset signal waveform and the original data of the controllable source, and subtracts the order harmonic related data from the original data of the controllable source to obtain the controllable source data after removing the order harmonics.

[0063] The amplitudes of the second, third, and fourth harmonic signals are obtained through the following steps: extracting the force signal data from the original data of the controllable source; transforming the force signal data to the time-frequency domain; determining the amplitude of each harmonic signal in the time-frequency domain data; and obtaining the amplitudes of the second, third, and fourth harmonic signals according to the order of the harmonic signal.

[0064] The preset signal waveform for each harmonic is obtained through the following steps: the scan start and end frequencies of each harmonic signal are designed based on the scan start and end frequencies of the fundamental wave; for each harmonic, the preset signal waveform for that harmonic is obtained based on the scan start and end frequencies and amplitude values ​​of that harmonic.

[0065] The harmonic correlation data is obtained through the following steps: the harmonic preset signal waveform is correlated with the original data of the controllable source; the correlated data is converted to the time-frequency domain; and the harmonic correlation data is extracted from the converted data.

[0066] Example 3

[0067] This disclosure provides an electronic device comprising: a memory storing executable instructions; and a processor executing the executable instructions in the memory to implement the aforementioned method for suppressing controllable source harmonic interference based on relevant characteristics.

[0068] An electronic device according to an embodiment of the present disclosure includes a memory and a processor.

[0069] This 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.

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

[0071] Those skilled in the art should understand that, in order to solve the technical problem of how to obtain 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.

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

[0073] Example 4

[0074] This disclosure provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for suppressing controllable source harmonic interference based on relevant characteristics.

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

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

[0077] 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 suppressing harmonic interference of a vibroseis based on a correlation characteristic, characterized in that, The method comprises the following steps: obtaining the amplitudes of the second-order harmonic, the third-order harmonic and the fourth-order harmonic based on the original data of the controllable source; obtaining the preset signal waveforms of each order harmonic based on the scanning start and end frequencies of the fundamental wave and the amplitudes of each order harmonic signal; for each order harmonic preset signal waveform, obtaining the order harmonic related data based on the order harmonic preset signal waveform and the original data of the controllable source, subtracting the original data of the controllable source from the order harmonic related data to obtain the controllable source data after removing the order harmonic; the amplitudes of the second-order harmonic, the third-order harmonic and the fourth-order harmonic are obtained by the following steps: extracting the force signal data of the original data of the controllable source; transforming the force signal data to the time-frequency domain; determining the amplitude of each harmonic signal in the time-frequency domain data; obtaining the amplitudes of the second-order harmonic, the third-order harmonic and the fourth-order harmonic according to the order of the harmonic signal; the preset signal waveforms of each order harmonic are obtained by the following steps: designing the scanning start and end frequencies of each order harmonic signal according to the scanning start and end frequencies of the fundamental wave; for each order harmonic, obtaining the order harmonic preset signal waveform according to the scanning start and end frequencies and the amplitude value of the order harmonic; the order harmonic related data is obtained by the following steps: correlation processing is performed on the order harmonic preset signal waveform and the original data of the controllable source; the data after correlation processing is converted to the time-frequency domain; extracting the order harmonic related data in the converted data.

2. An electronic device, comprising: The electronic device comprises: a memory storing executable instructions; a processor running the executable instructions in the memory to implement the method for suppressing harmonic interference of a controllable source based on correlation characteristics according to claim 1.

3. A computer-readable storage medium, characterized in that, The computer readable storage medium stores a computer program which is executed by the processor to implement the method for suppressing harmonic interference of a controllable source based on correlation characteristics according to claim 1.

4. A device for suppressing harmonic interference of a controllable source based on a correlation characteristic, characterized in that The method comprises the following steps: an amplitude acquisition module for obtaining the amplitudes of the second-order harmonic, the third-order harmonic and the fourth-order harmonic based on the original data of the controllable source; a harmonic preset signal establishment module for obtaining the preset signal waveforms of each order harmonic based on the scanning start and end frequencies of the fundamental wave and the amplitudes of each order harmonic signal; a harmonic removal module for obtaining the order harmonic related data based on the order harmonic preset signal waveform and the original data of the controllable source for each order harmonic preset signal waveform, and subtracting the original data of the controllable source from the order harmonic related data to obtain the controllable source data after removing the order harmonic; the amplitudes of the second-order harmonic, the third-order harmonic and the fourth-order harmonic are obtained by the following steps: extracting the force signal data of the original data of the controllable source; transforming the force signal data to the time-frequency domain; determining the amplitude of each harmonic signal in the time-frequency domain data; obtaining the amplitudes of the second-order harmonic, the third-order harmonic and the fourth-order harmonic according to the order of the harmonic signal; the preset signal waveforms of each order harmonic are obtained by the following steps: designing the scanning start and end frequencies of each order harmonic signal according to the scanning start and end frequencies of the fundamental wave; For each order harmonic, according to the scanning start and end frequencies and the amplitude value of the order harmonic, an order harmonic preset signal waveform is obtained; The order harmonic related data is obtained by the following steps: The order harmonic preset signal waveform is correlated with the original data of the controllable seismic source; The data after the correlation processing is converted to a time-frequency domain; The order harmonic related data is extracted from the converted data.

Images