Variable depth cable data processing method, apparatus, device, medium, and program product

Abstract

The application relates to a variable-depth cable data processing method, device, equipment, medium and program product, and relates to the technical field of oil and gas exploration. The method comprises the following steps: acquiring seismic wave data collected through a variable-depth cable technology; setting correction waves of multiple seismic waves based on offsets of the multiple seismic waves, paths of the multiple seismic waves and directions of the multiple seismic waves; determining a difference between travel times of the correction waves and travel times of reflected waves as a base surface correction formula of two-way travel times; and correcting two-way travel times of the multiple seismic waves based on the base surface correction formula. Through analysis of the seismic wave data collected through the variable-depth cable technology, the base surface correction formula of the marine variable-depth cable seismic data is obtained, so that the shot point and the geophone point are located on the same base surface, thereby improving the result accuracy and imaging quality of conventional processing of the variable-depth cable seismic data.

Description

Technical Field

[0001] This application relates to the field of oil and gas exploration technology, and in particular to a variable depth cable data processing method, apparatus, equipment, medium and program product. Background Technology

[0002] Seismic exploration is an important tool for oil and gas exploration, especially offshore oil and gas exploration. Setting appropriate seismic sources and tow cable immersion depths are crucial to ensuring the quality of the acquired data.

[0003] In related technologies, variable depth cable acquisition can significantly improve the resolution of acquired data compared to conventional flat cable acquisition, effectively improving the quality of seismic data and the accuracy of seismic interpretation.

[0004] However, because the immersion depth of the detector in the variable depth cable acquisition technology changes with the shot-receiver distance, the data source and receiver are not coplanar, and there are depth differences between the source and receiver nodes. This causes the propagation path of the seismic rays to become asymmetrical, which directly affects the accuracy of velocity analysis in the conventional processing of variable depth cable data and the final imaging quality. Summary of the Invention

[0005] This application provides a variable depth cable data processing method, apparatus, equipment, medium, and program product, which can correct the seismic source and receiving point to a unified reference plane, thereby improving the resolution of the results of variable depth cable seismic data processing and the final imaging quality; the technical solution is as follows.

[0006] According to one aspect of this application, a variable depth cable data processing method is provided, the method comprising:

[0007] Seismic wave data acquired using variable depth cable technology is obtained; the seismic wave data includes the two-way travel time of multiple seismic waves, the offset of the multiple seismic waves, the path of the multiple seismic waves, and the direction of the multiple seismic waves; the seismic waves include incident waves and reflected waves; the incident waves are seismic waves between the shot point and the reflection point, and the reflected waves are seismic waves between the reflection point and the receiver point, and the reflected waves are divided into up-going waves and down-going waves;

[0008] Based on the offset distance, path, and direction of the multiple seismic waves, a reference surface correction formula for the two-way travel time is obtained.

[0009] Based on the aforementioned reference surface correction formula, the two-way travel time of the multiple seismic waves is corrected.

[0010] According to one aspect of this application, a data processing apparatus for a variable depth cable is provided, the apparatus comprising:

[0011] The data acquisition module is used to acquire seismic wave data collected by variable depth cable technology. The seismic wave data includes the two-way travel time of multiple seismic waves, the offset of the multiple seismic waves, the path of the multiple seismic waves, and the direction of the multiple seismic waves. The seismic waves include incident waves and reflected waves. The incident wave is the seismic wave between the shot point and the reflection point, and the reflected wave is the seismic wave between the reflection point and the receiver point. The reflected wave is divided into up-going waves and down-going waves.

[0012] The setting module sets correction waves for the multiple seismic waves based on their offset distances, paths, and directions; the correction waves are axially symmetric to the incident waves about the normal to the reflection point.

[0013] The determination module is used to determine the calculation formula for the difference between the travel time of the corrected wave and the travel time of the reflected wave as the reference plane correction formula for the two-way travel time.

[0014] The correction module is used to correct the two-way travel time of the multiple seismic waves based on the reference surface correction formula.

[0015] In some embodiments, the determining module is used for,

[0016] In response to the reflected wave being an upward wave, the formula for correcting the reference surface is determined as follows:

[0017]

[0018] Where Δh is the depth of the detector point, V w Let V be the seawater velocity, X be the gun-receiver distance, T0(X) be the two-way reflected wave propagation time, and V be the velocity of the seawater. rms The velocity is the depositional velocity of the strata above the reflection point.

[0019] In some embodiments, the determining module is used for,

[0020] In response to the reflected wave being a downflow wave, the reference surface correction formula is determined as follows:

[0021]

[0022] Where Δh is the depth of the detector point, V w Let V be the seawater velocity, X be the offset distance, T0(X) be the two-way reflected wave propagation time, and V be the velocity of the seawater. rms The velocity is the depositional velocity of the strata above the reflection point.

[0023] In some embodiments, the correction module is used for,

[0024] Based on the reference plane correction formula, the time difference correction value for the round-trip travel time corresponding to the first offset distance is calculated; the first offset distance is any one of the plurality of offset distances;

[0025] The time difference correction value is applied to the round-trip travel time corresponding to the first offset using Sinc interpolation to correct the round-trip travel time corresponding to the first offset.

[0026] In some embodiments, the apparatus further includes: an inverse correction module, for,

[0027] Based on the corrected two-way travel time of the multiple seismic waves, corrected seismic wave data is generated.

[0028] The corrected seismic wave data is subjected to correlation processing to generate target seismic wave data; the target seismic wave data includes the target two-way travel time of multiple seismic waves, the target offset of the multiple seismic waves, the target path of the multiple seismic waves, and the target direction of the multiple seismic waves.

[0029] Based on the reference plane correction formula, the target offset corresponding to the target round-trip travel time is inversely corrected.

[0030] In some embodiments, the inverse correction module is used for,

[0031] Based on the reference plane correction formula, the negative of the time difference correction value of the round-trip travel time of the target corresponding to the first target offset is calculated; the first target offset is any one of the plurality of target offsets;

[0032] The target round-trip travel time corresponding to the first target offset is inversely corrected by applying the negative number to the target round-trip travel time corresponding to the first target offset using Sinc interpolation.

[0033] According to another aspect of this application, a computer device is provided, the computer device including a processor and a memory, the memory storing at least one instruction, at least one program, code set or instruction set, the at least one instruction, the at least one program, the code set or instruction set being loaded and executed by the processor to implement the variable depth cable data processing method as described above.

[0034] According to another aspect of this application, a computer-readable storage medium is provided, wherein at least one instruction, at least one program, code set, or instruction set is stored therein, wherein the at least one instruction, the at least one program, the code set, or the instruction set is loaded and executed by a processor to implement the variable depth cable data processing method as described above.

[0035] According to another aspect of this application, a computer program product is provided, the computer program product including computer instructions stored in a computer-readable storage medium, wherein a processor reads from the computer-readable storage medium and executes the computer instructions to implement the variable depth cable data processing method described above.

[0036] The technical solutions provided in this application embodiment may have the following beneficial effects:

[0037] This method acquires seismic wave data using variable depth cable technology, sets multiple correction waves for each seismic wave, and determines the reference plane correction formula for the two-way travel time by calculating the difference between the travel time of the correction wave and the travel time of the reflected wave. Based on this formula, the two-way travel times of multiple seismic waves are corrected. Since the shot point and receiver point in the seismic wave data acquired using variable depth cable technology are not coplanar, affecting processing accuracy and imaging quality, this method analyzes the seismic wave data acquired using this technology to derive a reference plane correction formula for the marine variable depth cable seismic data reference plane. This ensures that the shot point and receiver point are located on the same reference plane, thereby improving the accuracy and imaging quality of conventional processing of variable depth cable seismic data. Attached Figure Description

[0038] 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.

[0039] Figure 1 This is a schematic diagram of a conventional maritime flat cable technology provided in an exemplary embodiment of this application;

[0040] Figure 2 This is a schematic diagram of a marine variable depth cable technology provided in an exemplary embodiment of this application;

[0041] Figure 3 This is a flowchart of a variable depth cable data processing method provided in an exemplary embodiment of this application;

[0042] Figure 4 This is a schematic diagram of the geometric relationship for up-traveling wave reference plane correction provided in an exemplary embodiment of this application;

[0043] Figure 5 This is a schematic diagram of the geometric relationship for downlink reference plane correction provided in an exemplary embodiment of this application;

[0044] Figure 6 This is a flowchart of a variable depth cable data processing method provided in another exemplary embodiment of this application;

[0045] Figure 7 This is a flowchart of a variable depth cable data processing method provided in yet another exemplary embodiment of this application;

[0046] Figure 8 This is a flowchart of a variable depth cable data processing method provided in another exemplary embodiment of this application;

[0047] Figure 9 This is a schematic diagram of seismic data from a traveling wave model on a variable-depth cable in the ocean, provided in an exemplary embodiment of this application.

[0048] Figure 10 This is a schematic diagram of upflow reference plane correction for variable depth cable model data provided in an exemplary embodiment of this application;

[0049] Figure 11 This is a schematic diagram of up-wave reference plane correction + reverse reference plane correction for variable depth cable model data provided in an exemplary embodiment of this application;

[0050] Figure 12 This is a schematic diagram of seismic data from a marine variable-depth cable downlink wave model provided in an exemplary embodiment of this application;

[0051] Figure 13 This is a schematic diagram of downwave reference plane correction for variable depth cable model data provided in an exemplary embodiment of this application;

[0052] Figure 14 This is a schematic diagram of downwave reference plane correction + inverse reference plane correction for variable depth cable model data provided in an exemplary embodiment of this application;

[0053] Figure 15 This is a schematic diagram of up-wave reference plane correction + conventional dynamic correction for variable depth cable model data provided in an exemplary embodiment of this application;

[0054] Figure 16 This is a schematic diagram of downwave reference plane correction + conventional dynamic correction for variable depth cable model data provided in an exemplary embodiment of this application;

[0055] Figure 17 This is a schematic diagram of traveling wave model data on a variable depth cable plus conventional dynamic correction provided in an exemplary embodiment of this application;

[0056] Figure 18 This is a schematic diagram of variable-depth cable downlink wave model data plus conventional dynamic correction provided in an exemplary embodiment of this application;

[0057] Figure 19 This is a block diagram illustrating a data processing apparatus for a variable depth cable, as shown in an exemplary embodiment of this application.

[0058] Figure 20 This is a structural block diagram of a computer device provided in an exemplary embodiment of this application.

[0059] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application. Detailed Implementation

[0060] 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.

[0061] Exemplary embodiments will now be described in detail, examples of which are illustrated in the accompanying drawings. When the following description relates to the drawings, unless otherwise indicated, the same numbers in different drawings denote the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this application. Rather, they are merely examples of apparatuses and methods consistent with some aspects of this application as detailed in the appended claims.

[0062] The terminology used in this disclosure is for the purpose of describing particular embodiments only and is not intended to be limiting of the disclosure. The singular forms “a,” “the,” and “the” as used in this disclosure and the appended claims are also intended to include the plural forms unless the context clearly indicates otherwise. It should also be understood that the term “and / or” as used herein refers to and includes any and all possible combinations of one or more of the associated listed items.

[0063] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of related data must comply with the relevant laws, regulations and standards of the relevant countries and regions. For example, the attack operations and other target behaviors involved in this application were all obtained under full authorization.

[0064] It should be understood that although the terms first, second, etc., may be used in this disclosure to describe various information, such information should not be limited to these terms. These terms are only used to distinguish information of the same type from one another. For example, a first parameter may also be referred to as a second parameter without departing from the scope of this disclosure, and similarly, a second parameter may also be referred to as a first parameter. Depending on the context, the word "if" as used herein may be interpreted as "when," "when," or "in response to determination."

[0065] The following is a definition of some terms used in this application:

[0066] Normal Move Out (NMO) correction, also known as normal time difference correction, is used for multi-coverage seismic records. Horizontal stacking is performed at common-depth gathers. Due to the non-zero shot-receiver offset normal time difference, the time-distance curve of the common-depth reflector is a hyperbola. NMO corrects the arrival times of reflected waves from the same interface and point on traces with different shot-receiver offsets to the echo times at the common center point, ensuring that they can be stacked in phase during stacking, forming a stacked trace with prominent reflected wave energy (equivalent to a self-excited and self-received recording trace).

[0067] NMO processing requires velocity parameters. For horizontally layered media, if the selected velocity is correct, the hyperbola of reflections can be corrected to a straight line, allowing for in-phase superposition of traces. Using too high a velocity will result in insufficient correction; conversely, using too low a velocity will lead to overcorrection. Neither of these situations guarantees in-phase superposition. For single-cover records, dynamic correction can be used for shot gather records, directly yielding single-cover seismic profiles.

[0068] Sinc interpolation is a commonly used signal reconstruction method that interpolates discretely sampled signals back into continuous functions. Its core idea is to treat each sample point in the sampling sequence as a convolution of a Sinc function and a Dirac impulse function. By adjusting the factor of the Sinc function, different types of interpolation results can be obtained.

[0069] Sinc interpolation exhibits excellent frequency domain characteristics and theorem preservation, but this requires truncation of the function during interpolation. Furthermore, Sinc interpolation has high computational complexity, forces interpolation within a specific frequency band, and excessive distance between interpolation points and sampling points increases interpolation error. Therefore, in practical applications, optimization of interpolation parameters is necessary to ensure that the interpolation error remains within acceptable limits. Sinc interpolation has been widely applied in signal processing, image processing, and seismic exploration. In this scheme, since the correction value obtained from time difference correction is often not an integer number of samples, interpolation processing of the seismic traces is required, and Sinc interpolation can achieve very good results.

[0070] Snell's Law: In optics physics, Snell's Law is a formula describing the relationship between the angle of incidence and the angle of refraction of light or other waves as they travel from one medium to another. This law gives a constant ratio of the sine of the angle of incidence to the sine of the angle of refraction, and this constant value depends on the incident and refracting media.

[0071] Currently, marine variable-depth cable data acquisition technology has become increasingly mature. Acquisition practice has shown that, compared to conventional horizontal cable acquisition, variable-depth cable acquisition can significantly improve the resolution of acquired data, effectively enhancing seismic data quality and the accuracy of seismic interpretation. However, the key challenge remains: how to correctly process the acquired data.

[0072] Please refer to Figure 1 This illustrates a schematic diagram of a conventional maritime flat-line cable technology provided in an exemplary embodiment of this application; please refer to... Figure 2 This illustrates a schematic diagram of a marine variable-depth cable technology provided in an exemplary embodiment of this application. Figure 1 As shown, in conventional offshore flat-cable technology, the shot point and receiver point are almost coplanar, and the propagation paths of seismic rays are almost symmetrical; Figure 2 As shown, in the case of marine variable depth cable technology, the shot point and the receiver point are not coplanar, and the propagation path of seismic rays is no longer symmetrical.

[0073] Most processing modules are based on the assumption that the shot point and receiver point are located on the same horizontal plane. The most typical example is the common center point assumption of conventional horizontal cable seismic data. However, this assumption no longer exists in variable depth cable data. This directly means that the velocity analysis and imaging processing methods in conventional processing are no longer suitable for processing variable depth cable seismic data.

[0074] Because the immersion depth of the detector in variable depth cable acquisition technology changes with the shot-receiver distance, the data source and receiver are not coplanar, and there are depth differences between the source and receiver nodes. This causes the propagation path of seismic rays to become asymmetrical. Many theoretical assumptions in conventional processing procedures are invalid or algorithms are difficult to implement efficiently, which directly affects multiple processing stages in variable depth cable data processing, such as multiple wave suppression, velocity analysis, broadband processing, dynamic correction, and superposition quality control.

[0075] In response to the aforementioned problems in the processing of variable depth cable data, researchers urgently need a reference plane correction technique suitable for marine variable depth cable data to correct the seismic source and receiver points to a unified reference plane, thereby improving the resolution and final imaging quality of the results from conventional processing of variable depth cable seismic data.

[0076] After calibration to a unified reference plane, the shot point and receiver point are located on the same horizontal plane, equivalent to ordinary flat cable data. This allows conventional processing modules to handle variable-depth cable data effectively. Conversely, if reference plane calibration is not performed and conventional modules are used directly, the accuracy of the processing results is reduced because many theoretical assumptions in the conventional processing algorithm do not hold.

[0077] Please refer to Figure 3 This document illustrates a flowchart of a variable depth cable data processing method provided in an exemplary embodiment of this application. The method is executed by a computer device, such as... Figure 3As shown, the method may include steps 310, 320, 330 and 340.

[0078] Step 310: Acquire seismic wave data using variable depth cable technology; the seismic wave data includes the two-way travel time of multiple seismic waves, the offset of multiple seismic waves, the path of multiple seismic waves, and the direction of multiple seismic waves; the seismic waves include incident waves and reflected waves; the incident wave is the seismic wave between the shot point and the reflection point, and the reflected wave is the seismic wave between the reflection point and the receiver point, and the reflected wave is divided into up-going waves and down-going waves.

[0079] In this embodiment of the application, the computer device can obtain information through, for example... Figure 1 The image shows seismic wave data acquired using variable depth cable technology.

[0080] The seismic wave data can be a collection of data acquired multiple times from multiple variable-depth cables, including the two-way travel time, offset, path, and direction of multiple seismic waves. Based on the two-way travel time and offset of multiple seismic waves, a system can be constructed as follows: Figure 9 or Figure 12 The shown is the orthogonal artillery assembly record.

[0081] Seismic waves include incident waves and reflected waves; incident waves are seismic waves between the shot point and the reflection point, and reflected waves are seismic waves between the reflection point and the receiver point; according to the direction of the reflected waves, reflected waves are divided into ascending waves and descending waves.

[0082] Step 320: Based on the offset of multiple seismic waves, the path of multiple seismic waves, and the direction of multiple seismic waves, set correction waves for multiple seismic waves; the correction waves are axially symmetric to the incident waves about the normal of the reflection point.

[0083] In this embodiment of the application, the computer device can set multiple correction waves for the seismic waves based on the offset distance, path, and direction of the multiple seismic waves.

[0084] The correction wave is set on the side of the reflected wave, and the correction wave and the incident wave are axially symmetric about the normal of the reflection point, so that the shot point and the receiver point of the variable depth cable technology are located on the same reference plane, thereby eliminating the influence of the depth difference between the shot point and the receiver point.

[0085] Step 330: Determine the formula for calculating the difference between the travel time of the corrected wave and the travel time of the reflected wave as the reference surface correction formula for the two-way travel time.

[0086] In this embodiment of the application, the computer device can analyze the geometric relationship between the corrected wave path and the reflected wave path, derive a calculation formula for the difference between the travel time of the corrected wave and the travel time of the reflected wave, and determine the calculation formula as the reference surface correction formula for the two-way travel time.

[0087] The aforementioned reference plane correction formula may include the up-traveling wave reference plane correction formula and the down-traveling wave reference plane correction formula.

[0088] Step 340: Based on the datum correction formula, correct the two-way travel time of multiple seismic waves.

[0089] In this embodiment of the application, based on the reference plane correction formula obtained in step 330, the computer device can calculate the time difference correction value corresponding to multiple round-trip travel times, and correct the round-trip travel time of multiple seismic waves according to the time difference correction value.

[0090] Based on the time difference correction value derived from the reference plane correction formula, time difference correction is performed, and the point-by-point position and shape of the seismic wave data acquired by the variable depth cable technology can be adjusted accordingly, so that the seismic wave data acquired by the variable depth cable technology is basically consistent with the data acquired by the horizontal cable, thereby realizing that the shot point and receiver point are located on the same reference plane.

[0091] After calibration to a unified reference plane, the shot point and receiver point are located on the same horizontal plane, equivalent to seismic wave data acquired by a conventional horizontal cable. This allows conventional processing modules to handle seismic wave data acquired by variable-depth cables effectively. Conversely, if reference plane calibration is not performed and conventional modules are used directly, the accuracy of the processing results will be reduced because many theoretical assumptions in the conventional processing algorithm will not hold.

[0092] In summary, the scheme described in this application acquires seismic wave data using variable depth cable technology, sets multiple correction waves for the seismic waves, and determines the formula for calculating the difference between the travel time of the correction waves and the travel time of the reflected waves as the reference plane correction formula for the two-way travel time. Based on the reference plane correction formula, the two-way travel times of multiple seismic waves are corrected. Since the shot point and receiver point are not coplanar in the seismic wave data acquired using variable depth cable technology, this affects the processing accuracy and imaging quality. This scheme, through analysis of the seismic wave data acquired using variable depth cable technology, derives the reference plane correction formula for the reference plane of marine variable depth cable seismic data, ensuring that the shot point and receiver point are located on the same reference plane, thereby improving the accuracy and imaging quality of conventional processing of variable depth cable seismic data.

[0093] In some embodiments, based on Figure 3 In the embodiment shown, step 330 includes:

[0094] Since the reflected wave is an upward wave, the formula for datum plane correction is determined as follows:

[0095]

[0096] Where Δh is the depth of the detector point, V w Let V be the seawater velocity, X be the gun-receiver distance, T0(X) be the two-way reflected wave propagation time, and V be the velocity of the seawater. rms The velocity of the strata above the reflection point is the superposition velocity.

[0097] Please refer to Figure 4 This illustration shows a geometrical diagram of uptraveling wave reference plane correction provided in an exemplary embodiment of this application. Based on the uptraveling wave propagation ray path and direction information of the variable-depth cable, an accurate formula for uptraveling wave reference plane correction of marine variable-depth cable seismic data is derived:

[0098] In this embodiment of the application, for data with a zero offset, due to self-excitation and self-reception, the formula for calculating the difference between the travel time of the corrected wave and the travel time of the reflected wave is obviously:

[0099]

[0100] Where Δh is the depth of the detector point, V w The speed of the seawater.

[0101] like Figure 4 As shown, for data with a non-zero offset, the angle between the corrected wave and the normal at the receiver point is θ1. The formula for calculating the difference between the travel time of the corrected wave and the travel time of the reflected wave is approximated as follows:

[0102]

[0103] This is based on a certain degree of approximation. In actual marine data collection, since the seabed depth is much greater than the cable depth, the time difference between the corrected wave and the reflected wave can be approximated.

[0104] According to Snell's Law, we have:

[0105]

[0106] Where H is half the gun-receiver distance, that is, half of the gun-receiver distance X, and T ref (H) represents the travel time of the seismic reflected wave at a shot-receiver distance of H. It is the one-way travel time and half of the two-way reflection wave propagation time T0(X). rms The velocity of the strata above the reflection point is the superposition velocity.

[0107] Substituting formulas (A1) and (C1) into formula (B1), the time difference for variable depth cable correction is:

[0108]

[0109] Formula (D1) is the derived reference surface correction formula for up-wave ocean variable depth cable seismic data.

[0110] This application provides a reference plane correction formula when the reflected wave is an upgoing wave. By analyzing the geometric relationship of the propagation ray path of the upgoing wave in a variable depth cable seismic event, a reference plane correction formula for marine variable depth cable seismic data is derived so as to correct the source and receiver to a unified reference plane.

[0111] In some embodiments, based on Figure 3 In the embodiment shown, step 330 includes:

[0112] Since the reflected wave is a downflow wave, the formula for determining the reference plane correction is:

[0113]

[0114] Where Δh is the depth of the detector point, V w Let V be the seawater velocity, X be the offset distance, T0(X) be the two-way reflected wave propagation time, and V be the velocity of the seawater. rms The velocity is the depositional velocity of the strata above the reflection point.

[0115] Please refer to Figure 5 This illustration shows a geometrical diagram of downwave reference plane correction provided in an exemplary embodiment of this application. Based on the downwave propagation ray path and direction information of the variable-depth cable, an accurate formula for downwave reference plane correction of marine variable-depth cable seismic data is derived:

[0116] In this embodiment of the application, for data with a zero offset, due to self-excitation and self-reception, the formula for calculating the difference between the travel time of the corrected wave and the travel time of the reflected wave is obviously:

[0117]

[0118] Where Δh is the depth of the detector point, V w The speed of the seawater.

[0119] like Figure 5 As shown, for data with a non-zero offset, when obtaining the angle θ2 between the corrected wave and the normal of the receiver point, the formula for calculating the difference between the travel time of the corrected wave and the travel time of the reflected wave is approximated as follows:

[0120]

[0121] This is based on a certain degree of approximation. In actual marine data collection, since the seabed depth is much greater than the cable depth, the time difference between the corrected wave and the reflected wave can be approximated.

[0122] According to Snell's Law, we have:

[0123]

[0124] Where H is half the gun-receiver distance, that is, half of the gun-receiver distance X, and T ref (H) represents the travel time of the seismic reflected wave at a shot-receiver distance of H. It is the one-way travel time and half of the two-way reflection wave propagation time T0(X). rms The velocity of the strata above the reflection point is the superposition velocity.

[0125] Substituting formulas (A2) and (C2) into formula (B2), the time difference for variable depth cable correction is:

[0126]

[0127] Formula (D2) is the derived reference surface correction formula for down-wave ocean variable depth cable seismic data.

[0128] This application provides a reference plane correction formula when the reflected wave is a downflow wave. By analyzing the geometric relationship of the propagation ray path of the downflow wave in a variable depth cable seismic event, a reference plane correction formula for marine variable depth cable seismic data is derived so as to correct the source and receiver to a unified reference plane.

[0129] Please refer to Figure 6 This illustrates a flowchart of a variable depth cable data processing method provided in another exemplary embodiment of this application. The method is executed by a computer device, such as... Figure 6 As shown above, Figure 3 Step 340 in the illustrated embodiment can be implemented as steps 340a and 340b.

[0130] Step 340a: Based on the datum plane correction formula, calculate the time difference correction value of the round-trip travel time corresponding to the first offset distance; the first offset distance is any one of multiple offset distances.

[0131] In this embodiment of the application, based on the reference plane correction formula obtained in step 330, the computer device can sequentially calculate the time difference correction value for the round-trip travel time for each record.

[0132] For example, based on the first offset distance, i.e. keeping the offset distance parameter unchanged, the time difference correction value for the round-trip travel time corresponding to the first offset distance is calculated.

[0133] Step 340b: Apply time difference correction value to the round-trip travel time corresponding to the first offset using Sinc interpolation to correct the round-trip travel time corresponding to the first offset.

[0134] In this embodiment of the application, after step 340a calculates the time difference correction value for the round-trip travel time corresponding to the first offset, the computer device can apply the time difference correction value to the round-trip travel time corresponding to the first offset by Sinc interpolation to correct the round-trip travel time corresponding to the first offset.

[0135] Subsequently, based on steps 340a and 340b, the round-trip travel time corresponding to the second offset is corrected.

[0136] Similarly, time difference correction values ​​are calculated and applied to each seismic wave data acquired by the variable depth cable technology, so that the shot point and receiver point of the seismic wave data are located on the same reference plane, thus completing the reference plane correction of the variable depth cable data.

[0137] This application provides a technical solution for correcting the two-way travel time of multiple seismic waves based on a reference plane correction formula. The time difference is corrected sequentially for multiple two-way travel times with the same offset parameter in order to improve the efficiency of reference plane correction.

[0138] Please refer to Figure 7 This illustrates a flowchart of a variable depth cable data processing method provided in yet another exemplary embodiment of this application. The method is executed by a computer device, such as... Figure 7 As shown above, Figure 3 The method in the illustrated embodiment may further include steps 350, 360, and 370.

[0139] Step 350: Generate corrected seismic wave data based on the two-way travel time of the corrected multiple seismic waves.

[0140] In this embodiment of the application, after step 340 corrects the two-way travel time of multiple seismic waves based on the reference surface correction formula, the computer device can generate corrected seismic wave data so that the corrected seismic wave data can be processed by multiple wave suppression, velocity analysis, broadband processing, dynamic correction and superposition quality control using theoretical assumptions in conventional processing procedures.

[0141] Step 360: Perform relevant processing on the corrected seismic wave data to generate target seismic wave data; the target seismic wave data includes the target two-way travel time of multiple seismic waves, the target offset of multiple seismic waves, the target path of multiple seismic waves, and the target direction of multiple seismic waves.

[0142] In this embodiment of the application, the computer device can perform relevant processing on the corrected seismic wave data generated in step 350 to generate target seismic wave data, as needed.

[0143] The target seismic wave data can be similar to seismic wave data, including the target two-way travel time of multiple seismic waves, the target offset of multiple seismic waves, the target path of multiple seismic waves, and the target direction of multiple seismic waves.

[0144] Step 370: Based on the datum plane correction formula, inversely correct the target offset corresponding to the target's two-way travel time.

[0145] In this embodiment of the application, the computer device can perform inverse correction on the target seismic wave data based on the reference surface correction formula obtained in step 330.

[0146] Some processing modules can only be completed on data with a unified reference surface, such as the module for removing ocean multiple interference: the computer equipment can first perform time difference correction on the original seismic wave data according to the reference surface correction formula to obtain corrected seismic wave data; then, based on the module for removing ocean multiple interference, the interference of the corrected seismic wave data is removed on the unified reference surface; then, inverse correction is performed. Although the shape is similar to the original seismic wave data, the seismic wave data that has been processed by the module and then inversely corrected has removed the interference, which is more conducive to improving the signal-to-noise ratio of subsequent processing modules.

[0147] This application provides a technical solution for inverse correction based on the datum plane correction formula. Through inverse datum plane correction, the datum plane correction result can be restored to the original variable depth cable seismic wave data, realizing reversible correction processing and eliminating interference from the original variable depth cable seismic wave data.

[0148] Please refer to Figure 8 This illustrates a flowchart of a variable depth cable data processing method provided in another exemplary embodiment of this application. The method is executed by a computer device, such as... Figure 8 As shown above, Figure 7 Step 370 in the illustrated embodiment can be implemented as steps 370a and 370b.

[0149] Step 370a: Based on the reference plane correction formula, calculate the negative of the time difference correction value of the round-trip travel time of the target corresponding to the first target offset; the first target offset is any one of multiple target offsets.

[0150] In this embodiment of the application, after step 360 performs correlation processing on the corrected seismic wave data to generate the target seismic wave data, the computer device can refer to the above-mentioned reference surface correction formula obtained in step 330. Figure 6 The embodiment shown uses the negative number of the time difference correction value for calculating the round-trip travel time of the target corresponding to the first target offset.

[0151] Step 370b: Apply the inverse number to the target round-trip travel time corresponding to the first target offset using Sinc interpolation to inversely correct the target round-trip travel time corresponding to the first target offset.

[0152] In this embodiment of the application, after step 370a calculates the negative number of the time difference correction value of the target round-trip travel time corresponding to the first target offset, the computer device can apply the negative number to the round-trip travel time corresponding to the first target offset by Sinc interpolation to correct the target round-trip travel time corresponding to the first target offset.

[0153] Subsequently, based on steps 370a and 370b, the target round-trip travel time corresponding to the second target offset is corrected.

[0154] Similarly, the inverse of the time difference correction value is calculated for each seismic wave data acquired by the variable depth cable technology and applied to complete the reverse reference plane correction of the variable depth cable data.

[0155] This application provides a supplementary technical solution for inverse correction based on the datum plane correction formula. The time difference is inversely corrected sequentially for the two-way travel time of multiple targets with the same target offset, so as to improve the efficiency of inverse datum plane correction.

[0156] In summary, this application addresses the issues of low velocity analysis accuracy and low final imaging quality in conventional processing of variable-depth cable seismic data. By analyzing the geometric relationship of the propagation ray paths of uplink / downlink waves in variable-depth cable seismic data, an accurate reference plane correction formula for marine variable-depth cable seismic data is derived. This formula can distinguish between uplink and downlink waves in variable-depth cables and accurately correct the source and receiver points to a unified reference plane, thereby improving the resolution and final imaging quality of the processed variable-depth cable seismic data. Through inverse reference plane correction, the reference plane correction result can be restored to the original input variable-depth cable data, achieving a essentially reversible correction process.

[0157] Based on the methods shown in the above embodiments of this application, this application provides an implementation flow diagram of a variable depth cable data processing method, including the following steps.

[0158] S1: As Figure 9 As shown, the ocean variable-depth cable traveling wave seismic data S(t) is input to the forward model of the wave equation. i ,x j ).

[0159] Where t represents the two-way travel time of the reflected wave, x represents the offset distance; i represents the index of the sampling point along the time direction; j represents the index of the sampling point along the offset distance direction; i = 1, 2...n, n represents the total number of sampling points along the time direction; j = 1, 2...m, m represents the total number of sampling points along the offset distance direction.

[0160] In this embodiment, n = 6000, m = 120, meaning there are 6000 sampling points in the time direction and 120 seismic waves. Seawater velocity V w =1500m / s, the cable depth range of the variable depth cable is 20 meters to 200 meters, the detector spacing is 25 meters, the shot-detector distance range is 25 meters to 3000 meters, and the stacking velocity V rms ={(933ms, 1500m / s), (1733ms, 1740m / s), (2200ms, 2078m / s)}.

[0161] S2: For the up-wave variable depth cable seismic data input in S1 above, calculate the time difference correction value of the two-way travel time, and perform time difference correction on the two-way travel time of the above data so that the shot point and the receiver point are located on the same reference plane.

[0162] S21: Selecting up-wave variable depth cable seismic data S(t) i ,x j A record in ) where the number of sampling points increases while the offset parameter remains constant, for example, S(t i (x1), i = 1, 2…6000, j = 1, using the upward wave reference plane correction formula:

[0163]

[0164]

[0165] For S(t) i x1) The time difference correction value TX(t) is calculated for each sampling point according to the sampling time. i )=Δt(θ), i=1, 2...6000.

[0166] In the embodiments of this application, S(t) i x1) The depth of the variable depth cable ΔH = 20 meters, and the seawater velocity V w =1500m / s, offset distance X = 25m, select the first sample point T0(X) = 2ms, the superposition velocity at this point V rms The time difference correction TX(t1) at this point is obtained by linear interpolation of the time velocity {(933ms, 1500m / s), (1733ms, 1740m / s), (2200ms, 2078m / s)} and then substituting these values ​​into the up-traveling wave reference plane correction formula (G).

[0167] Similarly, the time difference correction values ​​TX(t2)...TX(t) at the 2nd to 6000th sampling points can be calculated respectively. 6000 ), that is, S(t) ix1) The time difference correction value TX(t) of this data. i ), i = 1, 2...6000.

[0168] S22: Using Sinc interpolation, apply the time difference correction value TX(t) to this data. i Time difference correction is performed to ensure that the shot detection point is located on the same reference plane.

[0169] S23: Repeat steps S21 and S22, selecting up-wave variable depth cable seismic data S(t) i ,x j The next record in ), such as S(t) i Given x2), i = 1, 2…6000, j = 2, calculate the time difference correction value TX(t) for each sample point of this data. i This is applied to ensure that the shot detector is located on the same reference plane; the operation is repeated 120 times to achieve S(t) i ,x i All 120 records in this data have been corrected.

[0170] S24: As Figure 10 As shown, the up-wave reference plane correction of the variable depth cable model data is completed through the above steps.

[0171] S25: The reverse reference plane correction of the traveling wave data on the variable depth cable is the reverse process of steps S21 to S23 above. The time difference value calculated in step S21 is negative, and steps S22 to S23 are repeated to complete the correction. Figure 11 As shown.

[0172] S3: As Figure 12 As shown, the input wave equation forward model uses the variable depth cable downlink seismic data S(T) of the ocean. I ,X J ), I=1, 2...N, J=1, 2...M.

[0173] T represents the two-way travel time of the reflected wave, X represents the offset distance, I represents the number of the sampling point along the time direction, J represents the number of the sampling point along the offset distance direction, i = 1, 2...N, N represents the total number of sampling points along the time direction, j = 1, 2...M, M represents the total number of sampling points along the offset distance direction.

[0174] In this embodiment of the application, N = 6000 and M = 120, meaning there are 6000 sampling points in the time direction and 120 seismic waves. Seawater velocity V w =1500m / s, the cable depth range of the variable depth cable is 20 meters to 200 meters, the detector spacing is 25 meters, the shot-detector distance range is 25 meters to 3000 meters, and the stacking velocity V rms={(933ms, 1500m / s), (1733ms, 1740m / s), (2200ms, 2078m / s)}.

[0175] S4: For the down-wave variable-depth cable seismic data input in S3 above, calculate the time difference correction value of the two-way travel time, and perform time difference correction on the two-way travel time of the above data so that the shot point and receiver point are located on the same reference plane.

[0176] S41: Select downlink wave depth-variable cable seismic data S(T) I ,T J A record in ) that is, the number of sampling points increases while the offset parameter remains unchanged, such as S(T) I X1), I = 1, 2…6000, J = 1, using the downlink reference plane correction formula:

[0177]

[0178]

[0179] For S(T) I X1) The time difference correction value TX(t) is calculated for each sampling point according to the sampling time. I )=Δt(θ), I=1, 2...6000.

[0180] In the embodiments of this application, S(T) I X1) The depth of the variable depth cable ΔH = 20 meters, and the seawater velocity V w =1500m / s, offset distance X = 25m, select the first sample point T0(X) = 2ms, the superposition velocity at this point V rms The time difference correction TX(t1) at this point can be obtained by linear interpolation of the time velocity on the basis of {(933ms, 1500m / s), (1733ms, 1740m / s), (2200ms, 2078m / s)}. Substituting the above values ​​into the downlink wave reference plane correction formula (H), the time difference correction at this point can be obtained.

[0181] Similarly, the time difference correction values ​​TX(t2)...TX(t) at the 2nd to 6000th sampling points can be calculated respectively. 6000 ), that is, the time difference correction value TX(t) for that channel. I ), I = 1, 2...6000.

[0182] S42: Using Sinc interpolation, apply the time difference correction value TX(t) to this data. I Time difference correction is performed to ensure that the shot detection point is located on the same reference plane.

[0183] S43: Repeat steps S41 and S42, selecting downlink variable depth cable seismic data S(T) I The next record in (X1), such as S(T) I Given X2), I = 1, 2…6000, J = 2, calculate the time difference correction value TX(t) for each sample point of this data. I This process is then applied to ensure that the shot and receiver points are located on the same reference plane. This operation is repeated M times to correct all M records of the data.

[0184] S44: As Figure 13 As shown, the downlink reference plane correction of the variable depth cable model data is completed through the above steps.

[0185] S45: The reverse reference plane correction of the traveling wave data on the variable depth cable is the reverse process of steps S41 to S43 above. The time difference value calculated in step S41 is negativeed, and steps S42 to S43 are repeated to complete the correction. Figure 14 As shown.

[0186] In this embodiment of the application, conventional dynamic correction was performed on each result data to verify the effect of the up-going wave reference plane correction and the down-going wave reference plane correction.

[0187] like Figure 15 As shown, after correcting the uptraveling wave model data of the variable depth cable using the uptraveling wave reference plane, the results of conventional dynamic correction are then performed; for example... Figure 16 As shown, after correcting the downlink reference plane of the downlink wave model data for the variable depth cable, conventional dynamic correction is then performed; for example... Figure 17 As shown, for the uptraveling wave model data of the variable depth cable, no uptraveling wave reference plane correction was performed; the conventional dynamic correction results were directly applied. Figure 18 As shown, the downwave model data of the variable depth cable was directly subjected to conventional dynamic correction without downwave reference plane correction.

[0188] contrast Figure 15 , Figure 16 , Figure 17 and Figure 18 As can be seen, due to the influence of the depth difference between the shot point and the receiver point, direct conventional dynamic correction cannot flatten the original uplink or downlink channel set of the variable depth cable model data. However, after the uplink reference plane correction processing and downlink reference plane correction processing of the embodiments of this application, the shot point and receiver point of the variable depth cable data are located on the same reference plane, effectively eliminating the influence of the depth difference between the shot and receiver points. Therefore, conventional dynamic correction can also make the uplink or downlink phase axis of the variable depth cable data correctly flattened.

[0189] For the case of reverse datum plane correction, compare with the attached... Figure 9 With appendix Figure 11Appendix Figure 12 With appendix Figure 14 It can be seen that by inputting variable depth cable data after up-wave reference plane correction and variable depth cable data after down-wave reference plane correction, and then performing up-wave reverse reference plane correction and down-wave reverse reference plane correction respectively, the original input variable depth cable data can be restored. Therefore, this reference plane correction technology is reversible.

[0190] Since some processing modules can only be completed on data with a unified reference plane, such as the module for removing ocean multiple wave interference, after removing the interference on the unified reference plane, it is back-corrected. Although the data is similar in form to the original data, the interference has been removed, which is more conducive to improving the signal-to-noise ratio of subsequent processing modules.

[0191] In summary, this application can accurately achieve reversible reference level correction for marine variable depth cable seismic data.

[0192] This application addresses the challenge of inconsistent shot and receiver planes in variable-depth cable seismic data acquisition, which negatively impacts processing accuracy and imaging quality during conventional processing. It introduces a reference plane correction method for marine variable-depth cable seismic data. This technical solution analyzes the ray path and direction information acquired by the variable-depth cable and derives reference plane correction formulas for both uplink and downlink waves. Utilizing information such as velocity, offset, water velocity, and cable depth, the variable-depth cable data can be accurately corrected to sea level, ensuring that the shot and receiver planes are located on the same reference plane. This improves the accuracy and imaging quality of conventional processing of variable-depth cable seismic data. Through reverse reference plane correction, the correction results can be restored to the original input variable-depth cable data, achieving essentially reversible correction processing.

[0193] Figure 19 This application illustrates a block diagram of a data processing apparatus for a variable depth cable, as shown in an exemplary embodiment, which can be used to perform tasks such as... Figure 3 , Figure 6 , Figure 7 or Figure 8 In the method shown, all or part of the steps performed by the computer device are as follows: Figure 19 As shown, the device includes:

[0194] The data acquisition module 1901 is used to acquire seismic wave data collected by variable depth cable technology. The seismic wave data includes the two-way travel time of multiple seismic waves, the offset of multiple seismic waves, the path of multiple seismic waves, and the direction of multiple seismic waves. The seismic waves include incident waves and reflected waves. The incident wave is the seismic wave between the shot point and the reflection point, and the reflected wave is the seismic wave between the reflection point and the receiver point. The reflected wave is divided into up-going waves and down-going waves.

[0195] The setting module 1902 is used to set correction waves for multiple seismic waves based on the offset of multiple seismic waves, the path of multiple seismic waves, and the direction of multiple seismic waves; the correction waves are axially symmetric to the incident waves about the normal of the reflection point.

[0196] The determination module 1903 is used to determine the calculation formula for the difference between the travel time of the corrected wave and the travel time of the reflected wave as the reference surface correction formula for the two-way travel time.

[0197] The correction module 1904 is used to correct the two-way travel time of multiple seismic waves based on the reference surface correction formula.

[0198] In some embodiments, the determining module 1903 is used for,

[0199] Since the reflected wave is an upward wave, the formula for datum plane correction is determined as follows:

[0200]

[0201] Where Δh is the depth of the detector point, V w Let V be the seawater velocity, X be the gun-receiver distance, T0(X) be the two-way reflected wave propagation time, and V be the velocity of the seawater. rms The velocity of the strata above the reflection point is the superposition velocity.

[0202] In some embodiments, the determining module 1903 is used for,

[0203] Since the reflected wave is a downflow wave, the formula for determining the reference plane correction is:

[0204]

[0205] Where Δh is the depth of the detector point, V w Let V be the seawater velocity, X be the offset distance, T0(X) be the two-way reflected wave propagation time, and V be the velocity of the seawater. rms The velocity is the depositional velocity of the strata above the reflection point.

[0206] In some embodiments, the correction module 1904 is used for,

[0207] Based on the datum correction formula, the time difference correction value for the round-trip travel time corresponding to the first offset is calculated; the first offset is any one of multiple offsets.

[0208] The round-trip travel time corresponding to the first offset is corrected by applying a time difference correction value to the round-trip travel time using Sinc interpolation.

[0209] In some embodiments, the apparatus further includes: an inverse correction module, for,

[0210] Based on the two-way travel time of multiple corrected seismic waves, corrected seismic wave data is generated.

[0211] The corrected seismic wave data is processed to generate target seismic wave data. The target seismic wave data includes the target two-way travel time of multiple seismic waves, the target offset of multiple seismic waves, the target path of multiple seismic waves, and the target direction of multiple seismic waves.

[0212] Based on the datum plane correction formula, the target offset corresponding to the target's two-way travel time is inversely corrected.

[0213] In some embodiments, the inverse correction module is used for,

[0214] Based on the reference plane correction formula, the negative of the time difference correction value of the round-trip travel time of the target corresponding to the first target offset is calculated; the first target offset is any one of multiple target offsets;

[0215] The target two-way travel time corresponding to the first target offset is inversely corrected by applying the negative number to the target two-way travel time corresponding to the first target offset using Sinc interpolation.

[0216] Figure 20 A structural block diagram of a computer device 2000 illustrating an exemplary embodiment of this application is shown. This computer device can be implemented as a server as described above in this application. The computer device 2000 includes a Central Processing Unit (CPU) 2001, a system memory 2004 including Random Access Memory (RAM) 2002 and Read-Only Memory (ROM) 2003, and a system bus 2005 connecting the system memory 2004 and the CPU 2001. The computer device 2000 also includes a mass storage device 2006 for storing an operating system 2009, application programs 2010, and other program modules 2011.

[0217] The mass storage device 2006 is connected to the central processing unit 2001 via a mass storage controller (not shown) connected to the system bus 2005. The mass storage device 2006 and its associated computer-readable media provide non-volatile storage for the computer device 2000. That is, the mass storage device 2006 may include computer-readable media (not shown), such as a hard disk or a compact disc read-only memory (CD-ROM) drive.

[0218] Without loss of generality, the computer-readable medium may include computer storage media and communication media. Computer storage media include volatile and non-volatile, removable and non-removable media implemented using any method or technology for storing information such as computer-readable instructions, data structures, program modules, or other data. Computer storage media include RAM, ROM, erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other solid-state storage technologies, CD-ROM, digital versatile disc (DVD) or other optical storage, magnetic tape cassettes, magnetic tape, disk storage, or other magnetic storage devices. Of course, those skilled in the art will recognize that the computer storage medium is not limited to the above-mentioned types. The system memory 2004 and mass storage device 2006 described above can be collectively referred to as memory.

[0219] According to various embodiments of this disclosure, the computer device 2000 can also be connected to a remote computer on a network, such as the Internet. That is, the computer device 2000 can be connected to a network 2008 via a network interface unit 2007 connected to the system bus 2005, or the network interface unit 2007 can be used to connect to other types of networks or remote computer systems (not shown).

[0220] The memory also includes at least one computer program stored in the memory, and the central processing unit 2001 executes the at least one computer program to implement all or part of the steps in the methods shown in the above embodiments.

[0221] In an exemplary embodiment, a chip is also provided, the chip including programmable logic circuitry and / or program instructions, which, when the chip is run on a computer device, are used to implement the variable depth cable data processing method described above.

[0222] In an exemplary embodiment, a computer program product is also provided, comprising 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 to implement the variable depth cable data processing method provided in the above-described method embodiments.

[0223] In an exemplary embodiment, a computer-readable storage medium is also provided, which stores a computer program that is loaded and executed by a processor to implement the variable depth cable data processing method provided in the above-described method embodiments.

[0224] 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.

[0225] Those skilled in the art will recognize that the functions described in the embodiments of this application in one or more of the above examples can be implemented using hardware, software, firmware, or any combination thereof. When implemented using software, these functions can be stored in a computer-readable medium or transmitted as one or more instructions or code on a computer-readable medium. Computer-readable media include computer storage media and communication media, wherein communication media include any medium that facilitates the transfer of a computer program from one place to another. Storage media can be any available medium that can be accessed by a general-purpose or special-purpose computer.

[0226] The above description is merely an optional embodiment of this application and is not intended to limit this application. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the protection scope of this application.

Claims

1. A method for processing variable depth cable data, characterized in that, The method includes: Seismic wave data acquired using variable depth cable technology is obtained; the seismic wave data includes the two-way travel time of multiple seismic waves, the offset of the multiple seismic waves, the path of the multiple seismic waves, and the direction of the multiple seismic waves; the seismic waves include incident waves and reflected waves; the incident waves are seismic waves between the shot point and the reflection point, and the reflected waves are seismic waves between the reflection point and the receiver point, and the reflected waves are divided into up-going waves and down-going waves; Based on the offset distance, path, and direction of the multiple seismic waves, a correction wave is set for the multiple seismic waves; the correction wave is axially symmetric to the incident wave about the normal to the reflection point. In response to the reflected wave being an upward wave, the reference plane correction formula for determining the two-way travel time is as follows: in, The depth of the detector point. For the speed of seawater, X This is the offset distance. The propagation time of the two-way reflected wave. The depositional velocity of the strata above the reflection point; In response to the reflected wave being a downflow wave, the reference plane correction formula for determining the two-way travel time is as follows: ； Based on the reference plane correction formula, the time difference correction value for the round-trip travel time corresponding to the first offset distance is calculated; the first offset distance is any one of the plurality of offset distances; The time difference correction value is applied to the round-trip travel time corresponding to the first offset using Sinc interpolation to correct the round-trip travel time corresponding to the first offset.

2. The method according to claim 1, characterized in that, The method further includes: Based on the corrected two-way travel time of the multiple seismic waves, corrected seismic wave data is generated. The corrected seismic wave data is subjected to correlation processing to generate target seismic wave data; the target seismic wave data includes the target two-way travel time of multiple seismic waves, the target offset of the multiple seismic waves, the target path of the multiple seismic waves, and the target direction of the multiple seismic waves. Based on the reference plane correction formula, the target offset corresponding to the target round-trip travel time is inversely corrected.

3. The method according to claim 2, characterized in that, The step of inversely correcting the target offset corresponding to the target round-trip travel time based on the reference plane correction formula includes: Based on the reference plane correction formula, the negative of the time difference correction value of the round-trip travel time of the target corresponding to the first target offset is calculated; the first target offset is any one of the plurality of target offsets; The target round-trip travel time corresponding to the first target offset is inversely corrected by applying the negative number to the target round-trip travel time corresponding to the first target offset using Sinc interpolation.

4. A data processing device for a variable depth cable, characterized in that, The device includes: The acquisition module is used to acquire seismic wave data collected by variable depth cable technology; the seismic wave data includes the two-way travel time of multiple seismic waves, the offset of the multiple seismic waves, the path of the multiple seismic waves, and the direction of the multiple seismic waves; the seismic waves include incident waves and reflected waves; the incident wave is the seismic wave between the shot point and the reflection point, and the reflected wave is the seismic wave between the reflection point and the receiver point, and the reflected wave is divided into up-going waves and down-going waves; The setting module sets correction waves for the multiple seismic waves based on their offset distances, paths, and directions; the correction waves are axially symmetric to the incident waves about the normal to the reflection point. The determination module, in response to the reflected wave being an upward wave, determines the reference plane correction formula for the two-way travel time as follows: in, The depth of the detector point. For the speed of seawater, X This is the offset distance. The propagation time of the two-way reflected wave. The depositional velocity of the strata above the reflection point; In response to the reflected wave being a downflow wave, the reference plane correction formula for determining the two-way travel time is as follows: ； The correction module is used to calculate the time difference correction value of the round-trip travel time corresponding to the first offset distance based on the reference plane correction formula; the first offset distance is any one of the plurality of offset distances; the time difference correction value is applied to the round-trip travel time corresponding to the first offset distance by Sinc interpolation to correct the round-trip travel time corresponding to the first offset distance.

5. A computer device, characterized in that, The computer device includes a processor and a memory, the memory storing at least one computer instruction, which is loaded and executed by the processor to implement the variable depth cable data processing method as described in any one of claims 1 to 3.

6. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one computer instruction, which is loaded and executed by a processor to implement the variable depth cable data processing method as described in any one of claims 1 to 3.

7. A computer program product, characterized in that, The computer program product includes computer instructions stored in a computer-readable storage medium; the computer instructions are read and executed by a processor of a computer device to implement the variable depth cable data processing method as described in any one of claims 1 to 3.

Images