A method and device for calculating a double-check matching operator of seismic data of a seabed node

Abstract

The application provides a kind of seabed node data double-detection matching operator calculation method and device, the method includes based on the water depth of excitation point, the propagation speed of seismic wave in water, the propagation speed of seismic wave in seabed and the offset of first arrival direct wave, the range of direct wave is calculated;Linear moveout correction and common receiver point stack are carried out to the hydrophone direct wave and the seismic detector direct wave in the range of direct wave respectively, to obtain hydrophone common receiver point stack data and seismic detector common receiver point stack data;Based on hydrophone common receiver point stack data and seismic detector common receiver point stack data, double-detection matching operator of each receiver point is calculated in the preset time window by least square method.The seabed node data double-detection matching operator calculation method provided by the application can use direct wave to obtain stable and reliable double-detection matching operator, and provide guarantee for subsequent data processing.

Description

Technical Field

[0001] This application relates to the field of oil and gas exploration, and in particular to a method and apparatus for calculating a dual-detection matching operator for seismic data of seafloor nodes. Background Technology

[0002] In the field of subsea oil and gas exploration, marine seismic data acquisition is indispensable. However, traditional methods such as towed cables or ocean bottom cables (OBCs) are difficult to implement when faced with complex operational obstacles, and their data acquisition precision is insufficient for achieving the desired subsurface imaging. Therefore, ocean bottom node (OBN) marine seismic data acquisition has become the primary method for marine oil and gas exploration in complex seabed conditions. OBN marine seismic data acquisition can achieve better azimuth and offset distribution, higher coverage times, and longer offsets. Furthermore, OBN acquisition systems can provide more subsurface information and obtain better reservoir imaging when acquiring four-component data.

[0003] The main challenges in seafloor node seismic exploration lie in data processing and imaging. Due to the ghost wave problem, up and down wavefield separation occurs during the detection process. This separation is achieved by merging hydrophone and seismic detector data in the vertical component of the seismic data, separating the wavefield into up and down waves. To address this issue, direct vertical stacking or using fixed matching coefficients for correction before merging hydrophone and seismic detector data can be used to suppress ghost waves and achieve wavefield separation. However, these methods are highly inaccurate and cannot adapt to variations in water depth and seabed medium, resulting in poor wavefield separation performance.

[0004] Current methods for summing and separating uplink and downlink wavefields using dual-detector systems include spatiotemporally variable matching algorithms and classical algorithms. Matching algorithms require calculating the amplitude matching coefficient for each sampling point at each detector, matching the seismic detector data, and merging it with the hydrophone data to separate the uplink and downlink wavefields. Classical algorithms first perform water depth correction on the superposition of hydrophone and seismic detector reflection data at each detector point. Then, they sum the superimposed data before and after correction for the hydrophone data and subtract the superimposed data before and after correction for the seismic detector data to obtain superimposed data from both the hydrophone and seismic detectors without downlink waves. Finally, the matching operator, i.e., interactive ghosting, is obtained using these two superimposed data to avoid interference from downlink waves of opposite polarity in the dual-detector matching operator calculation. However, methods that only perform amplitude matching have poor adaptability, and classical algorithms cannot avoid noise interference, thus requiring improvement in accuracy. Summary of the Invention

[0005] This application relates to a method and apparatus for calculating a dual-detection matching operator for seismic data of seafloor nodes, which can improve the accuracy and reliability of the dual-detection matching operator.

[0006] The technical solution is as follows:

[0007] On the one hand, a method for calculating a dual-detection matching operator for seabed node data is provided, the method comprising:

[0008] The range of the direct wave is calculated based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave.

[0009] Linear dynamic correction and common detector point superposition are performed on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range, respectively, to obtain the hydrophone common detector point superposition data and the seismic detector common detector point superposition data;

[0010] Based on the superimposed data of the hydrophone's common detector points and the superimposed data of the seismic detector's common detector points, the dual-detection matching operator for each detector point is calculated using the least squares method within a preset time window.

[0011] On the other hand, a calculation device for a dual-detection matching operator of seabed node data is provided, the device comprising:

[0012] The direct wave range calculation module is used to calculate the direct wave range based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave.

[0013] The data overlay module is used to perform linear dynamic correction and common detector point overlay on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range, respectively, to obtain the hydrophone common detector point overlay data and the seismic detector common detector point overlay data.

[0014] The dual-detection matching operator calculation module is used to calculate the dual-detection matching operator for each detector point within a preset time window using the least squares method, based on the superimposed data of the hydrophone's common detector points and the superimposed data of the seismic detector's common detector points.

[0015] On the other hand, a computer device is provided, which includes a processor and a memory. The memory stores at least one instruction, at least one program, code set, or instruction set. The processor can load and execute at least one instruction, at least one program, code set, or instruction set to implement the seabed node data dual-detection matching operator calculation method provided in the above-mentioned application embodiments.

[0016] On the other hand, a computer-readable storage medium is provided, which stores at least one instruction, at least one program, code set or instruction set, and a processor can load and execute at least one instruction, at least one program, code set or instruction set to implement the seabed node data dual detection matching operator calculation method provided in the above embodiments of this application.

[0017] On the other hand, a computer program product or computer program is provided, which includes computer program instructions stored in a computer-readable storage medium. A processor reads the computer instructions from the computer-readable storage medium and executes the computer instructions, causing the computer device to perform the seabed node data dual-detection matching operator calculation method described in any of the above embodiments.

[0018] The beneficial effects of the technical solution provided in this application include at least the following: This invention provides a method for calculating a dual-detection matching operator for seabed node data. The method includes calculating the direct wave range based on the excitation point depth, the propagation speed of seismic waves in water, the propagation speed of seismic waves on the seabed, and the offset of the first arriving direct wave; performing linear dynamic correction and co-detector point superposition on the hydrophone direct wave and the seismic detector direct wave within the direct wave range, respectively, to obtain hydrophone co-detector point superposition data and seismic detector co-detector point superposition data; and calculating the dual-detection matching operator for each detector point within a preset time window using the least squares method based on the hydrophone co-detector point superposition data and the seismic detector co-detector point superposition data. The method for calculating a dual-detection matching operator for seabed node data provided by this invention can obtain a stable and reliable dual-detection matching operator using direct waves, providing a guarantee for subsequent data processing. Attached Figure Description

[0019] To more clearly illustrate the technical solutions in the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0020] Figure 1 This diagram illustrates the calculation process of the matching operator in the prior art.

[0021] Figure 2 This invention provides a flowchart illustrating a method for calculating a dual-detection matching operator for seabed node data, as provided in an exemplary embodiment of this application.

[0022] Figure 3 This illustration shows a schematic diagram of the direct wave range in a dual-detection matching operator calculation method for seabed node data provided in an exemplary embodiment of this application;

[0023] Figure 4 This illustration shows a schematic diagram of the superposition effect of common detector points in a dual-detection matching operator calculation method for seabed node data provided in an exemplary embodiment of this application.

[0024] Figure 5 This illustration shows a schematic diagram of the calculation result of the dual-detection matching operator in a dual-detection matching operator calculation method for seabed node data provided in an exemplary embodiment of this application;

[0025] Figure 6 This illustration shows another flowchart of a method for calculating a dual-detection matching operator for seabed node data provided in an exemplary embodiment of this application;

[0026] Figure 7 This illustration shows a schematic diagram of the seismic detector data correction effect in a dual-detection matching operator calculation method for seabed node data provided in an exemplary embodiment of this application;

[0027] Figure 8 This paper presents another flowchart of a method for calculating a dual-detection matching operator for seabed node data, provided in an exemplary embodiment of this application.

[0028] Figure 9 This illustration shows a schematic diagram of the wavefield separation effect in a dual-detection matching operator calculation method for seabed node data provided in an exemplary embodiment of this application;

[0029] Figure 10 This illustration shows a comparison of the wavefield separation effect of a dual-detection matching operator calculation method for seabed node data provided by another exemplary embodiment of this application with other methods;

[0030] Figure 11 This invention provides a structural diagram of a dual-detection matching operator computing device for seabed node data, according to an exemplary embodiment of this application.

[0031] Figure 12 This invention provides another structural diagram of a submarine node data dual-detection matching operator computing device according to an exemplary embodiment of the present application;

[0032] Figure 13 This illustration shows a schematic diagram of the structure of a computer device for calculating a dual-detection matching operator for seabed node data, provided in an exemplary embodiment of this application. Detailed Implementation

[0033] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0034] The existing spatiotemporally variable matching algorithm was proposed by Barr Fred J in 1997, and this method is also the main OBN data wavefield separation algorithm currently used in Paradigm processing systems.

[0035] Specifically, the DEGHOST solution for the dual-detector field separation module in the Paradigm system includes:

[0036] For each hydrophone and seismic detector at each receiver point, for a time window centered at sampling point j, slide downwards with a step size of one sampling point, and calculate the average absolute amplitude of the hydrophone within the time window:

[0037]

[0038] in, Let be the average amplitude of the hydrophone data within the time window, N be the number of sampling points, and (t1-t2) be the calculation time window. For time t i The amplitude of the hydrophone data at the location.

[0039] Calculate the average absolute amplitude of the seismic detector within the time window:

[0040]

[0041] in, Let be the average amplitude of the seismic detector data within the time window, N be the number of sampling points, and (t1-t2) be the calculation time window. For time t i The amplitude of the seismic detector data at the location.

[0042] Calculate the wave field (upward wave) after suppressing the ghost wave:

[0043]

[0044] Where U represents the ascending wave field, H represents hydrophone data, and G represents seismic detector data.

[0045] Calculate the downwave field:

[0046]

[0047] Where D represents the downlink wave field, H represents hydrophone data, and G represents seismic detector data.

[0048] Matching algorithms can find spatiotemporally varying matching operators, but this scheme only performs amplitude matching and does not consider factors such as phase, frequency, and time difference. This scheme can only achieve good wavefield separation when the hydrophone and seismic detector have the same frequency and phase (upward wave) or opposite phase (downward wave), and no time difference. Since data meeting these stringent requirements is rare in actual seismic data acquisition, this scheme has poor adaptability.

[0049] The existing classic algorithm is the classic double-detection matching algorithm proposed by Soubaras in 1996. The classic algorithm can provide a stable method for obtaining matching operators and is suitable for processing extremely shallow water seismic data.

[0050] Figure 1 A schematic diagram of the process of calculating the matching operator using a classical algorithm in the prior art is shown.

[0051] See Figure 1 The hydrophone and seismic detector reflection data are superimposed at the same detector point to obtain... Figure 1 The data includes hydrophone data H and seismic detector data G. The seismic detector reflection data does not include the first arrival wave.

[0052] Water depth corrections were applied downwards to the superimposed data of the hydrophone and seismic detector reflection data at the same detector point, resulting in... Figure 1 The corrected hydrophone data ZH and the corrected seismic detector data ZG are shown in the figure.

[0053] The summation of the superimposed data before and after hydrophone correction yields (1+Z)H=H′;

[0054] The difference between the superimposed data before and after the seismic detector correction is calculated to obtain (1-Z)G=G′.

[0055] The matching operator is calculated using the least squares method:

[0056] f c =argmin‖H′-fG′‖2;

[0057] Among them, f c For the matching operator, f is the matching operator solved using least squares, H′ is the sum of the superimposed data before and after hydrophone correction, and G′ is the sum of the superimposed data before and after seismic detector correction.

[0058] Finally, wavefield separation is performed:

[0059]

[0060] Where U represents the ascending wave field, D represents the descending wave field, H represents hydrophone data, G represents seismic detector data, and f represents the descending wave field. c For matching operators.

[0061] Classical algorithms utilize superimposed data after interactive ghosting to derive the matching operator, thus avoiding interference from downlink waves of opposite polarity in the matching operator deriving process for dual-detection systems. While this method is applicable to various dual-detection combined wavefield separation scenarios and can provide a stable matching operator deriving method when the signal-to-noise ratio of the data is sufficient, even after interactive ghosting, the classical algorithm still cannot completely eliminate interference from downlink waves and other noise in the matching operator deriving process.

[0062] Example 1

[0063] Figure 2 A flowchart illustrating a method for calculating a dual-detection matching operator for seabed node data provided in an exemplary embodiment of this application is shown.

[0064] See Figure 2 The method for calculating the dual-detection matching operator for seabed node data provided in this application may include steps 101 to 103.

[0065] In seabed node data acquisition, if the water depth reaches a certain value, a stable direct wave signal with a certain offset distance can be obtained. Based on water depth data, the propagation speed of seismic waves on the seabed, and the propagation speed of seismic waves in water, the offset distance range of the direct wave signal can be obtained.

[0066] Among them, the difference wavelet of the wave field of the land-sea detector can be replaced by the difference wavelet of the direct wave.

[0067] A matching operator that directly calculates hydrophone and seismic detector data based on direct waves.

[0068] In marine exploration, if a certain water depth is reached, there will be a certain time interval between the direct wave and its corresponding seabed. At this time, selecting the direct wave with an hour window can avoid the cross-correlation problem of noisy signals in the classic algorithm.

[0069] Step 101: Calculate the range of the direct wave based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave.

[0070] In some embodiments, step 101 includes: calculating the offset range of the first arrival direct wave based on the offset calculation formula.

[0071] Specifically, the formula for calculating the offset includes:

[0072]

[0073] Where h is the water depth at the excitation point, and V w V is the speed at which seismic waves propagate in water. b denoted as , where is the propagation speed of the seismic wave on the seabed, and is the first arrival distance and the offset distance of the direct wave.

[0074] The offset range for the first arrival wave to be a direct wave is determined based on the offset calculation formula.

[0075] Furthermore, based on the initial arrival direct wave offset range, the range of the direct wave is determined to be 0 to off.

[0076] In this embodiment, since the direct wave is not contaminated by noise, not many seismic channels are needed for the superposition of the co-detectors.

[0077] Figure 3 This is a schematic diagram of the direct wave range, showing the direct wave offset distance calculated based on the water depth at the excitation point, the offset velocity of the seismic wave in the water, and the propagation velocity of the seismic wave on the seabed, and the direct wave range determined based on the direct wave offset distance.

[0078] Figure 3 In the diagram, the horizontal axis represents the channel number, and the vertical axis represents time, in seconds. A suitable direct wave range is selected from the common receiver gather based on the direct wave range. Figure 3 The part of the curve that is framed.

[0079] Step 102: Perform linear dynamic correction and common detector point superposition on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range respectively to obtain the hydrophone common detector point superposition data and the seismic detector common detector point superposition data.

[0080] In some embodiments, linear dynamic correction is performed on the direct waves from the hydrophone and the direct waves from the seismic detector based on the water velocity.

[0081] Specifically, linear dynamic correction uses a constant velocity to dynamically correct the target data, thereby smoothing out the signal at a certain velocity.

[0082] Optional, the water velocity is 1500m / s.

[0083] Furthermore, the direct waves from the hydrophone after linear dynamic correction and the direct waves from the seismic detector after linear dynamic correction are superimposed using common detector points.

[0084] Specifically, common detector point superposition means superimposing data from the same detector point, similar to common center point superposition.

[0085] In this embodiment, the superposition of common detector points can improve the signal-to-noise ratio of the data by superimposing multiple channels, and can also avoid occasional single-channel anomalies.

[0086] Figure 4 The data shown are hydrophone data and seismograph data after the superposition of common detector points.

[0087] Figure 4In the graph, the horizontal axis represents the detector point number, and the vertical axis represents time, in seconds. The left side shows the hydrophone data after stacking the common detector points, and the right side shows the seismic detector data after stacking the common detector points.

[0088] Step 103: Based on the superimposed data of the hydrophone co-detector points and the superimposed data of the seismic detector co-detector points, calculate the dual-detector matching operator for each detector point using the least squares method within a preset time window.

[0089] Step 103 includes: calculating the double-detection matching operator based on the double-detection matching operator calculation formula.

[0090] The calculation formula for the double-detection matching operator includes:

[0091] f c =argmin‖H1-fG1‖2

[0092] Among them, f c Here, H1 is the superimposed data of the common detector points of the hydrophone, G1 is the superimposed data of the common detector points of the seismic detector, and f is the preset matching operator in the least squares solution.

[0093] Figure 5 A schematic diagram of the dual-detection matching operator provided in this embodiment is shown.

[0094] Figure 5 In the graph, the horizontal axis represents the detector point number, and the total axis represents time, in seconds. Figure 5 In this process, each channel is a wavelet, which is used to perform data convolution with seismic detectors at the same receiver point to match the data.

[0095] The method for calculating the dual-detection matching operator for seafloor node seismic data provided in this embodiment of the invention is based on direct wave-based dual-detection matching operator calculation, which avoids noise interference and improves the accuracy and reliability of the dual-detection matching operator.

[0096] Example 2

[0097] The method for calculating the dual-detection matching operator for seabed node data provided in this application may include steps 101 to 103. Before step 103, it may also include a step of time window calculation.

[0098] Step 101: Calculate the range of the direct wave based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave.

[0099] In some embodiments, step 101 includes: calculating the offset range of the first arrival direct wave based on the offset calculation formula.

[0100] Specifically, the formula for calculating the offset includes:

[0101]

[0102] Where h is the water depth at the excitation point, and V w V is the speed at which seismic waves propagate in water. b denoted as , where is the propagation speed of the seismic wave on the seabed, and is the first arrival distance and the offset distance of the direct wave.

[0103] Furthermore, based on the initial arrival direct wave offset range, the range of the direct wave is determined to be from 0 to off.

[0104] In this embodiment, since the direct wave is not contaminated by noise, not many seismic traces are needed for co-detector point superposition.

[0105] Step 102: Perform linear dynamic correction and common detector point superposition on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range respectively to obtain the hydrophone common detector point superposition data and the seismic detector common detector point superposition data.

[0106] In some embodiments, linear dynamic correction is performed on the direct waves from the hydrophone and the direct waves from the seismic detector based on the water velocity.

[0107] Specifically, linear dynamic correction uses a constant velocity to dynamically correct the target data, thereby flattening the signal at a certain velocity.

[0108] Optional, the water velocity is 1500m / s.

[0109] Furthermore, the direct waves from the hydrophone after linear dynamic correction and the direct waves from the seismic detector after linear dynamic correction are superimposed using common detector points.

[0110] Specifically, common detector point superposition means superimposing data from the same detector point, similar to common center point superposition.

[0111] In this embodiment, the superposition of common detector points can improve the signal-to-noise ratio of the data by superimposing multiple channels, and can also avoid occasional single-channel anomalies.

[0112] Step 103: Based on the superimposed data of the hydrophone co-detector points and the superimposed data of the seismic detector co-detector points, calculate the dual-detector matching operator for each detector point using the least squares method within a preset time window.

[0113] Step 103 includes: calculating the double-detection matching operator based on the double-detection matching operator calculation formula.

[0114] The calculation formula for the double-detection matching operator includes:

[0115] f c =argmin‖H1-fG1‖2

[0116] Among them, fc Here, H1 represents the stacked data of the hydrophone's common detector points, G1 represents the stacked data of the seismic detector's common detector points, and f represents the preset matching operator in the least squares solution.

[0117] Before step 103, the method for calculating the dual-detection matching operator for seafloor node seismic data provided in this embodiment may further include:

[0118] The time at the center point of the time window is calculated based on the water depth at the receiving point and the propagation speed of seismic waves in the water.

[0119] The preset time window is determined based on the center point time of the time window and the preset time window length.

[0120] The formula for calculating the center point time of the time window includes:

[0121] t m =2h r / V w

[0122] Among them, t m h is the center point time of the time window. r V is the water depth at the receiving point. w This represents the speed at which seismic waves propagate in water.

[0123] Furthermore, the time at the center point of the time window is calculated based on the formula for calculating the time at the center point of the time window.

[0124] Specifically, the time window is t m -τ1 to t m +τ2, the time window length is τ1+τ2.

[0125] Optionally, the time window length can range from 100ms to 300ms.

[0126] In this embodiment, by selecting an hourly window, interference from multiple underwater noises and other noises can be avoided.

[0127] The dual-detection matching operator calculation method for seafloor node seismic data provided in this embodiment of the invention is based on direct wave-based dual-detection matching operator calculation, and further reduces noise interference and improves the accuracy and reliability of the dual-detection matching operator by selecting an appropriate time window.

[0128] Example 3

[0129] Figure 6 This paper presents another flowchart of a method for calculating a dual-detection matching operator for seabed node data provided in an exemplary embodiment of this application.

[0130] See Figure 6The method for calculating the dual-detection matching operator for seabed node data provided in this application may include steps 201 to 205.

[0131] Step 201: Calculate the range of the direct wave based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave.

[0132] In some embodiments, step 201 includes: calculating the offset range of the first arrival direct wave based on the offset calculation formula.

[0133] Specifically, the formula for calculating the offset includes:

[0134]

[0135] Where h is the water depth at the excitation point, and V w V is the speed at which seismic waves propagate in water. b denoted as , where is the propagation speed of the seismic wave on the seabed, and is the first arrival distance and the offset distance of the direct wave.

[0136] Furthermore, based on the initial arrival direct wave offset range, the range of the direct wave is determined to be from 0 to off.

[0137] In this embodiment, since the direct wave is not contaminated by noise, not many seismic channels are needed for the superposition of the co-detectors.

[0138] Step 202: Perform linear dynamic correction and common detector point superposition on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range respectively to obtain the hydrophone common detector point superposition data and the seismic detector common detector point superposition data.

[0139] In some embodiments, linear dynamic correction is performed on the direct waves from the hydrophone and the direct waves from the seismic detector based on the water velocity.

[0140] Specifically, linear dynamic correction uses a constant velocity to dynamically correct the target data, thereby flattening the signal at a certain velocity.

[0141] Optional, the water velocity is 1500m / s.

[0142] Furthermore, the direct waves from the hydrophone after linear dynamic correction and the direct waves from the seismic detector after linear dynamic correction are superimposed using common detector points.

[0143] Specifically, common detector point superposition means superimposing data from the same detector point, similar to common center point superposition.

[0144] In this embodiment, the superposition of common detector points can improve the signal-to-noise ratio of the data by superimposing multiple channels, and can also avoid occasional single-channel anomalies.

[0145] Step 203: Based on the superimposed data of the hydrophone co-detector points and the superimposed data of the seismic detector co-detector points, calculate the dual-detector matching operator for each detector point using the least squares method within a preset time window.

[0146] Step 203 includes: calculating the double-detection matching operator based on the double-detection matching operator calculation formula.

[0147] The calculation formula for the double-detection matching operator includes:

[0148] f c =argmin‖H1-fG1‖2

[0149] Among them, f c Here, H1 is the superimposed data of the common detector points of the hydrophone, G1 is the superimposed data of the common detector points of the seismic detector, and f is the preset matching operator in the least squares solution.

[0150] Before step 203, the method for calculating the dual-detection matching operator for seafloor node seismic data provided in this embodiment may further include:

[0151] The time at the center point of the time window is calculated based on the water depth at the receiving point and the propagation speed of seismic waves in the water.

[0152] The preset time window is determined based on the center point time of the time window and the preset time window length.

[0153] The formula for calculating the center point time of the time window includes:

[0154] t m =2h r / V w

[0155] Among them, t m h is the center point time of the time window. r V is the water depth at the receiving point. w This represents the speed at which seismic waves propagate in water.

[0156] Furthermore, the time at the center point of the time window is calculated based on the formula for calculating the time at the center point of the time window.

[0157] Specifically, the time window is t m -τ1 to t m +τ2, the time window length is τ1+τ2.

[0158] Optionally, the time window length can range from 100ms to 300ms.

[0159] In this embodiment, by selecting an hourly window, interference from multiple underwater noises and other noises can be avoided.

[0160] In some embodiments, after step 203, the method may further include steps 204 and 205 for data correction and wavefield separation based on the dual-detection matching operator.

[0161] Step 204: Acquire hydrophone data and seismic detector data, and correct the seismic detector data based on the dual-detection matching operator to obtain the corrected seismic detector data.

[0162] By calibrating the detector data, the data from the hydrophone and the seismic detector are matched.

[0163] Figure 7 This diagram illustrates the effect of the dual-detection matching operator correcting the seismic detector data in this embodiment. The left side shows the hydrophone and seismic detector data before matching, and the right side shows the hydrophone and seismic detector data after matching.

[0164] Figure 7 In the diagram, the horizontal axis represents the track number, and the vertical axis represents time, with the unit being seconds. The direct waves from the hydrophone and the direct waves from the seismic detector are arranged alternately.

[0165] See Figure 7 Before matching, there were significant amplitude and phase differences between the hydrophone data and the seismic detector data. After matching, the direct wave amplitude, frequency, phase, and time of the two detector data were consistent. The dual-detection matching operator obtained based on the seabed node seismic data dual-detection matching operator calculation method provided in this embodiment of the invention can effectively correct the seismic detector data and achieve matching between the hydrophone data and the seismic detector data.

[0166] Step 205: Based on the hydrophone data and the corrected seismic detector data, perform wavefield separation.

[0167] In some embodiments, wavefield separation is performed based on wavefield separation calculation formulas.

[0168] The wavefield separation calculation formula includes:

[0169]

[0170] Where U represents the ascending wave field, D represents the descending wave field, H represents hydrophone data, G represents seismic detector data, and f represents the descending wave field. c This is a matching operator derived based on direct waves.

[0171] In classic algorithms, after calculating the matching operator using reflection data, the hydrophone data and seismic detector data are summed to obtain the upgoing wave, and the difference is calculated to obtain the downgoing wave. However, in the dual-detection matching operator calculation method for seafloor node seismic data provided in this embodiment, after correcting the seismic detector data using the dual-detection matching operator, the difference between the hydrophone data and the seismic detector data is calculated to obtain the upgoing wave, and the sum is calculated to obtain the downgoing wave.

[0172] The method provided in this embodiment of the invention is applied to OBN acquisition, where the seabed depth reaches a point where direct waves can be selected to calculate the data matching operator for hydrophones and seismic detectors.

[0173] The method for calculating the dual-detection matching operator for seabed node data provided in this embodiment of the invention can use direct waves to obtain a stable dual-detection matching operator, thereby improving the reliability and accuracy of uplink and downlink wave field separation.

[0174] Example 4

[0175] Figure 8 This paper presents another flowchart of a method for calculating a dual-detection matching operator for seabed node data provided in an exemplary embodiment of this application.

[0176] The process for calculating the double-detection matching operator is as follows:

[0177] First, the direct waves from the hydrophone and the seismic detector are obtained. Then, linear dynamic correction and common detector point superposition are performed on the direct waves from the hydrophone and the seismic detector, respectively.

[0178] Based on the direct wave from the hydrophone and the direct wave from the seismic detector after correction and co-detector point superposition, the dual-detector matching operator f is calculated. c .

[0179] Seismic data from seismic detectors and hydrophones are preprocessed separately, and the preprocessed seismic data from seismic detectors is corrected based on a dual-detection matching operator.

[0180] Wavefield separation is performed based on the corrected seismic detector data and the preprocessed hydrophone seismic data to obtain the upgoing and downgoing waves.

[0181] Figure 9 This diagram illustrates another application effect of the dual-detection matching operator calculation method for seabed node data provided in this embodiment of the invention.

[0182] Figure 9 In the graph, the horizontal axis represents the detector point number, and the vertical axis represents time, in seconds.

[0183] Figure 9 In the middle, from left to right, are hydrophone data, seismic detector data, downwave records, and upwave records.

[0184] See Figure 9 The matching operator obtained by this invention can achieve good separation results for wave field separation.

[0185] Figure 10 This diagram illustrates another application effect of the dual-detection matching operator calculation method for seabed node data provided in this embodiment of the invention.

[0186] Figure 10 In the graph, the horizontal axis represents the track number, and the vertical axis represents time, with the unit being seconds.

[0187] Figure 10 The images, from left to right, show the wavefield separation results obtained by two methods: direct merging of detector data, the Paradigm processing method, the Soubaras method, and the seabed node dual-detector matching operator calculation method provided in this embodiment. In each set of images, the left side shows the dual-detector record, and the right side shows the uplink wavefield after wavefield separation. For the direct merging method, the dual-detector record is not matched; for the Paradigm processing method, the Soubaras method, and the method provided in this embodiment, the dual-detector record is matched.

[0188] See Figure 10 The wave field separation effect of directly merging the dual detection records is not good, and there is still obvious downflow wave energy in the upflow wave field.

[0189] For the Paradigm processing method, the energy recorded by the two detectors can be matched well. However, since this technique does not take into account the frequency and phase differences between the two detector data, the wave field separation effect is still not good.

[0190] For the Soubaras method, the amplitude, frequency, phase, and time difference of the two detector data were matched, and the wave field separation effect was good. The in-phase axis of the downwave in the upwave record was weak, but there is still room for improvement.

[0191] The method for calculating the dual-detection matching operator for seabed node data provided in this embodiment of the invention uses direct wave to calculate the matching operator for correction before performing dual-detection merging to obtain a clean upgoing wave field without observing downgoing wave energy.

[0192] Compared to the Paradim processing method, the method provided in this embodiment of the invention can not only correct the amplitude of dual-detection data, but also comprehensively consider the differences between the two detectors in amplitude, frequency, phase and time difference.

[0193] Compared to the Soubaras method, which is more suitable for extremely shallow seas, the method provided in this embodiment of the invention avoids the cross-correlation problem of noisy signals by selecting an hour window, making the dual-detection matching operator more stable and reliable. Furthermore, it can directly obtain the matching operator by superimposing the common detector points of the hydrophone and the seismic detector, without the need for the step of interactive ghosting, making the method simpler.

[0194] In summary, the dual-detection matching operator calculation method for seabed node data provided in this application can perform dual-detection matching operator calculation based on direct waves. By selecting an appropriate time window, noise interference is reduced, improving the accuracy and reliability of the dual-detection matching operator. Furthermore, the dual-detection matching operator is used to match seismic detector data, ultimately separating the upgoing and downgoing waves, providing a guarantee for subsequent data analysis.

[0195] Example 5

[0196] Figure 11 A schematic diagram of the computational device for the dual-detection matching operator at seabed nodes is shown. (See also...) Figure 11 The seabed node data dual-detection matching operator calculation device provided in this embodiment of the invention may include modules 301 to 303.

[0197] The direct wave range calculation module 301 is used to calculate the direct wave range based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave.

[0198] Data overlay module 302 is used to perform linear dynamic correction and common detector point overlay on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range, respectively, to obtain hydrophone common detector point overlay data and seismic detector common detector point overlay data;

[0199] The dual-detection matching operator calculation module 303 is used to calculate the dual-detection matching operator for each detector point within a preset time window using the least squares method, based on the superimposed data of the hydrophone common detector points and the superimposed data of the seismic detector common detector points.

[0200] In some embodiments, the direct wave range calculation module 301 is specifically used for:

[0201] The offset range of the first-arrival direct wave is calculated based on the offset calculation formula;

[0202] Based on the initial arrival direct wave offset range, the direct wave range is determined to be [0, off];

[0203] The formula for calculating the offset distance includes:

[0204]

[0205] Where h is the water depth at the excitation point, and V w V is the speed at which seismic waves propagate in water. b denoted as , where is the propagation speed of the seismic wave on the seabed, and 'off' is the offset distance of the first arriving wave.

[0206] In some embodiments, the dual-detection matching operator calculation module 303 is specifically used for:

[0207] The double-detection matching operator is calculated based on the double-detection matching operator calculation formula;

[0208] The calculation formula for the dual-detection matching operator includes: f c =argmin‖H1-fG1‖2;

[0209] Among them, f c Here, H1 is the superimposed data of the common detector points of the hydrophone, G1 is the superimposed data of the common detector points of the seismic detector, and f is the preset matching operator in the least squares solution.

[0210] In some embodiments, the seafloor node seismic data dual-detection matching operator calculation device may further include a time window calculation module, used for:

[0211] The time at the center point of the time window is calculated based on the water depth at the receiving point and the propagation speed of the seismic wave in the water.

[0212] The preset time window is determined based on the center point time of the time window and the preset time window length.

[0213] In some embodiments, the time window calculation module is specifically used for:

[0214] The time at the center point of the time window is calculated based on the formula for calculating the time at the center point of the time window.

[0215] The formula for calculating the time at the center point of the time window includes: t m =2h r / V w ;

[0216] Among them, t m h is the time at the center of the time window. r V is the water depth at the receiving point. w This represents the speed at which seismic waves propagate in water.

[0217] The seabed node data dual-detection matching operator calculation device provided in this embodiment of the invention calculates the dual-detection matching operator based on direct waves, which is simple to calculate and highly accurate.

[0218] Example 6

[0219] Figure 12 Another schematic diagram of the computational device for the dual-detection matching operator of seabed node data is shown. See also Figure 12 The seabed node data dual-detection matching operator calculation device provided in this embodiment of the invention may include modules 401 to 405.

[0220] The direct wave range calculation module 401 is used to calculate the direct wave range based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave.

[0221] Data overlay module 402 is used to perform linear dynamic correction and common detector point overlay on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range, respectively, to obtain hydrophone common detector point overlay data and seismic detector common detector point overlay data;

[0222] The dual-detection matching operator calculation module 403 is used to calculate the dual-detection matching operator for each detector point within a preset time window using the least squares method, based on the superimposed data of the hydrophone common detector points and the superimposed data of the seismic detector common detector points.

[0223] The correction module 404 is used to acquire hydrophone data and seismic detector data, and correct the seismic detector data based on the dual-detection matching operator to obtain corrected seismic detector data.

[0224] The wavefield separation module 405 performs wavefield separation based on hydrophone data and corrected seismic detector data.

[0225] In some embodiments, the direct wave range calculation module 401 is specifically used for:

[0226] The offset range of the first-arrival direct wave is calculated based on the offset calculation formula;

[0227] Based on the initial arrival direct wave offset range, the direct wave range is determined to be [0, off];

[0228] The formula for calculating the offset distance includes:

[0229]

[0230] Where h is the water depth at the excitation point, and V w V is the speed at which seismic waves propagate in water. b denoted as , where is the propagation speed of the seismic wave on the seabed, and 'off' is the offset distance of the first arriving wave.

[0231] In some embodiments, the dual-detection matching operator calculation module 403 is specifically used for:

[0232] The double-detection matching operator is calculated based on the double-detection matching operator calculation formula;

[0233] The calculation formula for the dual-detection matching operator includes: f c =argmin‖H1-fG1‖2;

[0234] Among them, f c Here, H1 is the superimposed data of the common detector points of the hydrophone, G1 is the superimposed data of the common detector points of the seismic detector, and f is the preset matching operator in the least squares solution.

[0235] In some embodiments, the seafloor node seismic data dual-detection matching operator calculation device may further include a time window calculation module, used for:

[0236] The time at the center point of the time window is calculated based on the water depth at the receiving point and the propagation speed of the seismic wave in the water.

[0237] The preset time window is determined based on the center point time of the time window and the preset time window length.

[0238] In some embodiments, the time window calculation module is specifically used for:

[0239] The time at the center point of the time window is calculated based on the formula for calculating the time at the center point of the time window.

[0240] The formula for calculating the time at the center point of the time window includes: t m =2h r / V w ;

[0241] Among them, t m h is the time at the center of the time window. r V is the water depth at the receiving point. w This represents the speed at which seismic waves propagate in water.

[0242] In some embodiments, the wavefield separation module 405 is specifically used for:

[0243] Wavefield separation is performed based on wavefield separation calculation formula;

[0244] The wavefield separation calculation formula includes:

[0245] Where U represents the ascending wave field, D represents the descending wave field, H represents hydrophone data, G represents seismic detector data, and f represents the descending wave field. c This is a double-check matching operator.

[0246] In summary, the dual-detection matching operator calculation device for seabed node data provided in this application can perform dual-detection matching operator calculations based on direct waves. By selecting an appropriate time window, noise interference is reduced, improving the accuracy and reliability of the dual-detection matching operator. Furthermore, the dual-detection matching operator is used to match seismic detector data, ultimately separating the upgoing and downgoing waves, providing a guarantee for subsequent data analysis.

[0247] It should be noted that the seabed node data dual detection matching operator device provided in the above embodiments is only an example of the division of the above functional modules. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above.

[0248] Example 7

[0249] Figure 13 This illustration shows a schematic diagram of a computer device for calculating a dual-detection matching operator for seabed node data according to an exemplary embodiment of this application. The computer device includes:

[0250] The processor 501 includes one or more processing cores. The processor 501 executes various functional applications and data processing by running software programs and modules.

[0251] The receiver 502 and transmitter 503 can be implemented as a communication component, which can be a communication chip. Optionally, this communication component can include signal transmission functionality. That is, the transmitter 503 can be used to transmit control signals to the image acquisition device and the scanning device, and the receiver 502 can be used to receive corresponding feedback commands.

[0252] The memory 504 is connected to the processor 501 via the bus 505.

[0253] The memory 504 can be used to store at least one instruction, and the processor 501 is used to execute the at least one instruction to implement steps 201 to 205 in the above method embodiments.

[0254] Step 201: Calculate the range of the direct wave based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave.

[0255] In some embodiments, step 201 includes: calculating the offset range of the first arrival direct wave based on the offset calculation formula.

[0256] Specifically, the formula for calculating the offset includes:

[0257]

[0258] Where h is the water depth at the excitation point, and V w V is the speed at which seismic waves propagate in water. b denoted as , where is the propagation speed of the seismic wave on the seabed, and is the first arrival distance and the offset distance of the direct wave.

[0259] Furthermore, based on the initial arrival direct wave offset range, the range of the direct wave is determined to be from 0 to off.

[0260] In this embodiment, since the direct wave is not contaminated by noise, not many seismic channels are needed for the superposition of the co-detectors.

[0261] Step 202: Perform linear dynamic correction and common detector point superposition on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range respectively to obtain the hydrophone common detector point superposition data and the seismic detector common detector point superposition data.

[0262] In some embodiments, linear dynamic correction is performed on the direct waves from the hydrophone and the direct waves from the seismic detector based on the water velocity.

[0263] Specifically, linear dynamic correction uses a constant velocity to dynamically correct the target data, thereby flattening the signal at a certain velocity.

[0264] Optional, the water velocity is 1500m / s.

[0265] Furthermore, the direct waves from the hydrophone after linear dynamic correction and the direct waves from the seismic detector after linear dynamic correction are superimposed using common detector points.

[0266] Specifically, common detector point superposition means superimposing data from the same detector point, similar to common center point superposition.

[0267] In this embodiment, the superposition of common detector points can improve the signal-to-noise ratio of the data by superimposing multiple channels, and can also avoid occasional single-channel anomalies.

[0268] Step 203: Based on the superimposed data of the hydrophone co-detector points and the superimposed data of the seismic detector co-detector points, calculate the dual-detector matching operator for each detector point using the least squares method within a preset time window.

[0269] Step 203 includes: calculating the double-detection matching operator based on the double-detection matching operator calculation formula.

[0270] The calculation formula for the double-detection matching operator includes:

[0271] f c =argmin‖H1-fG1‖2

[0272] Among them, f c Here, H1 is the superimposed data of the common detector points of the hydrophone, G1 is the superimposed data of the common detector points of the seismic detector, and f is the preset matching operator in the least squares solution.

[0273] Before step 203, the method for calculating the dual-detection matching operator for seafloor node seismic data provided in this embodiment may further include:

[0274] The time at the center point of the time window is calculated based on the water depth at the receiving point and the propagation speed of seismic waves in the water.

[0275] The preset time window is determined based on the center point time of the time window and the preset time window length.

[0276] The formula for calculating the center point time of the time window includes:

[0277] t m =2h r / V w

[0278] Among them, t m h is the center point time of the time window. r V is the water depth at the receiving point. w This represents the speed at which seismic waves propagate in water.

[0279] Furthermore, the time at the center point of the time window is calculated based on the formula for calculating the time at the center point of the time window.

[0280] Specifically, the time window is t m -τ1 to t m +τ2, the time window length is τ1+τ2.

[0281] Optionally, the time window length can range from 100ms to 300ms.

[0282] In this embodiment, by selecting an hourly window, interference from multiple underwater noises and other noises can be avoided.

[0283] In some embodiments, after step 203, the method may further include steps 204 and 205 for data correction and wavefield separation based on the dual-detection matching operator.

[0284] Step 204: Acquire hydrophone data and seismic detector data, and correct the seismic detector data based on the dual-detection matching operator to obtain the corrected seismic detector data.

[0285] By calibrating the seismic detector data, the data from the hydrophone and the seismic detector can be matched.

[0286] Step 205: Perform wavefield separation based on hydrophone data and corrected seismic detector data.

[0287] In some embodiments, wavefield separation is performed based on wavefield separation calculation formulas.

[0288] The wavefield separation calculation formula includes:

[0289]

[0290] Where U represents the ascending wave field, D represents the descending wave field, H represents the hydrophone record, G represents the seismograph record, and f represents the wave field. c This is a matching operator derived based on direct waves.

[0291] Example 8

[0292] This application also provides a computer-readable storage medium storing at least one instruction, at least one program, code set, or instruction set, which can be loaded and executed by a processor to implement the above-described method for calculating the dual-detection matching operator of seabed nodes.

[0293] Example 9

[0294] This application also provides a computer program product or computer program, which includes computer instructions stored in a computer-readable storage medium. A processor of a computer device reads the computer instructions from the computer-readable storage medium and executes the computer instructions, causing the computer device to perform the seabed node dual-detection matching operator calculation method described in any of the above embodiments.

[0295] Optionally, the computer-readable storage medium may include: read-only memory (ROM), random access memory (RAM), solid-state drives (SSDs), or optical discs, etc. The random access memory may include resistive random access memory (ReRAM) and dynamic random access memory (DRAM). The sequence numbers of the embodiments described above are for descriptive purposes only and do not represent the superiority or inferiority of the implementation.

[0296] Those skilled in the art will understand that all or part of the steps of the above embodiments can be implemented by hardware or by a program instructing related hardware. The program can be stored in a computer-readable storage medium, such as a read-only memory, a disk, or an optical disk.

[0297] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the above-described division of functional units and modules is merely an example. In practical applications, the above functions can be assigned to different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiments can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit. Furthermore, the specific names of the functional units and modules are only for easy differentiation and are not intended to limit the scope of protection of this application. The specific working process of the units and modules in the above system can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.

[0298] In the above embodiments, the descriptions of each embodiment have different focuses. For parts that are not described in detail or recorded in a certain embodiment, please refer to the relevant descriptions of other embodiments.

[0299] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0300] In the embodiments provided by this invention, it should be understood that the disclosed apparatus / terminal devices and methods can be implemented in other ways. For example, the apparatus / terminal device embodiments described above are merely illustrative. For instance, the division of modules or units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between devices or units may be electrical, mechanical, or other forms.

[0301] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0302] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0303] If the integrated module / unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, all or part of the processes in the methods of the above embodiments of the present invention can also be implemented by a computer program instructing related hardware. The computer program can be stored in a computer-readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the various embodiments of the present invention, and should all be included within the protection scope of the present invention.

Claims

1. A method for calculating a dual-detection matching operator for seismic data from seafloor nodes, characterized in that, The method includes: The range of the direct wave is calculated based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave. Linear dynamic correction and common detector point superposition are performed on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range, respectively, to obtain the hydrophone common detector point superposition data and the seismic detector common detector point superposition data; Based on the superimposed data of the hydrophone common detector points and the superimposed data of the seismic detector common detector points, the dual-detection matching operator for each detector point is calculated using the least squares method within a preset time window. The steps for calculating the range of the direct wave include: calculating the offset range of the first-arrival direct wave based on the offset calculation formula; Based on the initial arrival direct wave offset range, the direct wave range is determined to be [0, off]; The formula for calculating the offset distance includes: Where h is the water depth at the excitation point, This represents the speed at which seismic waves propagate in water. This represents the speed at which seismic waves propagate on the seabed. This represents the offset distance of the first direct wave.

2. The method according to claim 1, characterized in that, The method of calculating the dual-detector matching operator for each detector point using the least squares method within a preset time window, based on the superimposed data of the hydrophone's common detector points and the superimposed data of the seismic detector's common detector points, includes: The double-detection matching operator is calculated based on the double-detection matching operator calculation formula; The calculation formula for the dual-detection matching operator includes: ; in, For the double-check matching operator of least squares solution, This is the superimposed data of the common detector points of the hydrophone. This is the data stacked from the common detector points of the seismic detector. This is a preset matching operator.

3. The method according to claim 1, characterized in that, Before calculating the dual-detector matching operator for each detector point using the least squares method within a preset time window based on the superimposed data of the hydrophone's common detector points and the superimposed data of the seismic detector's common detector points, the method further includes: The time at the center point of the time window is calculated based on the water depth at the receiving point and the propagation speed of the seismic wave in the water. The preset time window is determined based on the center point time of the time window and the preset time window length.

4. The method according to claim 3, characterized in that, The calculation of the time window center point time based on the water depth at the receiving point and the propagation speed of the seismic wave in the water includes: The time at the center point of the time window is calculated based on the formula for calculating the time at the center point of the time window. The formula for calculating the time at the center point of the time window includes: ; in, The time is the center point of the time window. For the water depth at the receiving point, This represents the speed at which seismic waves propagate in water.

5. The method according to any one of claims 1 to 4, characterized in that, After calculating the dual-detector matching operator for each detector point using the least squares method within a preset time window based on the superimposed data of the hydrophone's common detector points and the superimposed data of the seismic detector's common detector points, the method further includes: Acquire hydrophone data and seismic detector data, correct the seismic detector data based on the dual-detection matching operator, and obtain corrected seismic detector data; Wavefield separation is performed based on the hydrophone data and the corrected seismic detector data.

6. The method according to claim 5, characterized in that, The wavefield separation based on the hydrophone data and the corrected seismic detector data includes: Wavefield separation is performed based on wavefield separation calculation formula; The wavefield separation calculation formula includes: ; ; Where U represents the ascending wave field, D represents the descending wave field, H represents hydrophone data, and G represents seismic detector data. This is a double-check matching operator.

7. A calculation device for a dual-detection matching operator for seafloor node seismic data, characterized in that, The device includes: The direct wave range calculation module is used to calculate the direct wave range based on the water depth at the excitation point, the propagation speed of the seismic wave in the water, the propagation speed of the seismic wave on the seabed, and the offset distance of the first arriving direct wave. The steps for calculating the range of the direct wave include: calculating the offset range of the first-arrival direct wave based on the offset calculation formula; Based on the initial arrival direct wave offset range, the direct wave range is determined to be [0, off]. The formula for calculating the offset distance includes: Where h is the water depth at the excitation point, This represents the speed at which seismic waves propagate in water. This represents the speed at which seismic waves propagate on the seabed. The offset distance of the first direct wave; The data overlay module is used to perform linear dynamic correction and common detector point overlay on the direct wave of the hydrophone and the direct wave of the seismic detector within the direct wave range, respectively, to obtain the hydrophone common detector point overlay data and the seismic detector common detector point overlay data. The dual-detection matching operator calculation module is used to calculate the dual-detection matching operator for each detector point within a preset time window using the least squares method, based on the superimposed data of the hydrophone's common detector points and the superimposed data of the seismic detector's common detector points.

8. A computer device, characterized in that, The computer device includes a processor and a memory, wherein the memory stores at least one instruction, at least one program, code set, or instruction set, and the at least one instruction, at least one program, code set, or instruction set is loaded and executed by the processor to implement the double-detection matching operator calculation method for seafloor node seismic data as described in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The readable storage medium stores at least one instruction, at least one program, code set, or instruction set, wherein the at least one instruction, at least one program, code set, or instruction set is loaded and executed by a processor to implement the double-detection matching operator calculation method for seafloor node seismic data as described in any one of claims 1 to 6.

Images