Method, device, equipment and medium for separating uplink and downlink wave fields of seabed node data

Abstract

The present invention relates to the field of geophysical signal processing, and proposes a method, device, equipment, and medium for separating upgoing and downgoing wave fields of seabed node data, wherein the method includes: based on the seabed node, first data and second data are collected by detectors in different media and respectively subjected to dual-tree complex wavelet transform to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum, and second complex wavelet coefficients and their corresponding second amplitude spectrum, wherein the signal-to-noise ratio of the first data is higher than the signal-to-noise ratio of the second data; calculating a matching coefficient that matches the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum; updating the second complex wavelet coefficient based on the matching coefficient, and performing a dual-tree complex wavelet inverse transform based on the updated second complex wavelet coefficient to obtain third data; calculating the upgoing wave field and the downgoing wave field based on the first data and the third data and separating them. The solution disclosed by the present invention improves the separation effect of the upgoing wave field and the downgoing wave field, and has a significant effect on suppressing ghost waves.

Description

Technical Field

[0001] The present invention relates to the field of geophysical signal processing, and in particular to a method, device, equipment and medium for separating uplink and downlink wave fields of seabed node data. Background Art

[0002] In ocean bottom node (OBN) seismic exploration, the sea surface is a strongly reflecting interface, which can generate ghost waves in the data received at the seafloor. Ghost wave interference creates false events in seismic records and creates notches in the spectrum, reducing the resolution of effective underground reflection seismic signals and affecting the accuracy of seismic data interpretation.

[0003] For seismic data collected by OBN (onshore navigation network) geophones and land geophones, the effective underground reflection signals have the same polarity in both the geophones, while ghost waves have opposite polarity. Ghost waves also cause complementary notches. Therefore, seismic data from both geophones are typically processed to generate upgoing and downgoing wavefields. Ghost waves are suppressed by separating these wavefields.

[0004] Among them, after matching and calibrating the seismic data received by the underwater geophone and the land geophone, the upgoing wave field and the downgoing wave field are calculated. The matching calibration factor is obtained within the window containing only the first reflection wave or the seabed refraction wave, which is not applicable to the overlapping of the first reflection wave and the ghost wave in shallow water; the matching calibration is performed according to the analytical relationship between the water detection pressure field and the land detection velocity field in the elastic model, which is not applicable to the case where the seawater density, velocity and detection point coupling are space-varying; the matching operator is obtained after cross-ghosting, and the ghost wave extension operator is derived based on the assumption of the horizontal seabed, which is not applicable to complex seabed conditions. Based on the problems existing in the above existing technologies, the seismic data received by the underwater geophone and the land geophone cannot be reasonably matched, resulting in the inability to effectively separate the upgoing and downgoing waves, which is not conducive to the suppression of ghost waves. Summary of the Invention

[0005] In view of this, the present invention proposes a method, device, equipment and medium for uplink and downlink wave field separation of seabed node data, which at least solves the problem in the related technologies of uplink and downlink wave field separation that the seismic data received by underwater detectors and land detectors cannot be reasonably matched, resulting in the inability to effectively separate the uplink and downlink waves, which is not conducive to the suppression of ghost waves.

[0006] Based on the above objectives, one aspect of an embodiment of the present invention provides a method for separating upgoing and downgoing wavefields of seabed node data, comprising: collecting first data and second data by detectors in different media based on the seabed node and performing dual-tree complex wavelet transforms on them respectively to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum, and second complex wavelet coefficients and their corresponding second amplitude spectrum, wherein the signal-to-noise ratio of the first data is higher than the signal-to-noise ratio of the second data; calculating a matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum; updating the second complex wavelet coefficients based on the matching coefficients, and performing inverse dual-tree complex wavelet transform on the updated second complex wavelet coefficients to obtain third data; and calculating and separating the upgoing and downgoing wavefields based on the first and third data.

[0007] In some embodiments, the step of calculating the matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum includes: calculating the amplitude value ratio of the first amplitude spectrum and the second amplitude spectrum at corresponding frequency positions, so as to obtain the matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum at the corresponding frequency position based on the amplitude value ratio; and traversing each corresponding frequency position of the first amplitude spectrum and the second amplitude spectrum to obtain each matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum at each corresponding frequency position.

[0008] In some embodiments, the step of updating the second complex wavelet coefficient based on the matching coefficient and performing a dual-tree complex wavelet inverse transform based on the updated second complex wavelet coefficient to obtain the third data includes: updating the second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum based on each matching coefficient; performing a dual-tree complex wavelet inverse transform based on the updated second complex wavelet coefficient corresponding to each corresponding frequency position to obtain the third data.

[0009] In some embodiments, the step of updating the second complex wavelet coefficient corresponding to the second amplitude spectrum at each corresponding frequency position based on each matching coefficient includes: multiplying the second complex wavelet coefficient corresponding to the second amplitude spectrum at each corresponding frequency position by each matching coefficient corresponding to each corresponding frequency position to obtain the updated second complex wavelet coefficient corresponding to each corresponding frequency position.

[0010] In some embodiments, the step of multiplying the second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum by each matching coefficient corresponding to each corresponding frequency position includes: multiplying the real part and imaginary part of the second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum by each matching coefficient corresponding to each corresponding frequency position, respectively.

[0011] In some embodiments, the step of calculating and separating an upgoing wavefield and a downgoing wavefield based on the first data and the third data includes: calculating an upgoing wavefield based on the sum of the first data and the third data; calculating a downgoing wavefield based on the difference between the first data and the third data; and separating the upgoing wavefield and the downgoing wavefield based on the upgoing wavefield and the downgoing wavefield.

[0012] In some embodiments, the step of collecting the first data and the second data by detectors in different media based on the seabed node and performing dual-tree complex wavelet transform respectively to obtain the first complex wavelet coefficient and its corresponding first amplitude spectrum and the second complex wavelet coefficient and its corresponding second amplitude spectrum includes: obtaining the first data collected by the underwater detector and the second data collected by the land detector based on the seabed node and performing dual-tree complex wavelet transform respectively to obtain the first complex wavelet coefficient and its corresponding first amplitude spectrum and the second complex wavelet coefficient and its corresponding second amplitude spectrum.

[0013] In another aspect of an embodiment of the present invention, a device for separating upgoing and downgoing wave fields of seabed node data is provided, comprising: a first module for collecting first data and second data by detectors in different media based on the seabed node and performing dual-tree complex wavelet transforms on them respectively to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum and second complex wavelet coefficients and their corresponding second amplitude spectrum, wherein the signal-to-noise ratio of the first data is higher than the signal-to-noise ratio of the second data; a second module for calculating a matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum; a third module for updating the second complex wavelet coefficients based on the matching coefficients and performing a dual-tree complex wavelet inverse transform on the updated second complex wavelet coefficients to obtain third data; and a fourth module for calculating and separating the upgoing and downgoing wave fields based on the first and third data.

[0014] Another aspect of an embodiment of the present invention provides an electronic device comprising at least one processor; and a memory storing computer instructions that can be run on the processor, wherein the instructions implement the steps of the above method when executed by the processor.

[0015] Another aspect of an embodiment of the present invention further provides a computer-readable storage medium storing a computer program that implements the above method steps when executed by a processor.

[0016] The present invention has at least the following beneficial effects: the present invention proposes a method for separating upgoing and downgoing wave fields of seabed node data, performs dual-tree complex wavelet transform on first data and second data collected by detectors in different media based on the seabed node, performs frequency division and multi-directional decomposition, and matches and calibrates the second data with reference to the amplitude spectrum corresponding to the first data with a better signal-to-noise ratio among the first data and the second data, so that the first data and the second data can be reasonably matched, calculates the upgoing wave field and the downgoing wave field based on the second data and the first data after matching and calibration, and separates them, thereby improving the separation effect of the upgoing wave field and the downgoing wave field, and significantly suppressing ghost waves. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other embodiments can be obtained based on these drawings without paying any creative work.

[0018] Figure 1 The figure shows a flow chart of a method for separating uplink and downlink wave fields of seabed node data provided by one embodiment of the present invention;

[0019] FIG2( a ) is a schematic diagram showing a forward transform of a dual-tree complex wavelet transform according to an embodiment of the present invention;

[0020] FIG2( b ) is a schematic diagram showing an inverse transform of a dual-tree complex wavelet transform according to an embodiment of the present invention;

[0021] FIG3( a ) is a schematic diagram showing a pulse a according to an embodiment of the present invention;

[0022] FIG3( b ) is a schematic diagram showing an amplitude spectrum b of a pulse a subjected to a three-level dual-tree complex wavelet transform according to an embodiment of the present invention;

[0023] Figure 4 The figure shows a waveform response diagram of a three-level decomposition of a point pulse dual-tree complex wavelet transform according to an embodiment of the present invention;

[0024] FIG5( a ) is a schematic diagram showing an amplitude spectrum corresponding to water inspection record data in the dual-tree complex wavelet domain according to an embodiment of the present invention;

[0025] FIG5( b ) is a schematic diagram showing the amplitude spectrum corresponding to the land inspection record data in the dual-tree complex wavelet domain provided by an embodiment of the present invention;

[0026] FIG5( c ) is a schematic diagram showing the amplitude spectrum of the land inspection record data in the dual-tree complex wavelet domain after matching and calibration according to an embodiment of the present invention;

[0027] Figure 6 The figure shows a comparison diagram of separating uplink and downlink channel gathers using a dual-tree complex wavelet domain of synthetic data provided by an embodiment of the present invention and a conventional cross-ghost method;

[0028] FIG7( a ) is a schematic diagram showing a superimposed cross-section of water inspection record data provided in one embodiment of the present invention;

[0029] FIG7( b ) is a schematic diagram showing a superimposed cross-section of land inspection record data provided by an embodiment of the present invention;

[0030] FIG7( c ) is a schematic diagram showing a superimposed cross-sectional view of separating upgoing waves using a cross-ghost wave method according to an embodiment of the present invention;

[0031] FIG7( d ) is a schematic cross-sectional view of upgoing wave superposition separated in the dual-tree complex wavelet domain according to an embodiment of the present invention;

[0032] Figure 8 A schematic diagram of a device for separating uplink and downlink wave fields of seabed node data provided by an embodiment of the present invention is shown;

[0033] Figure 9 Shown is a schematic diagram of an electronic device provided by an embodiment of the present invention;

[0034] Figure 10 Shown is a schematic diagram of a computer-readable storage medium provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0035] The following describes embodiments of the present invention. However, it is to be understood that the disclosed embodiments are merely examples and other embodiments can take various alternative forms.

[0036] Furthermore, it should be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but may also include elements not expressly listed or inherent to such processes, methods, articles, or apparatuses.

[0037] One or more embodiments of the present application will be described below with reference to the accompanying drawings.

[0038] Based on the above objectives, a first aspect of an embodiment of the present invention provides an embodiment of a method for uplink and downlink wavefield separation of seabed node data. Figure 1 FIG. 1 is a flow chart of a method for separating uplink and downlink wave fields of seabed node data provided by an embodiment of the present invention. Figure 1As shown, a method for separating uplink and downlink wave fields of seabed node data includes:

[0039] S1. Based on a seabed node, first data and second data are collected by detectors in different media and dual-tree complex wavelet transforms are performed on each of the data to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum, and second complex wavelet coefficients and their corresponding second amplitude spectrum, wherein the signal-to-noise ratio of the first data is higher than the signal-to-noise ratio of the second data;

[0040] S2. Calculating a matching coefficient between the amplitude value of the second amplitude spectrum and the amplitude value of the first amplitude spectrum;

[0041] S3. updating the second complex wavelet coefficient based on the matching coefficient, and performing a dual-tree complex wavelet inverse transform based on the updated second complex wavelet coefficient to obtain third data;

[0042] S4. Calculate an upgoing wavefield and a downgoing wavefield based on the first data and the third data and separate them.

[0043] Among them, the dual-tree complex wavelet transform (DTCWT) is implemented by two sets of parallel real discrete wavelet transforms, as shown in Figure 2(a). Figure 2(a) shows a schematic diagram of the forward transform of the dual-tree complex wavelet transform, and as shown in Figure 2(b), Figure 2(b) shows a schematic diagram of the inverse transform of the dual-tree complex wavelet transform. It can be seen that one set of real discrete wavelet transforms is completed in tree 1, and the other set of real discrete wavelet transforms is completed in tree 2. The subband coefficients generated in tree 1 and tree 2 are the real and imaginary parts of the second complex wavelet coefficients, respectively. Compared with the real discrete wavelet transform, the dual-tree complex wavelet transform has multi-scale frequency division characteristics, better translation invariance and more directional sorting.

[0044] The first data and the second data collected by detectors in different media based on the seabed node are subjected to dual-tree complex wavelet transform, and frequency division and multi-directional decomposition are performed. The amplitude spectrum corresponding to the first data with better signal-to-noise ratio is used as a reference to match and calibrate the second data, so that the first data and the second data can be reasonably matched. The upgoing wave field and the downgoing wave field are calculated based on the second data and the first data after matching and calibration, and are separated, thereby improving the separation effect of the upgoing wave field and the downgoing wave field, and significantly suppressing ghost waves.

[0045] According to several embodiments of the present invention, the step of calculating a matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum includes: calculating the ratio of the amplitude values ​​of the first amplitude spectrum and the second amplitude spectrum at corresponding frequency positions, and obtaining a matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum at the corresponding frequency position based on the amplitude value ratio; and traversing each corresponding frequency position of the first amplitude spectrum and the second amplitude spectrum to obtain each matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum at each corresponding frequency position. Because the signal-to-noise ratio of the first data is better than that of the second data, the amplitude spectrum corresponding to the second data is matched and calibrated using each matching coefficient to match the amplitude spectrum corresponding to the second data with the amplitude spectrum corresponding to the first data with a better signal-to-noise ratio, and accurately calibrate the amplitude and frequency components in the amplitude spectra of the two data, which is conducive to improving the accuracy and precision of subsequent uplink and downlink wavefield separation based on the new second data after matching and calibration.

[0046] According to several embodiments of the present invention, the steps of transforming the complex wavelet corresponding to the second amplitude spectrum based on the matching coefficients and performing an inverse dual-tree complex wavelet transform based on the updated second complex wavelet coefficients to obtain the third data include: updating the second complex wavelet coefficients corresponding to each corresponding frequency position of the second amplitude spectrum based on each matching coefficient; and performing an inverse dual-tree complex wavelet transform based on the updated second complex wavelet coefficients corresponding to each corresponding frequency position to obtain the third data. The third data is obtained by matching and calibrating the second data with the amplitude spectrum of the first data as a reference, and the third data can be reasonably matched with the first data.

[0047] According to some embodiments of the present invention, the step of updating the second complex wavelet coefficient corresponding to the second amplitude spectrum at each corresponding frequency position based on each matching coefficient includes: multiplying the second complex wavelet coefficient corresponding to the second amplitude spectrum at each corresponding frequency position by each matching coefficient corresponding to each corresponding frequency position to obtain an updated second complex wavelet coefficient corresponding to each corresponding frequency position.

[0048] According to some embodiments of the present invention, the step of multiplying the second complex wavelet coefficient corresponding to the second amplitude spectrum at each corresponding frequency position by each matching coefficient corresponding to each corresponding frequency position includes: multiplying the real part and imaginary part of the second complex wavelet coefficient corresponding to the second amplitude spectrum at each corresponding frequency position by each matching coefficient corresponding to each corresponding frequency position, respectively.

[0049] According to several embodiments of the present invention, the steps of calculating and separating the upgoing and downgoing wavefields based on the first and third data include: calculating the upgoing wavefield based on the sum of the first and third data; calculating the downgoing wavefield based on the difference between the first and third data; and separating the upgoing and downgoing wavefields. Because the third data reasonably matches the first data, calculating the upgoing and downgoing wavefields based on the third and first data allows for effective separation of the upgoing and downgoing wavefields, improving the effectiveness of suppressing ghost waves.

[0050] According to several embodiments of the present invention, the steps of collecting first data and second data by detectors in different media based on the seabed node and performing dual-tree complex wavelet transforms respectively to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum and second complex wavelet coefficients and their corresponding second amplitude spectrum include: obtaining first data collected by underwater detectors and second data collected by land detectors based on the seabed node and performing dual-tree complex wavelet transforms respectively to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum and second complex wavelet coefficients and their corresponding second amplitude spectrum.

[0051] The following is another embodiment of a method for uplink and downlink wavefield separation of seafloor node data provided by the present invention. The method involves collecting first and second data corresponding to a submerged geophone and a land geophone, respectively, at a seafloor node. The subsurface effective reflection signal has the same polarity in the data corresponding to the submerged geophone and the land geophone, while the ghost wave has opposite polarity in the data corresponding to the two geophones, resulting in a notch effect. Furthermore, due to the different physical mechanisms of the submerged geophone and the land geophone, the seismic signals collected by them differ in amplitude, frequency, and phase. For example, in terms of amplitude energy, the amplitude value recorded by the land geophone (referred to as the land geophone) is one to two energy levels weaker than that recorded by the submerged geophone (referred to as the submerged geophone). In terms of frequency, the land geophone has a higher sensitivity to seismic signals in the lower frequency band, and the submerged geophone records more high-frequency components. In terms of phase, theoretically, there is a 90° phase difference between the submerged geophone and the land geophone. However, in practice, the phase difference is more complex due to the influence of other factors.

[0052] Based on the aforementioned differences in amplitude, frequency, and phase between water inspection and road inspection data, the uplink and downlink wave fields are separated after effective and reasonable matching and calibration of the two. The specific steps include:

[0053] (1) The water inspection record data is transformed into a dual-tree complex wavelet transform and its amplitude spectrum is calculated.

[0054] (2) Perform dual-tree complex wavelet transform on the land inspection record data and calculate its amplitude spectrum.

[0055] (3) In the dual-tree complex wavelet domain, the amplitude spectrum of the water inspection record data is used as a reference to match and calibrate the amplitude spectrum of the land inspection record data. Specifically, the amplitude spectrum A_P corresponding to the water inspection record data P in the dual-tree complex wavelet domain and the amplitude spectrum A_Z corresponding to the land inspection record data Z in the dual-tree complex wavelet domain are calculated; the amplitude matching calibration coefficient coe is calculated to match the amplitude values ​​of the two at all corresponding positions of the amplitude spectrum, wherein the specific formula of the amplitude matching calibration coefficient is coe=A_P / A_Z; the amplitude matching calibration coefficient is applied to the real part R_Z and the imaginary part I_Z of the complex wavelet coefficient of the land inspection record data at the corresponding position, respectively, to obtain the real part R_Z_M=R_Z*coe and the imaginary part I_Z_M=I_Z*coe of the complex wavelet coefficient after amplitude matching calibration.

[0056] (4) Performing a dual-tree complex wavelet inverse transform on the complex wavelet coefficients of the land inspection record after amplitude matching calibration in the dual-tree complex wavelet domain to obtain the land inspection record data Z_M after matching calibration.

[0057] (5) The water inspection record data P and the land inspection record data Z_M after matching and calibration are added and subtracted respectively to obtain the upgoing wave field and the downgoing wave field. The specific calculation formulas for the upgoing wave field and the downgoing wave field are: upgoing wave field U1 = 0.5*(P+Z_M), downgoing wave field U2 = 0.5*(P-Z_M).

[0058] The signal is frequency-divided and multi-directionally decomposed through dual-tree complex wavelet transform, and the amplitude and frequency components of the water inspection record data and land inspection record data collected based on the seabed node are accurately calibrated in the dual-tree complex wavelet domain, which improves the separation effect of the upgoing wave field and the downgoing wave field and can effectively suppress ghost waves.

[0059] The following is another embodiment of a method for separating uplink and downlink wave fields of seabed node data provided by the present invention, and the specific steps include:

[0060] (1) Perform a dual-tree complex wavelet transform on the water inspection record data and calculate its amplitude spectrum, where the subband coefficients generated in Tree 1 and Tree 2 are the real and imaginary parts of the complex wavelet coefficients, respectively. Figure 3(a) shows a schematic diagram of pulse a provided by one embodiment of the present invention. As shown in Figure 3(a), pulse a has 128 rows and 128 columns. Figure 3(b) is a schematic diagram of the amplitude spectrum b of the three-level dual-tree complex wavelet transform of pulse a provided by an embodiment of the present invention. As shown in Figure 3(b), the number of rows of the amplitude spectrum after the three-level decomposition of the dual-tree complex wavelet transform of pulse a is nr=2*nr0=258, and the number of columns is nc0=2*nc0=258, wherein the quadrants q11, q12, and q13 are the first-level decomposition complex wavelet coefficient amplitude spectrum, the quadrants q21, q22, and q23 are the second-level decomposition complex wavelet coefficient amplitude spectrum, and the quadrants q31, q32, q33, and q34 are the last level, i.e., the third-level decomposition complex wavelet coefficient amplitude spectrum.

[0061] Figure 4 FIG. 1 is a schematic diagram of a waveform response of a point pulse dual-tree complex wavelet transform three-level decomposition according to an embodiment of the present invention. Figure 4 As shown, in order from top to bottom, the first row represents the waveform response corresponding to the real and imaginary parts of the first-level complex wavelet coefficients q11, q12, and q13 quadrants, the second row represents the waveform response corresponding to the real and imaginary parts of the second-level complex wavelet coefficients q21, q22, and q23 quadrants, and the third row represents the waveform response corresponding to the real and imaginary parts of the third-level complex wavelet coefficients q31, q32, and q33 quadrants.

[0062] It can be seen from the waveform responses corresponding to each row above that the dual-tree complex wavelet transform has a good effect on the frequency and direction decomposition of the data. Based on the dual-tree complex wavelet transform, the amplitude and frequency components of the water inspection record data and the land inspection record data can be accurately calibrated, which is conducive to the effective separation of the uplink and downlink wave fields.

[0063] (2) Perform dual-tree complex wavelet transform on the land inspection record data and calculate its amplitude spectrum.

[0064] (3) FIG5(a) shows a schematic diagram of the amplitude spectrum of the water inspection record data in the dual-tree complex wavelet domain provided by an embodiment of the present invention, and FIG5(b) shows a schematic diagram of the amplitude spectrum of the land inspection record data in the dual-tree complex wavelet domain provided by an embodiment of the present invention. In the dual-tree complex wavelet domain, the amplitude spectrum of the land inspection record data is matched and calibrated with reference to the amplitude spectrum of the water inspection record data. Specifically, the amplitude spectrum A_P corresponding to the water inspection record data P in the dual-tree complex wavelet domain and the amplitude spectrum A_Z corresponding to the land inspection record data Z in the dual-tree complex wavelet domain are calculated; the amplitudes of the two are calculated at all corresponding positions of the amplitude spectrum. The amplitude matching calibration coefficient coe is matched with the amplitude matching value, wherein the specific formula of the amplitude matching calibration coefficient is coe=A_P / A_Z; the amplitude matching calibration coefficient is applied to the real part R_Z and the imaginary part I_Z of the complex wavelet coefficient of the land inspection record data at the corresponding position, respectively, to obtain the real part R_Z_M=R_Z*coe and the imaginary part I_Z_M=I_Z*coe of the complex wavelet coefficient after amplitude matching calibration. FIG5(c) shows a schematic diagram of the amplitude spectrum corresponding to the land inspection record data after matching calibration in the dual-tree complex wavelet domain provided by one embodiment of the present invention.

[0065] (4) Performing a dual-tree complex wavelet inverse transform on the complex wavelet coefficients of the land inspection record after amplitude matching calibration in the dual-tree complex wavelet domain to obtain the land inspection record data Z_M after matching calibration.

[0066] (5) The water inspection record data P and the land inspection record data Z_M after matching and calibration are added and subtracted respectively to obtain the upgoing wave field and the downgoing wave field. The specific calculation formulas for the upgoing wave field and the downgoing wave field are: upgoing wave field U1 = 0.5*(P+Z_M), downgoing wave field U2 = 0.5*(P-Z_M). Figure 6 The figure shows a schematic diagram comparing the separation of upgoing and downgoing wave gathers using the dual-tree complex wavelet domain of synthetic data provided by an embodiment of the present invention and the conventional cross-ghost method. The two methods are compared from the perspectives of water inspection recording water channels, road inspection recording water channels, and separation of upgoing and downgoing wave fields. Figure 6 (a) in the figure is the water inspection record gather, Figure 6 (b) in the figure is the Lujian record gather. Figure 6 (c) is the conventional cross-ghost wave method to separate the upgoing wave field. Figure 6 (d) in the figure is the conventional cross-ghost wave method to separate the downlink wave field. Figure 6 (e) in the figure is the upgoing wave field separated by the dual-tree complex wavelet domain. Figure 6 Figure (f) shows the downgoing wavefield separated by the dual-tree complex wavelet domain. The above comparison shows that compared to the original water and land inspection data, the ghost waves in the upgoing wavefield are significantly suppressed, and the valid next reflection information in the downgoing wavefield is also effectively suppressed. Compared with the conventional cross-ghost method, the dual-tree complex wavelet domain separation of the upgoing and downgoing wavefields is more effective in suppressing ghost waves.

[0067] To further compare the effectiveness of separating upgoing and downgoing wavefields using the dual-tree complex wavelet domain and the conventional cross-ghost method, stacked profiles corresponding to separation of upgoing waves using the dual-tree complex wavelet domain and the conventional cross-ghost method were obtained and compared. Specifically, Figure 7(a) shows a schematic diagram of a stacked profile of water inspection record data according to an embodiment of the present invention, Figure 7(b) shows a schematic diagram of a stacked profile of land inspection record data according to an embodiment of the present invention, Figure 7(c) shows a schematic diagram of a stacked profile of upgoing waves separated using the cross-ghost method according to an embodiment of the present invention, and Figure 7(d) shows a schematic diagram of a stacked profile of upgoing waves separated using the dual-tree complex wavelet domain according to an embodiment of the present invention. As can be seen from the above, compared to the stacked profile corresponding to separation of upgoing and downgoing waves using the cross-ghost method, the effective reflection information in the stacked profile of upgoing waves separated using the dual-tree complex wavelet domain is significantly enhanced. Furthermore, the effective reflection events in the stacked profile of upgoing waves separated using the dual-tree complex wavelet domain are more continuous, resulting in a higher signal-to-noise ratio.

[0068] The second aspect of the embodiment of the present invention provides a device for separating uplink and downlink wave fields of seabed node data. Figure 8 FIG. 1 shows a schematic diagram of a device for separating uplink and downlink wave fields of seabed node data provided by an embodiment of the present invention. Figure 8 As shown, it includes: a first module 011, which is used to collect first data and second data by detectors in different media based on the seabed node and perform dual-tree complex wavelet transform on them respectively to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum and second complex wavelet coefficients and their corresponding second amplitude spectrum, wherein the signal-to-noise ratio of the first data is higher than the signal-to-noise ratio of the second data; a second module 012, which is used to calculate a matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum; a third module 013, which is used to update the second complex wavelet coefficients based on the matching coefficients, and perform dual-tree complex wavelet inverse transform on the updated second complex wavelet coefficients to obtain third data; a fourth module 014, which is used to calculate and separate the upgoing wave field and the downgoing wave field based on the first data and the third data.

[0069] According to a third aspect of the embodiments of the present invention, an electronic device is provided. Figure 9 FIG. 1 is a schematic diagram of an electronic device provided by an embodiment of the present invention. Figure 9 As shown, an electronic device provided by an embodiment of the present invention includes the following modules: at least one processor 021; and a memory 022, wherein the memory 022 stores computer instructions 023 that can be run on the processor 021, and the computer instructions 023 implement the steps of the method described above when executed by the processor 021.

[0070] The present invention also provides a computer-readable storage medium. Figure 10 FIG. 1 is a schematic diagram showing the structure of a computer-readable storage medium provided by an embodiment of the present invention. Figure 10 As shown, the computer readable storage medium 031 stores a computer program 032 for executing the steps of the method described above when executed by a processor. The method executed is the same as above.

[0071] Finally, it should be noted that those skilled in the art will understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program. The program for setting the system parameters can be stored in a computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above-mentioned methods. The storage medium of the program can be a magnetic disk, an optical disk, a read-only memory (ROM), or a random access memory (RAM). The above-mentioned computer program embodiments can achieve the same or similar effects as any of the corresponding aforementioned method embodiments.

[0072] In addition, the method disclosed in the embodiment of the present invention can also be implemented as a computer program executed by a processor, and the computer program can be stored in a computer-readable storage medium. When the computer program is executed by the processor, the above functions defined in the method disclosed in the embodiment of the present invention are performed.

[0073] In addition, the above method steps and system units can also be implemented using a controller and a computer-readable storage medium for storing a computer program that enables the controller to implement the above steps or unit functions.

[0074] It will also be appreciated by those skilled in the art that the various exemplary logic blocks, modules, circuits and algorithmic steps described in conjunction with the disclosure herein can be implemented as electronic hardware, computer software or a combination of the two. In order to clearly illustrate this interchangeability of hardware and software, a general description has been given of the functions of various schematic components, blocks, modules, circuits and steps. Whether this function is implemented as software or hardware depends on specific applications and the design constraints imposed on the entire system. Those skilled in the art can implement the function in various ways for each specific application, but this implementation decision should not be interpreted as causing a departure from the disclosed scope of the embodiments of the present invention.

[0075] In one or more exemplary designs, the function can be implemented in hardware, software, firmware or any combination thereof. If implemented in software, the function can be stored as one or more instructions or codes on a computer-readable medium or transmitted via a computer-readable medium. Computer-readable media include computer storage media and communication media, and the communication media include any media that helps to transmit a computer program from one location to another. The storage medium can be any available medium that can be accessed by a general or special-purpose computer. As an example and not limitation, the computer-readable medium can include RAM, ROM, EEPROM, CD-ROM or other optical disk storage devices, magnetic disk storage devices or other magnetic storage devices, or any other medium that can be used to carry or store the required program code in the form of an instruction or data structure and can be accessed by a general or special-purpose computer or a general or special-purpose processor. In addition, any connection can be appropriately referred to as a computer-readable medium. For example, if a coaxial cable, fiber optic cable, twisted pair, digital subscriber line (DSL) or wireless technologies such as infrared, radio and microwaves are used to send software from a website, server or other remote source, the above-mentioned coaxial cable, fiber optic cable, twisted pair, DOL or wireless technologies such as infrared, radio and microwaves are all included in the definition of medium. As used herein, disk and disc include compact disc (CD), laser disc, optical disc, digital versatile disc (DVD), floppy disk, and Blu-ray disc, where disks usually reproduce data magnetically, while discs reproduce data optically with lasers. Combinations of the above should also be included within the scope of computer-readable media.

[0076] The above are exemplary embodiments disclosed in the present invention, but it should be noted that various changes and modifications may be made without departing from the scope of the embodiments disclosed in the claims. The functions, steps and / or actions of the method claims according to the disclosed embodiments described herein do not need to be performed in any particular order. In addition, although the elements disclosed in the embodiments of the present invention may be described or required in individual form, they may also be understood as multiple unless expressly limited to the singular.

[0077] It should be understood that, as used herein, the singular forms "a" and "an" are intended to include the plural forms as well, unless the context clearly supports an exception. It should also be understood that, as used herein, "and / or" is intended to include any and all possible combinations of one or more of the associated listed items.

[0078] The serial numbers of the embodiments disclosed in the above embodiments of the present invention are only for description and do not represent the advantages or disadvantages of the embodiments.

[0079] Those skilled in the art will understand that all or part of the steps to implement the above embodiments may be accomplished by hardware, or may be accomplished by a program instructing the relevant hardware, and the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a disk, or an optical disk, etc.

[0080] Those skilled in the art should understand that the discussion of any of the above embodiments is merely illustrative and is not intended to imply that the scope of the disclosure of the embodiments of the present invention (including the claims) is limited to these examples. Within the spirit of the embodiments of the present invention, the technical features of the above embodiments or different embodiments may be combined, and there are many other variations of different aspects of the above embodiments of the present invention, which are not provided in detail for the sake of simplicity. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of the present invention should be included in the scope of protection of the embodiments of the present invention.

Claims

1. A method for separating uplink and downlink wave fields of seabed node data, characterized in that: include: Based on a seabed node, first data and second data are collected by detectors in different media and dual-tree complex wavelet transforms are performed on each of the data to obtain first complex wavelet coefficients and a corresponding first amplitude spectrum and second complex wavelet coefficients and a corresponding second amplitude spectrum, wherein a signal-to-noise ratio of the first data is higher than a signal-to-noise ratio of the second data; Calculating a matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum, comprising: calculating a ratio of the amplitude values ​​of the first amplitude spectrum and the second amplitude spectrum at corresponding frequency positions, and obtaining a matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum at the corresponding frequency position based on the amplitude value ratio; traversing each corresponding frequency position of the first amplitude spectrum and the second amplitude spectrum, and obtaining each matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum at each corresponding frequency position; Updating the second complex wavelet coefficient based on the matching coefficient, and performing a dual-tree complex wavelet inverse transform based on the updated second complex wavelet coefficient to obtain third data, which includes: updating the second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum based on each matching coefficient; performing a dual-tree complex wavelet inverse transform based on the updated second complex wavelet coefficient corresponding to each corresponding frequency position to obtain third data; An upgoing wavefield and a downgoing wavefield are calculated based on the first data and the third data and separated.

2. The method according to claim 1, characterized in that The step of updating the second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum based on each matching coefficient includes: The second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum is multiplied by each matching coefficient corresponding to each corresponding frequency position to obtain an updated second complex wavelet coefficient corresponding to each corresponding frequency position.

3. The method according to claim 2, characterized in that The step of multiplying the second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum by each matching coefficient corresponding to each corresponding frequency position includes: The real part and the imaginary part of the second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum are respectively multiplied by each matching coefficient corresponding to each corresponding frequency position.

4. The method according to claim 1, wherein The step of calculating and separating the upgoing wavefield and the downgoing wavefield based on the first data and the third data includes: Calculating an upgoing wavefield based on the sum of the first data and the third data; calculating a downlink wavefield based on a difference between the first data and the third data; The upgoing wavefield and the downgoing wavefield are separated based on the upgoing wavefield and the downgoing wavefield.

5. The method according to claim 1, wherein The steps of collecting first data and second data by detectors in different media based on the seabed node and performing dual-tree complex wavelet transform respectively to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum and second complex wavelet coefficients and their corresponding second amplitude spectrum include: Based on the seabed node, first data collected by the underwater detector and second data collected by the land detector are obtained and double-tree complex wavelet transform is performed on them respectively to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum and second complex wavelet coefficients and their corresponding second amplitude spectrum.

6. A device for separating uplink and downlink wave fields of seabed node data, characterized in that: include: The first module is configured to collect first data and second data using detectors in different media based on a seabed node and perform dual-tree complex wavelet transforms on each of the data to obtain first complex wavelet coefficients and their corresponding first amplitude spectrum, and second complex wavelet coefficients and their corresponding second amplitude spectrum, wherein the signal-to-noise ratio of the first data is higher than the signal-to-noise ratio of the second data; A second module is used to calculate a matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum; a third module, configured to update the second complex wavelet coefficient based on the matching coefficient, and perform a dual-tree complex wavelet inverse transform based on the updated second complex wavelet coefficient to obtain third data; A fourth module is configured to calculate and separate an upgoing wavefield and a downgoing wavefield based on the first data and the third data; The second module is further configured to calculate an amplitude ratio between the first amplitude spectrum and the second amplitude spectrum at corresponding frequency positions, and obtain a matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum at the corresponding frequency position based on the amplitude ratio; and traverse each corresponding frequency position of the first amplitude spectrum and the second amplitude spectrum to obtain each matching coefficient for matching the amplitude value of the second amplitude spectrum with the amplitude value of the first amplitude spectrum at each corresponding frequency position; The fourth module is also used to update the second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum based on each matching coefficient; and perform a dual-tree complex wavelet inverse transform based on the updated second complex wavelet coefficient corresponding to each corresponding frequency position to obtain third data.

7. An electronic device, characterized in that: include: at least one processor; as well as A memory storing computer instructions executable on the processor, wherein the computer instructions, when executed by the processor, implement the steps of the method according to any one of claims 1 to 5.

8. A computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the method according to any one of claims 1 to 5.

Images