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

Abstract

The invention relates to the field of geophysical signal processing, and provides a method, a device, equipment and a medium for separating uplink and downlink wave fields of seabed node data, and the method comprises the steps: collecting first data and second data through detectors in different mediums based on seabed nodes, and carrying out the dual-tree complex wavelet transformation, a first complex wavelet coefficient and a first amplitude spectrum corresponding to the first complex wavelet coefficient as well as a second complex wavelet coefficient and a second amplitude spectrum corresponding to the second complex wavelet coefficient are obtained, and the signal-to-noise ratio of the first data is higher than that of the second data; calculating a matching coefficient between the amplitude value of the second amplitude spectrum and the amplitude value of the first amplitude spectrum; updating the second complex wavelet coefficient based on the matching coefficient, and performing dual-tree complex wavelet inverse transformation based on the updated second complex wavelet coefficient to obtain third data; and calculating an uplink wave field and a downlink wave field based on the first data and the third data and separating the uplink wave field and the downlink wave field. According to the scheme disclosed by the invention, the separation effect of the uplink wave field and the downlink wave field is improved, and the ghost wave suppression effect is obvious.

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 Nodes (OBN) seismic exploration, ghost waves are generated in the data received from the seabed because the sea surface is a strong reflection interface. Ghost wave interference will form false phase axes in seismic records and cause notch effects in the spectrum, thereby reducing the resolution of effective underground reflection seismic signals and affecting the accuracy of seismic data interpretation.

[0003] For the seismic data collected by OBN and received by underwater geophones and land geophones, the underground effective reflection signals have the same polarity in the seismic data received by the underwater geophones and land geophones, while the ghost waves have the opposite polarity, and the ghost waves will cause complementary notches. Therefore, the seismic data received by the underwater geophones and land geophones are generally processed to obtain the upgoing wave field and the downgoing wave field, and the ghost waves are suppressed by separating the upgoing and downgoing wave fields.

[0004] Among them, the upgoing wave field and the downgoing wave field are calculated after matching and calibrating the seismic data received by the underwater geophone and the land geophone. The matching calibration factor is obtained in 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-variable; the matching operator is obtained after cross-ghosting, and the ghost wave extension operator is derived based on the assumption of seabed horizontality, which is not applicable to complex seabed conditions. Based on the above problems existing in the prior art, 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 separating the uplink and downlink wave fields of seabed node data, which at least solves the problem that in the related technologies of uplink and downlink wave field separation, the seismic data received by the underwater detector and the land detector 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 purpose, an aspect of an embodiment of the present invention provides a method for separating upgoing and downgoing wave fields of seabed node data, including: based on the seabed node, first data and second data are collected by detectors in different media and dual-tree complex wavelet transforms are performed 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 matching coefficient is calculated so that the amplitude value of the second amplitude spectrum matches the amplitude value of the first amplitude spectrum; the second complex wavelet coefficient is updated based on the matching coefficient, and a dual-tree complex wavelet inverse transform is performed based on the updated second complex wavelet coefficient to obtain third data; and the upgoing and downgoing wavefields are calculated based on the first data and the third data and separated.

[0007] In some embodiments, the step of calculating the matching coefficient of the amplitude value of the second amplitude spectrum to match 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 of the amplitude value of the second amplitude spectrum to match 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 of the amplitude value of the second amplitude spectrum to match 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 an 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 the second amplitude spectrum at each corresponding frequency position by each matching coefficient corresponding to the each corresponding frequency position includes: multiplying the real part and the 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 the 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 comprises: calculating an upgoing wavefield based on the sum of the first data and the third data; calculating a downgoing wavefield based on a 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: based on the seabed node, obtaining the first data collected by the underwater detector and the second data collected by the land detector 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 uplink and downlink 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 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, 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 coefficient based on the matching coefficient, and performing dual-tree complex wavelet inverse transform based on the updated second complex wavelet coefficient to obtain third data; and a fourth module, for calculating and separating the uplink and downlink wave fields based on the first data and the third data.

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

[0015] According to another aspect of an embodiment of the present invention, a computer-readable storage medium is provided, which stores 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, matches and calibrates the second data with reference to the amplitude spectrum corresponding to the first data with 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 obviously 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 drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings in the following description 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 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] FIG. 2( a ) is a schematic diagram showing a forward transformation of a dual-tree complex wavelet transform according to an embodiment of the present invention;

[0020] FIG. 2( b ) is a schematic diagram showing an inverse transformation 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 provided in one 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 schematic diagram of a waveform response of a three-level decomposition of a point pulse dual-tree complex wavelet transform provided by an embodiment of the present invention;

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

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

[0026] FIG5(c) is a schematic diagram showing an amplitude spectrum corresponding to the land inspection record data in the dual-tree complex wavelet domain after matching and calibration provided in 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 wave method;

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

[0029] FIG. 7( b ) is a schematic diagram showing a superimposed cross-section of land inspection record data provided in 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 diagram showing a cross-section of upgoing wave superposition separated in a 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] Fig. 9 Shown is a schematic diagram of an electronic device provided by an embodiment of the present invention;

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

[0035] The following describes embodiments of the present invention. However, it should 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 "includes," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that includes a series of elements includes not only those elements, but may also include elements not explicitly listed or inherent to these 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 separating uplink and downlink wave fields 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 the seabed node, the first data and the second data are collected by the detectors in different media and double-tree complex wavelet transform is performed respectively to obtain the first complex wavelet coefficient and the first amplitude spectrum corresponding thereto and the second complex wavelet coefficient and the second amplitude spectrum corresponding thereto, 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 for matching the amplitude value of the second amplitude spectrum with 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 realized 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, 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 a multi-scale frequency division characteristic, and also has 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 second data is matched and calibrated with reference to the amplitude spectrum corresponding to the first data with 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. The upgoing wave field and the downgoing wave field are calculated and separated based on the second data and the first data after matching and calibration, thereby improving the separation effect of the upgoing wave field and the downgoing wave field, and obviously suppressing the ghost wave.

[0045] According to some embodiments of the present invention, the step of calculating the matching coefficient of the amplitude value of the second amplitude spectrum to match the amplitude value of the first amplitude spectrum includes: calculating the amplitude value ratio of the first amplitude spectrum to the second amplitude spectrum at the corresponding frequency position, so as to obtain the matching coefficient of the second amplitude spectrum to match 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 of the second amplitude spectrum to match the amplitude value of the first amplitude spectrum at each corresponding frequency position. Since 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 by each matching coefficient, so as 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, which is conducive to improving the accuracy and precision of the subsequent uplink and downlink wave field separation based on the new second data after matching and calibration.

[0046] According to some embodiments of the present invention, the steps of transforming the complex wavelet corresponding to the second amplitude spectrum based on the matching coefficient and performing a dual-tree complex wavelet inverse 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; performing a dual-tree complex wavelet inverse 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 the 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 some embodiments of the present invention, the steps of calculating and separating the upgoing wavefield and the downgoing wavefield based on the first data and the third data include: calculating the upgoing wavefield based on the sum of the first data and the third data; calculating the downgoing wavefield based on the difference between the first data and the third data; and separating the upgoing wavefield and the downgoing wavefield. Since the third data is data that can reasonably match the first data, the upgoing wavefield and the downgoing wavefield are calculated based on the third data and the first data, so that the upgoing wavefield and the downgoing wavefield can be effectively separated, thereby improving the suppression effect on ghost waves.

[0050] According to several embodiments of the present invention, based on the seabed node, the first data and the second data are collected by detectors in different media and dual-tree complex wavelet transforms are performed 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. The step includes: based on the seabed node, the first data collected by the underwater detector and the second data collected by the land detector are obtained and dual-tree complex wavelet transforms are performed 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.

[0051] The following is another embodiment of a method for separating uplink and downlink wave fields of seafloor node data provided by the present invention, based on collecting first data and second data corresponding to the underwater geophone and the land geophone respectively at the seafloor node, wherein the underground effective reflection signal has the same polarity in the data corresponding to the underwater geophone and the land geophone respectively, while the ghost wave has opposite polarity in the data corresponding to the two respectively, and a notch effect will be generated. At the same time, due to the different physical mechanisms of the underwater geophone and the land geophone, the seismic signals collected by them have differences in amplitude, frequency and phase, for example, in terms of amplitude energy, the amplitude value recorded by the land geophone (referred to as land detection) is weaker than the amplitude value recorded by the underwater geophone (referred to as water detection) by more than 1 to 2 energy levels; in terms of frequency, the land geophone has a higher sensitivity to seismic signals in the lower frequency band range, and the seismic signal recorded by the underwater geophone has more high-frequency components; in terms of phase, theoretically, there is a 90° phase difference between the underwater geophone and the land geophone, but in reality, affected by other factors, the phase difference between the two is more complicated.

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

[0053] (1) Perform dual-tree complex wavelet transform on the water inspection record data and calculate its amplitude spectrum.

[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 and calibration in the dual-tree complex wavelet domain to obtain the land inspection record data Z_M after matching and 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 calculating the upgoing wave field and the downgoing wave field are as follows: the upgoing wave field U1 = 0.5*(P+Z_M), and the downgoing wave field U2 = 0.5*(P-Z_M).

[0058] The signal is divided into frequencies and decomposed in multiple directions through dual-tree complex wavelet transform. The amplitude and frequency components of the water inspection record data and land inspection record data collected by the seabed nodes 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 uplink and downlink wave field separation 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 part and imaginary part of the complex wavelet coefficients, respectively. FIG3(a) shows a schematic diagram of a pulse a provided by an embodiment of the present invention. As shown in FIG3(a), the number of data rows of pulse a is nr0=128, and the number of columns is nc0=128. Figure 3(b) is a schematic diagram of the amplitude spectrum b of the three-level dual-tree complex wavelet transform of the pulse a provided by an embodiment of the present invention. As shown in Figure 3(b), the number of rows of the amplitude spectrum of pulse a after the three-level decomposition of the dual-tree complex wavelet transform 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 spectra, the quadrants q21, q22, and q23 are the second-level decomposition complex wavelet coefficient amplitude spectra, 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, 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, q23 quadrants, and the third row represents the waveform response corresponding to the real and imaginary parts of the tertiary complex wavelet coefficients q31, q32, 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 beneficial 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) FIG. 5(a) is a schematic diagram showing an amplitude spectrum of water inspection record data in the dual-tree complex wavelet domain provided by an embodiment of the present invention, and FIG. 5(b) is a schematic diagram showing an amplitude spectrum of 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 that matches the value is obtained, 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 in an 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 and calibration in the dual-tree complex wavelet domain to obtain the land inspection record data Z_M after matching and 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 calculating the upgoing wave field and the downgoing wave field are as follows: the upgoing wave field U1 = 0.5*(P+Z_M), and the downgoing wave field U2 = 0.5*(P-Z_M). Figure 6 The figure shows a schematic diagram of the comparison between the dual-tree complex wavelet domain of synthetic data provided by an embodiment of the present invention and the conventional cross ghost wave method for separating upgoing and downgoing wave channels. The two are compared from the perspective of the effects of recording water channels for water inspection and road inspection, and separating upgoing and downgoing wave fields. Figure 6 (a) 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 (f) in the figure is the dual-tree complex wavelet domain separation of the downgoing wave field. From the above comparison, it can be seen that compared with the original water inspection record data and land inspection record data, the ghost waves in the upgoing wave field are significantly suppressed, and the effective next reflection information in the downgoing wave field is also effectively suppressed; compared with the conventional cross ghost wave method, the dual-tree complex wavelet domain separation of the upgoing and downgoing wave fields has a more significant effect in suppressing ghost waves.

[0067] In order to further compare the effects of separating the upgoing and downgoing wave fields in the dual-tree complex wavelet domain and the conventional cross ghost wave method, the stacked sections corresponding to the separation of the upgoing waves based on the dual-tree complex wavelet domain and the conventional cross ghost wave method are obtained and compared. Specifically, FIG7(a) shows a schematic diagram of the stacked section of the water inspection record data provided by an embodiment of the present invention, FIG7(b) shows a schematic diagram of the stacked section of the land inspection record data provided by an embodiment of the present invention, FIG7(c) shows a schematic diagram of the stacked section of the upgoing waves separated by the cross ghost wave method provided by an embodiment of the present invention, and FIG7(d) shows a schematic diagram of the stacked section of the upgoing waves separated by the dual-tree complex wavelet domain provided by an embodiment of the present invention. It can be seen from the above that compared with the stacked section corresponding to the separation of the upgoing and downgoing waves by the cross ghost wave method, the effective reflection information in the stacked section of the upgoing waves separated by the dual-tree complex wavelet domain is significantly enhanced. At the same time, the effective reflection event axis in the stacked section of the upgoing waves separated by the dual-tree complex wavelet domain is more continuous, and the data signal-to-noise ratio is higher.

[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 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. 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 a seabed node and perform dual-tree complex wavelet transform respectively to obtain a first complex wavelet coefficient and a first amplitude spectrum corresponding to it and a second complex wavelet coefficient and a second amplitude spectrum corresponding to it, 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 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 014, which is used to calculate and separate an upgoing wave field and a downgoing wave field based on the first data and the third data.

[0069] According to a third aspect of an embodiment of the present invention, an electronic device is provided. Fig. 9 FIG. 1 is a schematic diagram of an electronic device provided by an embodiment of the present invention. Fig. 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 executed 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. Fig.10 FIG. 1 is a schematic diagram showing the structure of a computer-readable storage medium provided by an embodiment of the present invention. Fig.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 a person of ordinary skill in the art can 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, and the program of the method for setting system parameters can be stored in a computer-readable storage medium, and when the program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, the storage medium of the program can be a disk, an optical disk, a read-only storage memory (ROM) or a random access memory (RAM), etc. The above-mentioned computer program embodiments can achieve the same or similar effects as the corresponding above-mentioned arbitrary 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 may 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 various exemplary logic blocks, modules, circuits and algorithm 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 to 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 on a computer-readable medium or transmitted by a computer-readable medium as one or more instructions or codes. Computer-readable media include computer storage media and communication media, and the communication media include any media that helps to transfer 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 computer. As an example and not limiting, the computer-readable medium may 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 computer or a general or special processor. In addition, any connection can be appropriately referred to as a computer-readable medium. For example, if a coaxial cable, optical fiber 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, optical fiber cable, twisted pair, DOL or wireless technologies such as infrared, radio and microwaves are all included in the definition of the medium. As used herein, disk and disc include compact disc (CD), laser disc, optical disc, digital versatile disc (DVD), floppy disk, Blu-ray disc, wherein disks usually reproduce data magnetically, while optical 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 disclosed in the embodiments of the present invention as defined 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 explicitly limited to the singular.

[0077] It should be understood that, as used herein, the singular forms "a", "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" refers to any and all possible combinations including 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] A person skilled in the art will understand that all or part of the steps to implement the above embodiments may be accomplished by hardware or by instructing related hardware through a program, 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] A person skilled in the art should understand that the discussion of any of the above embodiments is only exemplary 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; under the concept of the embodiments of the present invention, the technical features in the above embodiments or different embodiments can also be combined, and there are many other changes in different aspects of the embodiments of the present invention as above, 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 protection scope 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 the seabed node, first data and second data are collected by detectors in different media and double-tree complex wavelet transforms are performed 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 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; An upgoing wavefield and a downgoing wavefield are calculated and separated based on the first data and the third data.

2. The method according to claim 1, characterized in that 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 comprises: Calculating an amplitude value ratio of the first amplitude spectrum to the second amplitude spectrum at a corresponding frequency position, so as to obtain a matching coefficient for matching the amplitude value of the second amplitude spectrum to the amplitude value of the first amplitude spectrum at the corresponding frequency position based on the amplitude value ratio; Each corresponding frequency position of the first amplitude spectrum and the second amplitude spectrum is traversed to obtain each matching coefficient of the amplitude value of the second amplitude spectrum matching the first amplitude spectrum at each corresponding frequency position.

3. The method according to claim 2, characterized in that 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 comprises: Based on each matching coefficient, updating the second complex wavelet coefficient corresponding to each corresponding frequency position of the second amplitude spectrum; A dual-tree complex wavelet inverse transform is performed based on the updated second complex wavelet coefficient corresponding to each corresponding frequency position to obtain third data.

4. The method according to claim 3, 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 comprises: 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.

5. The method according to claim 4, 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 comprises: 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.

6. The method according to claim 1, characterized in that The step of calculating and separating an upgoing wavefield and a downgoing wavefield based on the first data and the third data comprises: Calculate an upgoing wavefield based on the sum of the first data and the third data; Calculating a downlink wave field 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.

7. The method according to claim 1, characterized in that 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 the first amplitude spectrum corresponding thereto and the second complex wavelet coefficient and the second amplitude spectrum corresponding thereto comprises: Based on the seabed node, the first data collected by the underwater detector and the second data collected by the land detector are obtained and double-tree complex wavelet transform is performed 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.

8. A device for separating uplink and downlink wave fields of seabed node data, characterized in that: include: The first module 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 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 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 is used 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; The fourth module is used to calculate and separate the upgoing wavefield and the downgoing wavefield based on the first data and the third data.

9. 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 instructions, when executed by the processor, implement the steps of the method according to any one of claims 1 to 7.

10. A computer-readable storage medium storing a computer program for implementing the steps of the method according to any one of claims 1 to 7 when the computer-readable storage medium is executed by a processor.

Images