Seismic data crosstalk noise suppression method and device based on POCS-TVM algorithm
By converting seismic data to the frequency-wavenumber domain and performing zero-wavenumber processing using the POCS-TVM algorithm, the problem of crosstalk noise in marine seismic exploration is solved, noise suppression and effective signal recovery are achieved, and the quality of seismic data and exploration results are improved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- CHINA NAT PETROLEUM CORP
- Filing Date
- 2024-11-15
- Publication Date
- 2026-05-15
AI Technical Summary
In marine seismic exploration, crosstalk noise caused by poor electrical insulation during towed cable acquisition affects the signal-to-noise ratio and data analysis of seismic data, thus impacting exploration results.
The POCS-TVM algorithm is used to transform seismic data from the time-space domain to the frequency-wavenumber domain, generate a sampling matrix and reconstruct zero-wavenumber data, and finally transform it back to the time-space domain to remove crosstalk noise.
It effectively removes crosstalk noise, improves the signal-to-noise ratio of seismic data, eliminates artifacts in seismic data processing and interpretation, and improves the efficiency and economy of marine exploration.
Smart Images
Figure CN122043579A_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of oil and gas exploration technology, and in particular to a method and apparatus for suppressing crosstalk noise in seismic data based on the POCS-TVM algorithm. Background Technology
[0002] During marine seismic exploration, towed cable acquisition mainly involves using an air gun to generate seismic waves near the sea surface, while a geophone is connected by a cable and submerged at a depth of several to tens of meters underwater to receive seismic signals, thereby achieving marine seismic data acquisition.
[0003] During offshore towed seismic data acquisition, the variety of equipment involved, the high uncertainty of sea conditions, and various situations that may arise during construction can sometimes lead to cable leakage. This leakage is primarily caused by poor electrical insulation between different parts of the towed cable, resulting in water ingress and electrical noise, commonly known as "crosstalk noise" (similar situations can also occur in onshore acquisition). Crosstalk noise distorts seismic data and reduces the signal-to-noise ratio, making it difficult to separate the effective signal from the noise. If crosstalk noise is not suppressed, it may affect subsequent processing, interfere with the analysis and interpretation of seismic data, and ultimately impact the effectiveness of marine exploration.
[0004] Therefore, there is an urgent need for a method to suppress crosstalk noise generated during marine exploration towed cable acquisition, so as to avoid affecting the analysis results of seismic data. Summary of the Invention
[0005] This application provides a method and apparatus for suppressing crosstalk noise in seismic data based on the POCS-TVM algorithm, which can effectively remove crosstalk noise and improve the signal-to-noise ratio of seismic data.
[0006] The first aspect of this application provides a method for suppressing crosstalk noise in seismic data based on the POCS-TVM algorithm, including:
[0007] The original seismic data containing crosstalk noise is transformed from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data;
[0008] Based on the first seismic data, a sampling matrix is generated;
[0009] Multiply the first seismic data by the sampling matrix to obtain the second seismic data;
[0010] The second seismic data was reconstructed using the POCS-TVM algorithm with zero wavenumber data to obtain the third seismic data;
[0011] The third seismic data is converted from the frequency-wavenumber domain to the time-space domain to obtain seismic data after crosstalk noise suppression.
[0012] Optionally, the step of converting the original seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data includes:
[0013] The original seismic data containing crosstalk noise is transformed in one dimension in the time direction to obtain the frequency domain data volume.
[0014] A two-dimensional transformation is performed on the frequency domain data volume along the spatial direction to obtain the first seismic data.
[0015] Optionally, the step of converting the original seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data includes:
[0016] By using Fourier transform, Hartley transform, or discrete Hartley transform, the original seismic data containing crosstalk noise is transformed from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data.
[0017] Optionally, based on the first seismic data, a sampling matrix is generated, including:
[0018] Based on the number of seismic traces in the first seismic data, a unit matrix is generated;
[0019] The midpoint of the unit matrix is set to 0 to obtain the sampling matrix.
[0020] Optionally, generating a unit matrix based on the number of seismic traces from the first seismic data includes:
[0021] If the number of seismic traces in the first seismic data is odd, generate a unit matrix with the same size as the number of seismic traces in the first seismic data.
[0022] Optionally, generating a unit matrix based on the number of seismic traces from the first seismic data includes:
[0023] If the number of seismic traces in the first seismic data is even, an empty trace is added to the first seismic data to obtain a new first seismic data, and a unit matrix with the same size as the number of seismic traces in the new first seismic data is generated.
[0024] Based on the same inventive concept, a second aspect of this application provides a seismic data crosstalk noise suppression device based on the POCS-TVM algorithm, comprising:
[0025] The seismic data transformation module is used to transform the raw seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data.
[0026] A sampling matrix generation module is used to generate a sampling matrix based on the first seismic data;
[0027] The zero-wavenumber data elimination module is used to multiply the first seismic data with the sampling matrix to obtain the second seismic data;
[0028] The effective signal recovery module is used to reconstruct the second seismic data from the third seismic data using the POCS-TVM algorithm with zero wavenumber data;
[0029] The seismic data inverse transformation module is used to transform the third seismic data from the frequency-wavenumber domain to the time-space domain to obtain seismic data after crosstalk noise suppression.
[0030] Based on the same inventive concept, a third aspect of the present application provides a storage medium storing machine-executable instructions, which, when executed by a processor, implement the seismic data crosstalk noise suppression method based on the POCS-TVM algorithm proposed in the first aspect of the present application.
[0031] A fourth aspect of this application provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When executed by the processor, the program implements the seismic data crosstalk noise suppression method based on the POCS-TVM algorithm proposed in the first aspect of this application.
[0032] Compared with the prior art, this application has at least the following advantages:
[0033] The seismic data crosstalk noise suppression method based on the POCS-TVM algorithm provided in this application first transforms the original seismic data containing crosstalk noise into the frequency-wavenumber domain to obtain first seismic data. Then, based on the first seismic data, a sampling matrix is generated, and the first seismic data is multiplied by the sampling matrix to obtain second seismic data with zero wavenumber data eliminated. The second seismic data is then reconstructed using the POCS-TVM algorithm to recover the effective seismic signals lost during zero wavenumber processing, resulting in third seismic data. Finally, the third seismic data is inversely transformed into the time-space domain to obtain the seismic data with crosstalk noise suppression. Thus, by transforming the seismic data into the frequency-wavenumber domain and performing zero wavenumber processing, noise in the seismic data is suppressed, effectively removing crosstalk noise, improving the signal-to-noise ratio of the seismic data, and eliminating interpretation artifacts in seismic data processing. Attached Figure Description
[0034] Figure 1 This is a flowchart of a seismic data crosstalk noise suppression method based on the POCS-TVM algorithm in one embodiment of this application;
[0035] Figure 2 This is a schematic diagram of seismic data recorded by a single shot before crosstalk noise suppression in one embodiment of this application;
[0036] Figure 3 This is a schematic diagram of seismic data after crosstalk noise suppression from a single shot record in one embodiment of this application;
[0037] Figure 4 This is a schematic diagram of the suppression crosstalk noise recorded by a single gun in one embodiment of this application;
[0038] Figure 5 This is a schematic diagram of the functional modules of a seismic data crosstalk noise suppression device based on the POCS-TVM algorithm in one embodiment of this application;
[0039] Figure 6 This is a schematic diagram of the structure of an electronic device according to one embodiment of this application. Detailed Implementation
[0040] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.
[0041] This application proposes a method for suppressing crosstalk noise in seismic data based on the POCS-TVM algorithm. By utilizing the zero wavenumber characteristics of crosstalk noise, the noise suppression problem is transformed into a zero wavenumber reconstruction problem, thereby suppressing noise in seismic data. This method can effectively remove crosstalk noise, improve the signal-to-noise ratio of seismic data, and eliminate artifacts in seismic data processing and interpretation.
[0042] Please refer to Figure 1 , Figure 1 This is a flowchart illustrating a seismic data crosstalk noise suppression method based on the POCS-TVM algorithm, as proposed in an embodiment of this application. Figure 1 As shown, the method includes the following steps:
[0043] S101: Transform the raw seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data.
[0044] In this embodiment, the seismic data is a set of seismic waveform data arranged in time and space. This data typically contains multiple seismic traces, each corresponding to a seismic waveform record at a spatial location.
[0045] The wavenumber domain is a mathematical tool for describing wave phenomena, where the wavenumber (also known as spatial frequency) represents the number of wavelengths per unit distance. This corresponds to frequency (time frequency), but describes spatial variation rather than temporal variation. This application converts raw seismic data (including noise) to the wavenumber domain, essentially analyzing the spatial frequency distribution of seismic signals. Crosstalk noise signals will exhibit zero or very low values in the wavenumber domain.
[0046] It should be noted that the size of the first seismic data obtained after domain transformation is the same as that of the original seismic data. Domain transformation only changes the representation of the seismic data and does not change the size or integrity of the data itself. For example, if the original seismic data is an m×n matrix, the size of the first seismic data after domain transformation will still be m×n.
[0047] Furthermore, in one feasible implementation, the original seismic data containing crosstalk noise is transformed from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data, including: performing a one-dimensional transformation of the original seismic data in the time direction to obtain a frequency domain data volume; and performing a two-dimensional transformation of the frequency domain data volume along the spatial direction to obtain the first seismic data.
[0048] In this embodiment, the seismic data is processed in one-dimensional and two-dimensional directions along the time and space directions, respectively, to transform the original seismic data from the time-space domain to the frequency-wavenumber domain. This facilitates further analysis and interpretation of the seismic data and helps to better understand the propagation characteristics of seismic waves and the structure of the marine subsurface medium.
[0049] In addition, there are many methods to convert earthquake data from the time-space domain (tx domain) to the frequency-wavenumber domain (fk domain), including the Fourier transform, Hartley transform, DHT transform, etc. This application does not impose any specific restrictions on these methods.
[0050] S102: Generate a sampling matrix based on the first seismic data.
[0051] In this embodiment, the sampling matrix is a square matrix with the same size as the number of seismic traces in the first seismic data.
[0052] Furthermore, the process of generating the sampling matrix mainly includes:
[0053] S102-1: Generate a unit matrix based on the number of seismic traces in the first seismic data.
[0054] In this embodiment, it is assumed that the first seismic data is a matrix of size m×n, where m represents the number of receivers and n represents the number of time sampling points. For example, for the same seismic source, if it is received by one receiver and the receiver's time sampling frequency is once every 2 milliseconds, then the receiver will record 2000 data points in 4 seconds, and the size of the first seismic data is 1×2000. If it is received by five receivers and the receivers' time sampling frequency is once every 2 milliseconds, then each receiver will record 2000 data points in 4 seconds, and the size of the first seismic data is 5×2000.
[0055] The identity matrix is a square matrix where all elements on the diagonal from the top left to the bottom right (called the main diagonal) are 1s, and all other elements are 0s. A key characteristic of the identity matrix is that the product of any matrix and the identity matrix is the same as the matrix itself.
[0056] In this embodiment, the size of the identity matrix generated based on the first seismic data is the same as the number of seismic traces in the first seismic data. The number of seismic traces refers to the number of receivers in the same arrangement when vibrations are generated by the same source (single shot). That is, for the first seismic data of size m×n, the number of seismic traces refers to the number of receivers, which is m. Therefore, the size of the identity matrix generated based on the first seismic data is m×m.
[0057] S102-2: Set the midpoint of the unit matrix to 0 to obtain the sampling matrix.
[0058] In this embodiment, when the number of seismic traces in the first seismic data is odd, that is, when m is odd, the midpoint of the m×m identity matrix is directly set to 0, that is, the position of (m / 2+0.5, m / 2+0.5) is set to 0, thereby obtaining the sampling matrix.
[0059] Additionally, if the number of seismic traces in the first earthquake data is even, an extra trace needs to be added to obtain the adjusted first earthquake data, making the number of seismic traces in the adjusted first earthquake data odd, thus generating an identity matrix with an odd number of columns (or rows). Then, similarly, the midpoint of the identity matrix is set to 0 to obtain the sampling matrix.
[0060] S103: Multiply the first seismic data with the sampling matrix to obtain the second seismic data.
[0061] In this embodiment, the second seismic data is the data after zero wavenumber data has been removed. The first seismic data is multiplied by the sampling matrix to eliminate zero wavenumber data in the wavenumber domain, resulting in second seismic data without zero wavenumber information. During processing, since crosstalk noise has zero wavenumber characteristics, eliminating zero wavenumbers also means eliminating crosstalk noise. Therefore, by processing zero wavenumbers in the wavenumber domain, crosstalk noise is suppressed.
[0062] S104: The third earthquake data is obtained by reconstructing the second earthquake data using the POCS-TVM algorithm with zero wavenumber data.
[0063] In this embodiment, since some effective seismic signals may also be zero wavenumber in the wavenumber domain, that is, the zero wavenumber data in the first seismic data may still contain some effective seismic signals, the elimination of zero wavenumber data in the first seismic data in the manner described in step S103 above is often accompanied by the elimination of some effective signals. Therefore, the second seismic data obtained may be missing or incomplete.
[0064] To address this, this application performs zero-wavenumber data reconstruction on the second seismic data before inversely transforming it to the time-space domain, in order to recover the lost effective signals, ensure the integrity and continuity of the data, and thus improve the quality and accuracy of the seismic data.
[0065] In practice, the POCS-TVM algorithm is used to reconstruct the zero wavenumber data from the second seismic data.
[0066] S105: Convert the third earthquake data from the frequency-wavenumber domain to the time-space domain to obtain earthquake data after crosstalk noise suppression.
[0067] In this embodiment, by reconstructing the second seismic data after eliminating zero wavenumber data, more accurate and complete seismic data can be obtained. Then, by inversely transforming the reconstructed third seismic data to the time-space domain, higher quality and more effective seismic data can be obtained, which can then be used for further seismic data analysis. This can reduce the workload of subsequent geological interpretation and oil and gas resource exploration, and improve the efficiency and economy of marine seismic exploration.
[0068] This application transforms the noise suppression problem into a zero-wavenumber data reconstruction problem, effectively removing the impact of crosstalk noise on the effective signal and restoring the effective signal while removing noise, thereby improving the accuracy and completeness of seismic data. Furthermore, the method is unaffected by the magnitude or amount of crosstalk noise, effectively suppressing all types of crosstalk noise, demonstrating strong applicability, and effectively improving the signal-to-noise ratio of seismic data, eliminating artifacts in seismic data processing and interpretation.
[0069] The above-mentioned seismic data reconstruction method is illustrated below with a specific embodiment:
[0070] 1. Use the Hartley transform to convert the seismic data D containing crosstalk noise to the wavenumber domain to obtain the data H;
[0071] 2. Based on the transformed data H in step 1, generate an identity matrix I. The size of identity matrix I is m (m is an odd number) and is consistent with the number of traces in the seismic data. Set the midpoint position (m / 2+0.5, m / 2+0.5) of identity matrix I to 0 to obtain the sampling matrix I'.
[0072] 3. Multiply the wavenumber domain data H from step 1 with the sampling matrix I' generated in step 2 to obtain data H1 that does not contain 0 wavenumber information;
[0073] 4. Use the POCS-TVM algorithm to reconstruct the zero-wavenumber data from data H1 to obtain data H2;
[0074] 5. Perform the corresponding inverse transformation on the data H2 to obtain the seismic data D' after suppressing crosstalk noise.
[0075] For example, please refer to Figures 2-4 ,in, Figure 2 This is a schematic diagram of seismic data before crosstalk noise suppression from a single-shot record. Figure 3 This is a schematic diagram of seismic data after crosstalk noise suppression from a single shot record. Figure 4 It is suppressed crosstalk noise. A comparison shows that after processing by the crosstalk noise suppression method proposed in this application, Figure 2 The crosstalk noise data in the original seismic data shown is... Figure 3 The middle part has been eliminated. Figure 3 The processed seismic data shown is cleaner and smoother.
[0076] Please refer to Figure 5 Based on the same inventive concept, a second aspect of this application provides a seismic data crosstalk noise suppression device based on the POCS-TVM algorithm, comprising:
[0077] The seismic data transformation module 501 is used to transform the original seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data.
[0078] The sampling matrix generation module 502 is used to generate a sampling matrix based on the first seismic data;
[0079] The zero wavenumber data elimination module 503 is used to multiply the first seismic data with the sampling matrix to obtain the second seismic data;
[0080] The effective signal recovery module 504 is used to reconstruct the second seismic data using the POCS-TVM algorithm with zero wavenumber data to obtain the third seismic data.
[0081] The seismic data inverse transformation module 505 is used to transform the third seismic data from the frequency-wavenumber domain to the time-space domain to obtain seismic data after crosstalk noise suppression.
[0082] Optionally, the seismic data transformation module 501, which transforms the original seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data, is used for:
[0083] The original seismic data containing crosstalk noise is transformed in one dimension in the time direction to obtain the frequency domain data volume; along the spatial direction, the frequency domain data volume is transformed in two dimensions to obtain the first seismic data.
[0084] Optionally, the seismic data transformation module 501, which transforms the original seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data, is used for:
[0085] By using Fourier transform, Hartley transform, or discrete Hartley transform, the original seismic data containing crosstalk noise is transformed from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data.
[0086] Optionally, the sampling matrix generation module 502, based on the first seismic data, generates a sampling matrix for:
[0087] Based on the number of seismic traces in the first seismic data, a unit matrix is generated; the midpoint of the unit matrix is set to 0 to obtain a sampling matrix.
[0088] Optionally, the sampling matrix generation module 502, which generates a unit matrix based on the number of seismic traces in the first seismic data, is used for:
[0089] If the number of seismic traces in the first seismic data is odd, generate a unit matrix with the same size as the number of seismic traces in the first seismic data.
[0090] Optionally, the sampling matrix generation module 502, which generates a unit matrix based on the number of seismic traces in the first seismic data, is used for:
[0091] If the number of seismic traces in the first seismic data is even, an empty trace is added to the first seismic data to obtain a new first seismic data, and a unit matrix with the same size as the number of seismic traces in the new first seismic data is generated.
[0092] As the device embodiment is basically similar to the method embodiment, the description is relatively simple, and relevant parts can be found in the description of the method embodiment.
[0093] Based on the same inventive concept, embodiments of this application provide a storage medium storing machine-executable instructions. When the machine-executable instructions are executed by a processor, they implement the seismic data crosstalk noise suppression method based on the POCS-TVM algorithm proposed in the first aspect of this application.
[0094] It should be noted that the specific implementation of the storage medium in this application embodiment refers to the specific implementation of the seismic data crosstalk noise suppression method based on the POCS-TVM algorithm proposed in the foregoing embodiments of this application, and will not be repeated here.
[0095] Based on the same inventive concept, and referring to Figure 6 This application provides an electronic device 600, including a processor 601 and a memory 602. The memory 602 stores machine-executable instructions that can be executed by the processor 601. The processor 601 executes the machine-executable instructions to implement the seismic data crosstalk noise suppression method based on the POCS-TVM algorithm proposed in this application.
[0096] It should be noted that the specific implementation of the electronic device 600 in this application embodiment refers to the specific implementation of the seismic data crosstalk noise suppression method based on the POCS-TVM algorithm proposed in the foregoing embodiments of this application, and will not be repeated here.
[0097] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other.
[0098] Those skilled in the art will understand that embodiments of this application can be provided as methods, apparatus, or computer program products. Therefore, embodiments of this application can take the form of entirely hardware embodiments, entirely software embodiments, or embodiments combining software and hardware aspects. Furthermore, embodiments of this application can take the form of computer program products implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0099] This application describes embodiments with reference to flowchart illustrations and / or block diagrams of methods, terminal devices (systems), and computer program products according to embodiments of this application. It should be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing terminal device to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal device, generate instructions for implementing the flowchart illustrations. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.
[0100] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing terminal device to operate in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.
[0101] These computer program instructions can also be loaded onto a computer or other programmable data processing terminal equipment, causing a series of operational steps to be performed on the computer or other programmable terminal equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable terminal equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.
[0102] Although preferred embodiments of the present application have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of the present application.
[0103] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or terminal device that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or terminal device. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or terminal device that includes said element.
[0104] The crosstalk noise suppression method, apparatus, medium, and device provided in this application have been described in detail above. Specific examples have been used to illustrate the principles and implementation methods of this application. The description of the above embodiments is only for the purpose of helping to understand the method and core ideas of this application. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of this application. Therefore, the content of this specification should not be construed as a limitation of this application.
Claims
1. A method for suppressing crosstalk noise in seismic data based on the POCS-TVM algorithm, characterized in that, include: The original seismic data containing crosstalk noise is transformed from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data; Based on the first seismic data, a sampling matrix is generated; Multiply the first seismic data by the sampling matrix to obtain the second seismic data; The second seismic data was reconstructed using the POCS-TVM algorithm with zero wavenumber data to obtain the third seismic data; The third seismic data is converted from the frequency-wavenumber domain to the time-space domain to obtain seismic data after crosstalk noise suppression.
2. The method according to claim 1, characterized in that, The process of converting the raw seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data includes: The original seismic data containing crosstalk noise is transformed in one dimension in the time direction to obtain the frequency domain data volume. A two-dimensional transformation is performed on the frequency domain data volume along the spatial direction to obtain the first seismic data.
3. The method according to claim 1, characterized in that, The process of converting the raw seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data includes: By using Fourier transform, Hartley transform, or discrete Hartley transform, the original seismic data containing crosstalk noise is transformed from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data.
4. The method according to claim 1, characterized in that, The step of generating a sampling matrix based on the first seismic data includes: Based on the number of seismic traces in the first seismic data, a unit matrix is generated; The midpoint of the unit matrix is set to 0 to obtain the sampling matrix.
5. The method according to claim 4, characterized in that, The generation of a unit matrix based on the number of seismic traces from the first seismic data includes: If the number of seismic traces in the first seismic data is odd, generate a unit matrix with the same size as the number of seismic traces in the first seismic data.
6. The method according to claim 4, characterized in that, The generation of a unit matrix based on the number of seismic traces from the first seismic data includes: If the number of seismic traces in the first seismic data is even, an empty trace is added to the first seismic data to obtain a new first seismic data, and a unit matrix with the same size as the number of seismic traces in the new first seismic data is generated.
7. A seismic data crosstalk noise suppression device based on the POCS-TVM algorithm, characterized in that, include: The seismic data transformation module is used to transform the raw seismic data containing crosstalk noise from the time-space domain to the frequency-wavenumber domain to obtain the first seismic data. A sampling matrix generation module is used to generate a sampling matrix based on the first seismic data; The zero-wavenumber data elimination module is used to multiply the first seismic data with the sampling matrix to obtain the second seismic data; The effective signal recovery module is used to reconstruct the second seismic data from the third seismic data using the POCS-TVM algorithm with zero wavenumber data; The seismic data inverse transformation module is used to transform the third seismic data from the frequency-wavenumber domain to the time-space domain to obtain seismic data after crosstalk noise suppression.
8. The apparatus according to claim 7, characterized in that, The sampling matrix generation module, based on the first seismic data, generates a sampling matrix for the following purposes: Based on the number of seismic traces in the first seismic data, a unit matrix is generated; the midpoint of the unit matrix is set to 0 to obtain a sampling matrix.
9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the seismic data crosstalk noise suppression method based on the POCS-TVM algorithm as described in any one of claims 1 to 6.
10. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the seismic data crosstalk noise suppression method based on the POCS-TVM algorithm as described in any one of claims 1 to 6.