Crosstalk noise suppression method and device, medium and equipment
By converting seismic data to the frequency-wavenumber domain to eliminate zero wavenumber data and reconstructing effective signals, the problem of crosstalk noise in marine seismic exploration is solved, thereby improving the signal-to-noise ratio and exploration effect of seismic data.
Patent Information
- Application Number
- CN202410992662.5
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-07-23
- Publication Date
- 2026-01-23
AI Technical Summary
In marine seismic exploration, crosstalk noise caused by cable leakage affects the signal-to-noise ratio of seismic data, leading to data distortion and difficulty in signal separation, which in turn affects subsequent processing and exploration results.
The original seismic data is transformed from the time-space domain to the frequency-wavenumber domain, a sampling matrix is generated and zero wavenumber data is eliminated, the effective signal is reconstructed through Shearlet transform, and finally the data is inversely transformed back to the time-space domain to suppress crosstalk noise.
It effectively removes the influence of crosstalk noise, improves the signal-to-noise ratio of seismic data, eliminates artifacts in seismic data processing and interpretation, and improves exploration efficiency and accuracy.
Smart Images
Figure CN121386006A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of oil and gas exploration, and particularly relates to a crosstalk noise suppression method and device, a medium and equipment. BACKGROUND
[0002] In the process of marine seismic exploration, the towed cable acquisition mainly uses air guns to excite seismic waves near the sea surface, and at the same time, the geophones are connected by cables and sunk to a position several meters to tens of meters underwater to receive seismic signals, thereby realizing offshore seismic data acquisition.
[0003] In the process of offshore towed cable acquisition, due to the variety of equipment involved in seismic data acquisition, the great uncertainty of marine conditions, and various situations that may be encountered during construction, the cable may sometimes have a leakage situation. This leakage situation is mainly caused by water entering due to poor electrical insulation between parts of the towed cable, thereby causing electrical noise, which is usually referred to as "crosstalk noise" (in fact, similar situations may also occur in land acquisition). Crosstalk noise can cause distortion of seismic data and reduce the signal-to-noise ratio, making it difficult to separate effective signals from noise. If crosstalk noise is not suppressed, it may affect the normal progress of subsequent processing, interfere with the analysis and interpretation of seismic data, and ultimately affect the effect of marine exploration.
[0004] Therefore, there is an urgent need for a method to suppress the crosstalk noise generated in the process of marine exploration towed cable acquisition to avoid affecting the analysis results of seismic data. SUMMARY
[0005] The present application provides a crosstalk noise suppression method, device, medium and equipment, which can effectively remove the influence of crosstalk noise on effective signals and improve the signal-to-noise ratio of seismic data.
[0006] The first aspect of the embodiment of the present application provides a crosstalk noise suppression method, and the method comprises:
[0007] Pretreating original seismic data containing crosstalk noise to obtain first seismic data; the first seismic data is data obtained after converting the original seismic data from the time-space domain to the frequency-wave number domain;
[0008] Generating a sampling matrix based on the first seismic data;
[0009] Processing the first seismic data based on the sampling matrix to obtain second seismic data; the second seismic data is data after eliminating zero wave number data;
[0010] Reconverting the second seismic data to the time-space domain to obtain target seismic data after crosstalk noise suppression.
[0011] Optionally, based on the first seismic data, a sampling matrix is generated, comprising:
[0012] Based on the first seismic data, a unit matrix is generated; the size of the unit matrix is the same as the number of seismic traces in the first seismic data;
[0013] The midpoint position of the unit matrix is set to 0 to obtain the sampling matrix.
[0014] Optionally, based on the first seismic data, a unit matrix is generated, comprising:
[0015] In the case where the number of seismic traces in the first seismic data is even, a new empty trace data is added to obtain adjusted first seismic data;
[0016] Based on the adjusted first seismic data, the unit matrix is generated.
[0017] Optionally, based on the sampling matrix, the first seismic data is processed to obtain second seismic data, comprising:
[0018] The first seismic data is multiplied by the sampling matrix to obtain the second seismic data.
[0019] Optionally, the second seismic data is inverse transformed to the time-space domain to obtain target seismic data after cross-talk noise suppression, comprising:
[0020] The second seismic data is zero wave number data reconstructed to obtain reconstructed second seismic data;
[0021] The reconstructed second seismic data is inverse transformed to the time-space domain to obtain the target seismic data after cross-talk noise suppression.
[0022] Optionally, the original seismic data containing cross-talk noise is preprocessed to obtain first seismic data, comprising:
[0023] The original seismic data is one-dimensional Hartley transformed in the time direction to obtain a frequency domain data body;
[0024] Along the spatial direction, the frequency domain data body is two-dimensional Hartley transformed to obtain the first seismic data.
[0025] Based on the same inventive concept, the second aspect of the embodiments of the present application provides a cross-talk noise suppression device, the device comprising:
[0026] A wave number domain transformation module is configured to preprocess original seismic data containing cross-talk noise to obtain first seismic data; the first seismic data is data obtained by converting the original seismic data from the time-space domain to the frequency-wave number domain.
[0027] a sampling module, configured to generate a sampling matrix based on the first seismic data;
[0028] a cancellation module, configured to process the first seismic data based on the sampling matrix to obtain second seismic data; the second seismic data is data after zero wavenumber data is eliminated;
[0029] a reverse transformation module, configured to reversely transform the second seismic data to a time-space domain to obtain target seismic data after crosstalk noise is suppressed.
[0030] Based on the same inventive concept, a storage medium is provided in a third aspect of embodiments of the present application. The storage medium stores machine executable instructions. The machine executable instructions are executed by a processor to implement the crosstalk noise suppression method according to the first aspect of the present application.
[0031] A fourth aspect of embodiments of the present application provides an electronic device. The electronic device includes a memory, a processor, and a computer program stored in the memory and executable on the processor. The processor implements the crosstalk noise suppression method according to the first aspect of the present application when executed.
[0032] Compared with the prior art, the present application has the following advantages:
[0033] The crosstalk noise suppression method provided by the embodiments of the present application first pre-processes original seismic data containing crosstalk noise to obtain first seismic data. The first seismic data is data obtained after the original seismic data is converted from a time-space domain to a frequency-wavenumber domain. Then, based on the first seismic data, a sampling matrix is generated, and the first seismic data is processed based on the sampling matrix to obtain second seismic data. The second seismic data is data after zero wavenumber data is eliminated. Finally, the second seismic data is reversely transformed to a time-space domain to obtain target seismic data after crosstalk noise is suppressed. Thus, by converting seismic data to a frequency-wavenumber domain for zero wavenumber processing, the suppression of noise in seismic data is realized, which can effectively remove the influence of crosstalk noise on effective signals, improve the signal-to-noise ratio of seismic data, and eliminate false images in seismic data processing and interpretation. BRIEF DESCRIPTION OF DRAWINGS
[0034] Figure 1 is a flowchart of a crosstalk noise suppression method in an embodiment of the present application;
[0035] Figure 2 is a schematic diagram of seismic data before crosstalk noise suppression for a single shot record in an embodiment of the present application;
[0036] Figure 3 is a schematic diagram of seismic data after crosstalk noise suppression for a single shot record in an embodiment of the present application;
[0037] Figure 4 is a functional module schematic diagram of a crosstalk noise suppression device in an embodiment of the present application;
[0038] Figure 5 is a structural schematic diagram of an electronic device in an embodiment of the present application. DETAILED DESCRIPTION
[0039] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the scope of the present application.
[0040] Ocean seismic exploration plays an important role in oil and gas exploration. The towed cable acquisition technology is mature and relatively low in cost, and is the main acquisition method for offshore exploration at present. Towed cable acquisition mainly uses air guns to excite seismic waves near the sea surface, and simultaneously connects the geophones through cables to be sunk to a position several meters to tens of meters underwater to receive seismic signals, so as to realize offshore seismic data acquisition.
[0041] In the process of offshore towed cable acquisition, due to the variety of equipment involved in seismic data acquisition, the great uncertainty of ocean conditions, and various situations that may be encountered in the process of construction, the cable may sometimes have a leakage situation. This leakage situation is mainly caused by water entering due to poor electrical insulation between parts of the towed cable, thereby causing electrical noise, which is usually referred to as "crosstalk noise" (in fact, similar situations may also be encountered in land acquisition). This noise will propagate along the cable at the speed of light into other channels, and present one or more flat events on different offsets of the received seismic single-shot record. The amplitude energy and phase frequency of these events are very similar, and if not suppressed, will interfere with subsequent seismic data processing, easily causing processing artifacts, directly affecting the analysis and judgment of the underground formation by the seismic interpretation personnel. If the noise is serious, it may even affect the final exploration results such as well location determination. Therefore, it is necessary to suppress the crosstalk noise generated in the towed cable acquisition of ocean exploration.
[0042] In view of this, the present application proposes a crosstalk noise suppression method, which uses the zero wavenumber characteristics of crosstalk noise to convert the noise suppression problem into a zero wavenumber reconstruction problem, to realize the suppression of noise in seismic data, effectively remove the influence of crosstalk noise on effective signals, improve the signal-to-noise ratio of seismic data, and eliminate the interpretation artifacts of seismic data processing.
[0043] Please refer to Figure 1 , Figure 1is a flow chart of a crosstalk noise suppression method according to an embodiment of the present application. As shown in Figure 1 the method comprises the following steps:
[0044] S101: preprocessing the original seismic data containing crosstalk noise to obtain first seismic data.
[0045] In the embodiment, the seismic data is a set of seismic waveform data arranged in time and space. These data usually contain multiple seismic traces, each corresponding to a seismic waveform record at a spatial location. Preprocessing the original seismic data containing crosstalk noise means converting the original seismic data from the time-space domain (t-x domain) to the frequency-wavenumber domain (f-k domain).
[0046] The wavenumber domain is a mathematical tool for describing wave phenomena, in which the wavenumber (also known as spatial frequency) represents the number of wavelengths per unit distance. This corresponds to frequency (temporal frequency), but describes spatial variation rather than temporal variation. The present application converts the original seismic data (including noise) to the wavenumber domain, which is actually an analysis of the spatial frequency distribution of the seismic signal. For crosstalk noise signals, they will exhibit zero or very low values in the wavenumber domain.
[0047] The first seismic data obtained by preprocessing in the embodiment is the data obtained after converting the original seismic data from the time-space domain to the frequency-wavenumber domain. 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 form of the seismic data, and does not change the size and integrity of the data itself. For example, if the original seismic data is an m x n matrix, the size of the first seismic data after domain transformation is still m x n.
[0048] Further, in a feasible implementation, preprocessing the original seismic data containing crosstalk noise to obtain first seismic data comprises: performing one-dimensional Hartley transform in the time direction on the original seismic data to obtain a frequency domain data body; and performing two-dimensional Hartley transform on the frequency domain data body along the spatial direction to obtain the first seismic data.
[0049] In the embodiment, two-dimensional processing of the seismic data is performed along the time direction and the spatial direction respectively using Hartley transform to convert the original seismic data from the time-space domain to the frequency-wavenumber domain, thereby facilitating further analysis and interpretation of the seismic data, and helping to better understand the propagation characteristics of seismic waves and the structure of the ocean underground medium.
[0050] In addition, there are many methods for converting the time-space domain to the frequency-wavenumber domain, including Fourier transform, Hartley transform, DHT transform, etc., and the present application does not make specific limitations thereto.
[0051] S102: generating a sampling matrix based on the first seismic data.
[0052] In the present embodiment, the sampling matrix is a square matrix with the same size as the number of seismic traces in the first seismic data.
[0053] Further, the process of generating the sampling matrix mainly includes:
[0054] S102-1: generating an identity matrix based on the first seismic data; the identity matrix has the same size as the number of seismic traces in the first seismic data.
[0055] In the present embodiment, it is assumed that the first seismic data is a matrix with a size of m x n, where m represents the number of geophones, and n represents the number of time sampling points. For example, for the same seismic source, if it is received by 1 geophone, and the time sampling frequency of the geophone is 2 milliseconds, then within 4 seconds, the geophone will record 2000 data, and at this time, the size of the first seismic data is 1 x 2000. If it is received by 5 geophones, and the time sampling frequency of the geophone is 2 milliseconds, then within 4 seconds, each geophone will record 2000 data, and at this time, the size of the first seismic data is 5 x 2000.
[0056] The identity matrix is a square matrix, and the elements on the diagonal line from the upper left corner to the lower right corner (referred to as the main diagonal line) are all 1, and all other elements are 0. The identity matrix has the characteristic that any matrix multiplied by the identity matrix is equal to itself.
[0057] The identity matrix generated based on the first seismic data in the present embodiment has the same size as the number of seismic traces in the first seismic data. The number of seismic traces refers to the number of geophones in the same arrangement when the same seismic source (single shot) is excited to vibrate. That is, for the first seismic data with a size of m x n described above, the number of seismic traces refers to the number of geophones, that is, m. Therefore, the size of the identity matrix generated based on the first seismic data is m x m.
[0058] S102-2: setting the midpoint position of the identity matrix to 0 to obtain the sampling matrix.
[0059] In the present embodiment, in the case where the number of seismic traces in the first seismic data is odd, that is, m is odd, the midpoint position of the identity matrix with a size of m x m 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.
[0060] In addition, in the case that the number of seismic traces in the first seismic data is even, a new empty trace data is added to obtain adjusted first seismic data, so that the number of seismic traces in the adjusted first seismic data is odd, thereby generating a unit matrix with an odd number of columns (or rows). Then, the midpoint position of the unit matrix is set to 0 to obtain the sampling matrix.
[0061] S103: processing the first seismic data based on the sampling matrix to obtain second seismic data.
[0062] In the embodiment, the second seismic data is data after the zero wave number data is eliminated. Processing the first seismic data based on the sampling matrix means multiplying the first seismic data by the sampling matrix, thereby eliminating the zero wave number data in the wave number domain to obtain the second seismic data without zero wave number information. During the processing, since the cross-talk noise has the zero wave number characteristic, the elimination of the zero wave number also means the elimination of the cross-talk noise, thereby realizing the suppression of the cross-talk noise by processing the zero wave number in the wave number domain.
[0063] S104: inverse transforming the second seismic data to the time-space domain to obtain target seismic data after the cross-talk noise is suppressed.
[0064] In the embodiment, after obtaining the second seismic data after the zero wave number data is eliminated, the second seismic data is inverse transformed from the frequency-wave number domain to the time-space domain by the method in step S101, thereby obtaining the effective seismic data after the cross-talk noise is suppressed, which is used for marine geological structure analysis and the like.
[0065] Further, the process mainly includes:
[0066] S104-1: reconstructing the zero wave number data of the second seismic data to obtain reconstructed second seismic data.
[0067] In the embodiment, since some effective seismic signals may also exhibit zero wave number in the wave number domain, that is, the zero wave number data in the first seismic data may also contain some effective seismic signals, and thus the elimination of the zero wave number data in the first seismic data by the method in step S103 often accompanies the elimination of some effective signals, and thus the obtained second seismic data may be incomplete data.
[0068] In this regard, the application reconstructs the zero wave number data of the second seismic data before inverse transforming the second seismic data to the time-space domain, so as to recover the lost effective signals, ensure the integrity and continuity of the data, and further improve the quality and accuracy of the seismic data.
[0069] In a specific implementation, the missing effective data can be reconstructed by using a Shearlet transform method. The Shearlet transform can decompose the seismic data into Shearlet functions of different scales and directions, each of which is described by a scale parameter and a direction parameter. The scale parameter determines the size of the Shearlet function, and the direction parameter determines the direction thereof. Then, the seismic data can be reconstructed by linearly combining the Shearlet functions.
[0070] S104-2: inverse transform the reconstructed second seismic data to the time-space domain to obtain target seismic data after cross-talk noise suppression.
[0071] In this embodiment, the second seismic data after the elimination of the zero wave number data can be reconstructed to obtain more accurate and complete seismic data. Then, the reconstructed second seismic data can be inverse transformed to the time-space domain to obtain high-quality and effective seismic data, which can be used for further seismic data analysis, thereby reducing the workload of subsequent geological interpretation and oil and gas resource exploration and improving the efficiency and economy of marine seismic exploration.
[0072] The present application can effectively remove the influence of cross-talk noise on effective signals by converting the noise suppression problem into a zero wave number data reconstruction problem, and can restore the effective signals while removing the noise, so as to improve the accuracy and completeness of the seismic data. Moreover, the above method can effectively suppress cross-talk noise regardless of the size and amount of the cross-talk noise, has strong applicability, can effectively improve the signal-to-noise ratio of the seismic data, and can eliminate the false image of seismic data processing and interpretation.
[0073] The above seismic data reconstruction method will be described below in combination with a specific embodiment.
[0074] 1. The seismic data D containing cross-talk noise is converted to the wave number domain by using the Hartley transform to obtain data H.
[0075] 2. According to the transformed data H in 1, a unit matrix I is generated, the size of the unit matrix I is m (m is an odd number), and the number of traces of the seismic data is consistent, and the midpoint position (m / 2+0.5, m / 2+0.5) of the unit matrix I is set to 0 to obtain a sampling matrix I';
[0076] 3. The wave number domain data H in 1 is multiplied by the sampling matrix I' generated in 2 to obtain data H1 containing no zero wave number information.
[0077] 4. The data H1 is reconstructed by using a Shearlet data reconstruction method to obtain data H2.
[0078] 5. Perform corresponding inverse transform on the data H2 to obtain the crosstalk noise suppressed seismic data D'.
[0079] For example, refer to Figure 2 and Figure 3 wherein, Figure 2 is a schematic diagram of the crosstalk noise suppressed seismic data of a single shot record. Figure 3 is a schematic diagram of the crosstalk noise suppressed seismic data of a single shot record. It can be found by comparison that the crosstalk noise data concentrated at about 900 ms in the original seismic data shown in Figure 2 is eliminated in Figure 3 , and Figure 3 the processed seismic data shown in is cleaner and smoother.
[0080] For example, refer to Figure 4 Based on the same inventive concept, a second aspect of the embodiments of the present application provides a crosstalk noise suppression device. The crosstalk noise suppression device 400 comprises:
[0081] A wave number domain transform module 401 is configured to perform preprocessing on original seismic data containing crosstalk noise to obtain first seismic data. The first seismic data is data obtained by converting the original seismic data from a time-space domain to a frequency-wave number domain.
[0082] A sampling module 402 is configured to generate a sampling matrix based on the first seismic data.
[0083] An elimination module 403 is configured to process the first seismic data based on the sampling matrix to obtain second seismic data. The second seismic data is data after the zero wave number data is eliminated.
[0084] An inverse transform module 404 is configured to perform inverse transform on the second seismic data to a time-space domain to obtain target seismic data after crosstalk noise suppression.
[0085] Optionally, the sampling module 402 comprises:
[0086] A unit matrix generation submodule is configured to generate a unit matrix based on the first seismic data. The size of the unit matrix is the same as the number of seismic traces in the first seismic data.
[0087] A sampling matrix generation submodule is configured to set the midpoint position of the unit matrix to 0 to obtain the sampling matrix.
[0088] Optionally, the unit matrix generation submodule is specifically configured to:
[0089] In the case that the number of seismic traces in the first seismic data is even, a new empty trace data is added to obtain adjusted first seismic data.
[0090] generate a unit matrix based on the adjusted first seismic data.
[0091] Optionally, the eliminating module 403 is specifically configured to:
[0092] multiply the first seismic data with the sampling matrix to obtain second seismic data.
[0093] Optionally, the inverse transform module 404 comprises:
[0094] a reconstruction sub-module, configured to perform zero wavenumber data reconstruction on the second seismic data to obtain reconstructed second seismic data;
[0095] an inverse transform sub-module, configured to inverse transform the reconstructed second seismic data to a time-space domain to obtain target seismic data after crosstalk noise suppression.
[0096] Optionally, the wavenumber domain transform module 401 comprises:
[0097] a first transform sub-module, configured to perform one-dimensional Hartley transform in a time direction on original seismic data to obtain a frequency domain data body;
[0098] a second transform sub-module, configured to perform two-dimensional Hartley transform on the frequency domain data body along a space direction to obtain the first seismic data.
[0099] For the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the related parts refer to the part of the method embodiment.
[0100] In a third aspect, based on the same inventive concept, the embodiments of the present application provide a storage medium, the storage medium storing machine executable instructions, the machine executable instructions being executed by a processor to implement the crosstalk noise suppression method according to the first aspect of the present application.
[0101] It should be noted that the specific implementation of the storage medium of the embodiments of the present application is described with reference to the specific implementation of the crosstalk noise suppression method according to the first aspect of the embodiments of the present application, which will not be repeated here.
[0102] In a fourth aspect, based on the same inventive concept, with reference to Figure 5 The embodiments of the present application provide an electronic device 500, comprising a processor 501 and a memory 502; the memory 502 stores machine executable instructions executable by the processor 501, and the processor 501 is configured to execute the machine executable instructions to implement the crosstalk noise suppression method according to the first aspect of the present application.
[0103] It should be noted that the specific implementation of the electronic device 500 in the embodiments of the present application is described with reference to the specific implementation of the crosstalk noise suppression method in the first aspect of the embodiments of the present application, which will not be described herein again.
[0104] Each of the embodiments in the present specification is described in a progressive manner, and each embodiment focuses on the difference from other embodiments. The same or similar parts between the embodiments can be referred to each other.
[0105] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, device, or computer program product. Therefore, the embodiments of the present application can be in the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the embodiments of the present application can be in the form of a computer program product 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.
[0106] The embodiments of the present application are described with reference to flowcharts and / or block diagrams according to the method, terminal device (system), and computer program product of the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of the flows and / or blocks in the flowcharts and / or block diagrams can be implemented by computer program instructions. These computer program instructions can be provided to a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing terminal device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal device produce a device for implementing the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one flow or multiple flows and / or blocks Figure 1 The functions specified in one flow or multiple flows and / or blocks
[0107] These computer program instructions can also be stored in a computer readable storage medium that can guide the computer or other programmable data processing terminal device to work in a specific way, so that the instructions stored in the computer readable storage medium produce a product including instruction devices that implement the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one flow or multiple flows and / or blocks Figure 1 The functions specified in one flow or multiple flows and / or blocks
[0108] These computer program instructions can also be loaded into a computer or other programmable data processing terminal device, so that a series of operation steps are performed on the computer or other programmable terminal device to produce a computer implemented process, so that the instructions executed on the computer or other programmable terminal device provide a device for implementing the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one flow or multiple flows and / or blocks Figure 1steps of the method.
[0109] While the preferred embodiments of the application have been described above, it should be understood that many modifications and adaptations will occur to those skilled in the art upon the reading and understanding of the foregoing description. For example, the embodiments of the application can be applied to any type of computerized system, not just the one described above. Therefore, it is intended that the claims be construed as including all such modifications and adaptations as fall within the scope of the embodiments of the application.
[0110] Finally, it should be noted that, in the description above, relative terms such as first and second, etc. are used merely to distinguish one entity or action from another, without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises", "comprising", or any other variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but can include other elements not expressly listed or inherent to such process, method, article, or apparatus. An element proceeded by "comprises... a" does not, without more constraints, exclude the existence of additional identical elements in the process, method, article, or apparatus that comprises the element.
[0111] The above provides a crosstalk noise suppression method, device, medium and equipment, and the principles and implementation manners of the application are described by using specific examples. The above description of the embodiments is only used to help understand the method and core idea of the application; meanwhile, for those skilled in the art, according to the idea of the application, the specific implementation manners and application scope can be changed, and the above description should not be understood as limiting the application.
Claims
1. A method of crosstalk noise suppression, the method comprising: The method comprises: preprocessing original seismic data containing crosstalk noise to obtain first seismic data; the first seismic data is data obtained after converting the original seismic data from a time-space domain to a frequency-wave number domain; generating a sampling matrix based on the first seismic data; processing the first seismic data based on the sampling matrix to obtain second seismic data; the second seismic data is data after zero wave number data is eliminated; inverse transforming the second seismic data to a time-space domain to obtain target seismic data after crosstalk noise suppression.
2. The method of claim 1, wherein, The method comprises: generating a unit matrix based on the first seismic data; the size of the unit matrix is the same as the number of seismic traces in the first seismic data; setting the midpoint position of the unit matrix to 0 to obtain the sampling matrix.
3. The method of claim 2, wherein, The method comprises: in the case that the number of seismic traces in the first seismic data is even, adding a new empty trace data to obtain adjusted first seismic data; generating the unit matrix based on the adjusted first seismic data.
4. The method of claim 1, wherein, The method comprises: multiplying the first seismic data and the sampling matrix to obtain the second seismic data.
5. The method of claim 1, wherein, The method comprises: reconstructing zero wave number data of the second seismic data to obtain reconstructed second seismic data; inverse transforming the reconstructed second seismic data to a time-space domain to obtain the target seismic data after crosstalk noise suppression.
6. The method of claim 1, wherein, The method comprises: performing one-dimensional Hartley transform in the time direction on the original seismic data to obtain a frequency domain data body; performing two-dimensional Hartley transform on the frequency domain data body along the space direction to obtain the first seismic data.
7. A crosstalk noise suppression device characterized by comprising: The device comprises: a wave number domain conversion module configured to preprocess original seismic data containing crosstalk noise to obtain first seismic data; the first seismic data is data obtained after converting the original seismic data from a time-space domain to a frequency-wave number domain; a sampling module configured to generate a sampling matrix based on the first seismic data; an elimination module configured to process the first seismic data based on the sampling matrix to obtain second seismic data; the second seismic data is data after zero wave number data is eliminated; an inverse transformation module configured to inverse transform the second seismic data to a time-space domain to obtain target seismic data after crosstalk noise suppression.
8. The apparatus of claim 7, wherein, The sampling module comprises: a unit matrix generation submodule configured to generate a unit matrix based on the first seismic data; the size of the unit matrix is the same as the number of seismic traces in the first seismic data; a sampling matrix generation submodule configured to set the midpoint position of the unit matrix to 0 to obtain the sampling matrix.
9. A computer readable storage medium having stored thereon a computer program, characterized in that, The computer program, when executed by a processor, implements the crosstalk noise suppression method according to any one of claims 1 to 7.
10. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor executes the computer program to implement the crosstalk noise suppression method according to any one of claims 1 to 7.