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

Abstract

The invention relates to a variable-depth cable data processing method, device and equipment, a medium and a program product, and relates to the technical field of oil-gas exploration. The method comprises the following steps: acquiring seismic wave data acquired through a variable depth cable technology; setting correction waves of the plurality of seismic waves based on the offset distances of the plurality of seismic waves, the paths of the plurality of seismic waves and the directions of the plurality of seismic waves; determining a calculation formula of a difference value between the travel time of the correction wave and the travel time of the reflection wave as a reference surface correction formula of the two-way travel time; and correcting the two-way travel time of the plurality of seismic waves based on a reference surface correction formula. According to the method, the seismic wave data acquired by the variable-depth cable technology is analyzed to obtain the reference surface correction formula of the ocean variable-depth cable seismic data, so that the shot point and the detection point are positioned on the same reference surface, and the result precision and the imaging quality of conventional processing of the variable-depth cable seismic data are improved.

Description

Technical Field

[0001] The present application relates to the technical field of oil and gas exploration, and in particular to a variable depth cable data processing method, device, equipment, medium and program product. Background Art

[0002] Seismic exploration is an important means of oil and gas exploration, especially offshore oil and gas exploration. Setting a reasonable source and towline sinking depth can ensure that the collected data has good quality.

[0003] In the relevant technology, variable depth cable acquisition can greatly improve the resolution of acquired data compared with conventional flat cable acquisition, and effectively improve the quality of seismic data and the accuracy of seismic interpretation.

[0004] However, since the placement depth of the detector of the variable depth cable acquisition technology will change with the change of the shot offset, the data shot and the detection point are not coplanar, and there is a depth difference between the source and the receiving node, which makes the propagation path of the seismic ray no longer symmetrical, directly affecting the velocity analysis accuracy in the conventional processing of the variable depth cable data and the final imaging quality. Summary of the invention

[0005] The present application provides a variable depth cable data processing method, device, equipment, medium and program product, which can correct the 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 content of the technical solution is as follows.

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

[0007] Acquire seismic wave data collected by variable depth cable technology; the seismic wave data includes round-trip travel time of multiple seismic waves, offset distances of the multiple seismic waves, paths of the multiple seismic waves, and directions of the multiple seismic waves; the seismic waves include incident waves and reflected waves; the incident waves are seismic waves between a shot point and a reflection point, the reflected waves are seismic waves between the reflection point and a detection point, and the reflected waves are divided into upgoing waves and downgoing waves;

[0008] Acquire a reference correction formula of the two-way travel time based on the offset distances of the multiple seismic waves, the paths of the multiple seismic waves, and the directions of the multiple seismic waves;

[0009] Based on the reference surface correction formula, the round-trip travel times of the plurality of seismic waves are corrected.

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

[0011] A data acquisition module, used to acquire seismic wave data collected by the variable depth cable technology; the seismic wave data includes round-trip travel times of multiple seismic waves, offset distances of the multiple seismic waves, paths of the multiple seismic waves, and directions of the multiple seismic waves; the seismic waves include incident waves and reflected waves; the incident waves are seismic waves between a shot point and a reflection point, the reflected waves are seismic waves between the reflection point and a detection point, and the reflected waves are divided into upgoing waves and downgoing waves;

[0012] A setting module, which sets correction waves of the multiple seismic waves based on the offset distances of the multiple seismic waves, the paths of the multiple seismic waves, and the directions of the multiple seismic waves; the correction waves are axially symmetric with the incident waves about the normal line where the reflection point is located;

[0013] A determination module, used to determine the calculation formula of the difference between the travel time of the correction wave and the travel time of the reflected wave as the reference plane correction formula of the round-trip travel time;

[0014] A 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 to:

[0016] In response to the reflected wave being an upgoing wave, the reference plane correction formula is determined as:

[0017]

[0018] Where Δh is the depth of the detection point, V w is the seawater velocity, X is the offset, T 0 (X) is the two-way reflected wave propagation time, V rms is the stacking velocity of the strata above the reflection point.

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

[0020] In response to the reflected wave being a downlink wave, the reference plane correction formula is determined as:

[0021]

[0022] Where Δh is the depth of the detection point, V w is the seawater velocity, X is the offset, T 0 (X) is the two-way reflected wave propagation time, V rms is the stacking velocity of the strata above the reflection point.

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

[0024] Based on the reference surface correction formula, calculating a time difference correction value of a round trip travel time corresponding to a first offset; the first offset is any one of the multiple offsets;

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

[0026] In some embodiments, the device further comprises: an inverse correction module, configured to:

[0027] generating corrected seismic wave data based on the corrected two-way travel times of the plurality of seismic waves;

[0028] Correlation processing is performed on the corrected seismic wave data to generate target seismic wave data; the target seismic wave data includes target two-way travel times of multiple seismic waves, target offsets of the multiple seismic waves, target paths of the multiple seismic waves, and target directions of the multiple seismic waves;

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

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

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

[0032] The inverse number is applied to the target round trip travel time corresponding to the first target offset by means of Sinc interpolation to inversely correct the target round trip travel time corresponding to the first target offset.

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

[0034] According to another aspect of the present application, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores at least one instruction, at least one program, a code set or an instruction set, and 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 the present application, a computer program product is provided, the computer program product comprising computer instructions, the computer instructions being stored in a computer-readable storage medium, and the processor reading and executing the computer instructions from the computer-readable storage medium to implement the variable depth cable data processing method as described above.

[0036] The technical solution provided by the embodiments of the present application may have the following beneficial effects:

[0037] The seismic wave data collected by the variable depth cable technology is obtained, and correction waves of multiple seismic waves are set. The calculation formula of the difference between the travel time of the correction wave and the travel time of the reflection wave is determined as the datum correction formula of the two-way travel time. Based on the datum correction formula, the two-way travel time of multiple seismic waves is corrected. Since the shot points and the detection points of the seismic wave data collected by the variable depth cable technology are not coplanar, the processing accuracy and imaging quality are affected. This scheme analyzes the seismic wave data collected by the variable depth cable technology and obtains the datum correction formula of the marine variable depth cable seismic data datum, so that the shot points and the detection points are located on the same datum, thereby improving the accuracy and imaging quality of the conventional processing of variable depth cable seismic data. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0039] Figure 1 is a schematic diagram of a conventional offshore flat cable technology provided by an exemplary embodiment of the present application;

[0040] Figure 2 is a schematic diagram of an offshore variable depth cable technology provided by an exemplary embodiment of the present application;

[0041] Figure 3 is a flow chart of a variable depth cable data processing method provided by an exemplary embodiment of the present application;

[0042] Figure 4 is a schematic diagram of the geometric relationship of the upgoing wave reference plane correction provided by an exemplary embodiment of the present application;

[0043] Figure 5 It is a schematic diagram of the geometric relationship of the downlink wave reference plane correction provided by an exemplary embodiment of the present application;

[0044] Figure 6 is a flow chart of a variable depth cable data processing method provided by another exemplary embodiment of the present application;

[0045] Figure 7 is a flow chart of a variable depth cable data processing method provided by another exemplary embodiment of the present application;

[0046] Figure 8 is a flow chart of a variable depth cable data processing method provided by yet another exemplary embodiment of the present application;

[0047] Fig. 9 It is a schematic diagram of seismic data of an ocean variable depth cable up-travel wave model provided by an exemplary embodiment of the present application;

[0048] Fig.10 is a schematic diagram of an upgoing wave reference surface correction of variable depth cable model data provided by an exemplary embodiment of the present application;

[0049] Fig.11 It is a schematic diagram of upgoing wave reference plane correction + reverse reference plane correction of variable depth cable model data provided by an exemplary embodiment of the present application;

[0050] Fig.12 This is a schematic diagram of seismic data of a down-going wave model of an ocean variable depth cable provided by an exemplary embodiment of the present application;

[0051] Fig.13 is a schematic diagram of downlink wave reference surface correction of variable depth cable model data provided by an exemplary embodiment of the present application;

[0052] Fig.14 It is a schematic diagram of downlink wave reference plane correction + reverse reference plane correction of variable depth cable model data provided by an exemplary embodiment of the present application;

[0053] Fig.15 It is a schematic diagram of upgoing wave reference surface correction + conventional dynamic correction of variable depth cable model data provided by an exemplary embodiment of the present application;

[0054] Fig.16 It is a schematic diagram of downlink wave reference surface correction + conventional dynamic correction of variable depth cable model data provided by an exemplary embodiment of the present application;

[0055] Fig.17 It is a schematic diagram of traveling wave model data + conventional dynamic correction on a variable depth cable provided by an exemplary embodiment of the present application;

[0056] Fig.18 It is a schematic diagram of variable depth cable downlink wave model data + conventional dynamic correction provided by an exemplary embodiment of the present application;

[0057] Fig.19 is a block diagram of a data processing device for a variable depth cable shown in an exemplary embodiment of the present application;

[0058] Fig. 20 It is a structural block diagram of a computer device provided by an exemplary embodiment of the present application.

[0059] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application. DETAILED DESCRIPTION

[0060] In order to make the objectives, technical solutions and advantages of the present application clearer, the implementation methods of the present application will be further described in detail below with reference to the accompanying drawings.

[0061] Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The implementations described in the following exemplary embodiments do not represent all implementations consistent with the present application. Instead, they are merely examples of devices and methods consistent with some aspects of the present application as detailed in the appended claims.

[0062] The terms used in this disclosure are for the purpose of describing specific embodiments only and are not intended to limit the disclosure. The singular forms of "a", "said" and "the" used in this disclosure and the appended claims are also intended to include plural forms unless the context clearly indicates otherwise. It should also be understood that the term "and / or" used herein refers to and includes any or all possible combinations of one or more 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 relevant data must comply with the relevant laws, regulations and standards of the relevant countries and regions. For example, the object behaviors such as attack operations involved in this application are all obtained with full authorization.

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

[0065] The following are some definitions of terms involved in this application:

[0066] Conventional dynamic correction: also known as normal move out correction (NMO). For multiple coverage seismic records, horizontal stacking is performed on the common depth point data set. Due to the existence of normal move out of non-zero shot offset, the time-distance curve of the reflection wave at the common depth point is a hyperbola. NMO is to correct the arrival time of the reflection waves from the same interface and the same point on each channel with different shot offsets to the echo time at the common center point after normal move out correction, so as to ensure that they can be superimposed in phase during superposition, forming a superposition channel with prominent reflection wave energy (equivalent to a self-excited and self-received recording channel).

[0067] Velocity parameters are required in NMO processing. For horizontal layered media, if the velocity is selected correctly, the reflection hyperbola can be corrected to a straight line, and each trace can be stacked in phase during stacking. If the velocity used is too large, the correction will be insufficient; conversely, if the velocity used is too small, the correction will be excessive. In both cases, it is impossible to guarantee the same-phase stacking during stacking. For single-cover records, dynamic correction can be used for shot gather records to directly obtain single-cover seismic profiles.

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

[0069] The Sinc interpolation method has good frequency domain characteristics and theorem preservation characteristics, but this requires certain truncation operations on the function during interpolation. At the same time, the Sinc interpolation method has a high computational complexity and is forced to interpolate within a certain frequency band. If the distance between the interpolation point and the sampling point is too large, the interpolation error will increase. Therefore, in practical applications, it is necessary to optimize the interpolation parameters to ensure that the interpolation error is within an acceptable range. The Sinc interpolation method has been widely used in signal processing, image processing, seismic exploration and other fields. In this scheme, since the correction amount obtained by time difference correction is often not an integer sample point, it is necessary to interpolate the seismic trace, and the use of Sinc interpolation can obtain very good results.

[0070] Snell's law: In optics, Snell's law is a formula that describes the relationship between the angle of incidence and the angle of refraction when light or other waves enter from one medium into another. This law states that the ratio of the sine of the angle of incidence to the sine of the angle of refraction is a certain value, and this certain value is related to the incident and refractive media.

[0071] At present, offshore variable depth cable data acquisition technology has become increasingly mature. Acquisition practice shows that variable depth cable acquisition can greatly 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. However, the key issue currently faced is how to correctly process the acquired data.

[0072] Please refer to Figure 1 , which shows a schematic diagram of a conventional offshore flat cable technology provided by an exemplary embodiment of the present application; please refer to Figure 2 , which shows a schematic diagram of an offshore variable depth cable technology provided by an exemplary embodiment of the present application. Figure 1 As shown in Figure 2, the shot points and receiver points of conventional offshore flat cable technology are almost coplanar, and the propagation paths of seismic rays are almost symmetrical; Figure 2 As shown in the figure, the shot points and receiver points of the offshore variable depth cable technology are not coplanar, and the propagation path of the seismic rays is no longer symmetrical.

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

[0074] Since the immersion depth of the detector of the variable depth cable acquisition technology will change with the change of the shot offset, the data shot and the detection point are not coplanar, and there is a depth difference between the source and the receiving node, the propagation path of the seismic ray is no longer symmetrical. Many theoretical assumptions in the conventional processing flow are not valid or the algorithm is difficult to implement efficiently, which directly affects multiple wave suppression, velocity analysis, broadband processing, dynamic correction and stacking quality control in the variable depth cable data processing.

[0075] In view of the above problems in variable depth cable data processing, processing personnel urgently need a reference correction technology suitable for marine variable depth cable data to correct the source and receiving point to a unified reference plane, thereby improving the resolution of the results of conventional processing of variable depth cable seismic data and the final imaging quality.

[0076] Among them, after correction to a unified reference plane, the shot points and the detection points are located on the same horizontal plane, which is equivalent to ordinary flat cable data, so that the conventional processing module can also process variable depth cable data well. On the contrary, if the reference plane correction is not performed and the conventional module is used directly for processing, the accuracy of the processing results will be reduced because many theoretical assumptions in the conventional processing algorithm are not valid.

[0077] Please refer to Figure 3 , which shows a flow chart of a variable depth cable data processing method provided by an exemplary embodiment of the present application. The method is executed by a computer device, such as Figure 3As shown, the method may include step 310 , step 320 , step 330 , and step 340 .

[0078] Step 310: Acquire seismic wave data collected by variable depth cable technology; the seismic wave data include round-trip travel times of multiple seismic waves, offset distances of multiple seismic waves, paths of multiple seismic waves, and directions of multiple seismic waves; seismic waves include incident waves and reflected waves; incident waves are seismic waves between shot points and reflection points, reflected waves are seismic waves between reflection points and detection points, and reflected waves are divided into upgoing waves and downgoing waves.

[0079] In the embodiment of the present application, the computer device can obtain Figure 1 Seismic wave data collected using variable depth cable technology is shown.

[0080] The seismic wave data may be a collection of multiple acquisition data from multiple variable depth cables, including the round trip travel time of multiple seismic waves, the offset distance of multiple seismic waves, the paths of multiple seismic waves, and the directions of multiple seismic waves. Based on the round trip travel time of multiple seismic waves and the offset distance of multiple seismic waves, the following can be constructed: Fig. 9 or Fig.12 The forward shot gather record is shown.

[0081] Among them, 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 detection point; according to the direction of the reflected wave, the reflected wave is divided into upgoing waves and downgoing waves.

[0082] Step 320: Based on the offset distances, paths and directions of the multiple seismic waves, correction waves of the multiple seismic waves are set; the correction waves are axially symmetric with the incident waves about the normal line where the reflection point is located.

[0083] In an embodiment of the present application, the computer device can set correction waves for the multiple seismic waves according to the offset distances of the multiple seismic waves, the paths of the multiple seismic waves, and the directions of the multiple seismic waves.

[0084] Among them, the correction wave is set on one side of the reflected wave, and the correction wave and the incident wave are axially symmetric about the normal line of the reflection point, so that the shot point and the detection point of the variable depth cable technology are located on the same reference plane to eliminate the influence of the depth difference between the shot point and the detection point.

[0085] Step 330: Determine the calculation formula for the difference between the travel time of the correction wave and the travel time of the reflected wave as the reference plane correction formula for the round-trip travel time.

[0086] In an embodiment of the present application, the computer device can analyze the geometric relationship between the correction wave path and the reflected wave path, derive a calculation formula for the difference between the travel time of the correction wave and the travel time of the reflected wave, and determine the calculation formula as the reference plane correction formula for the round-trip travel time.

[0087] The above-mentioned reference plane correction formula may include an upgoing wave reference plane correction formula and a downgoing wave reference plane correction formula.

[0088] Step 340: Correct the two-way travel times of multiple seismic waves based on the reference correction formula.

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

[0090] The time difference correction is performed according to the time difference correction value obtained from the reference plane correction formula. The point-by-point position shape of the seismic wave data collected by the variable depth cable technology can be adjusted accordingly, so that the seismic wave data collected by the variable depth cable technology is basically consistent with the data collected by the flat cable, thereby realizing that the shot point and the detection point are located on the same reference plane.

[0091] After correction to a unified reference plane, the shot points and the detection points are located on the same horizontal plane, which is equivalent to the seismic wave data collected by ordinary flat cables, so that the conventional processing module can also better process the seismic wave data collected by variable depth cables. On the contrary, if the reference plane correction is not performed and the conventional module is directly used for processing, the accuracy of the processing results will be reduced because many theoretical assumptions in the conventional processing algorithm are not valid.

[0092] In summary, the scheme shown in the embodiment of the present application obtains seismic wave data collected by the variable depth cable technology, sets correction waves for multiple seismic waves, determines the calculation formula of the difference between the travel time of the correction wave and the travel time of the reflection wave as the reference plane correction formula for the two-way travel time, and corrects the two-way travel time of multiple seismic waves based on the reference plane correction formula. Since the shot points and the detection points of the seismic wave data collected by the variable depth cable technology are not coplanar, the processing accuracy and imaging quality are affected; this scheme analyzes the seismic wave data collected by the variable depth cable technology to derive the reference plane correction formula for the reference plane of the marine variable depth cable seismic data, so that the shot points and the detection points are located on the same reference plane, thereby improving the accuracy and imaging quality of the conventional processing of variable depth cable seismic data.

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

[0094] In response to the reflected wave being an upgoing wave, the reference plane correction formula is determined as:

[0095]

[0096] Where Δh is the depth of the detection point, V w is the seawater velocity, X is the offset, T 0 (X) is the two-way reflected wave propagation time, V rms is the stacking velocity of the strata above the reflection point.

[0097] Please refer to Figure 4 , which shows a schematic diagram of the geometric relationship of the upgoing wave reference correction provided by an exemplary embodiment of the present application. According to the upgoing wave propagation ray path and direction information of the variable depth cable, an accurate upgoing wave reference correction formula for the marine variable depth cable seismic data is derived:

[0098] In the embodiment of the present application, for data with zero offset, due to self-excitation and self-reception, it is obvious that the calculation formula for the difference between the travel time of the correction wave and the travel time of the reflection wave is:

[0099]

[0100] Where Δh is the depth of the detection point, V w is the sea water speed.

[0101] like Figure 4 As shown, for data with non-zero offset, the angle between the correction wave and the normal line of the detection point is θ 1 , the calculation formula for the difference between the travel time of the correction wave and the travel time of the reflected wave is approximated as:

[0102]

[0103] Among them, this is based on a certain degree of approximate assumption. When actually collecting data at sea, since the seabed depth is much greater than the cable depth, the time difference between the correction wave and the reflected wave can be approximately calculated.

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

[0105]

[0106] Where H is half offset, which is half of the offset X, T ref (H) is the travel time of the seismic reflection wave when the offset is H, which is the one-way travel time and the two-way reflection wave propagation time T 0 Half of (X), V rms is the stacking velocity of the strata above the reflection point.

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

[0108]

[0109] Formula (D1) is the derived datum correction formula for upgoing ocean variable depth cable seismic data.

[0110] An embodiment of the present 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 variable depth cable seismic upgoing wave, a reference plane correction formula for marine variable depth cable seismic data is derived to correct the source and the receiving point to a unified reference plane.

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

[0112] In response to the reflected wave being a downlink wave, the reference plane correction formula is determined as:

[0113]

[0114] Where Δh is the depth of the detection point, V w is the seawater velocity, X is the offset, T 0 (X) is the two-way reflected wave propagation time, V rms is the stacking velocity of the strata above the reflection point.

[0115] Please refer to Figure 5 , which shows a schematic diagram of the geometric relationship of the downgoing wave reference correction provided by an exemplary embodiment of the present application. According to the path and direction information of the downgoing wave propagation ray of the variable depth cable, an accurate downgoing wave reference correction formula for the marine variable depth cable seismic data is derived:

[0116] In the embodiment of the present application, for data with zero offset, due to self-excitation and self-reception, it is obvious that the calculation formula for the difference between the travel time of the correction wave and the travel time of the reflection wave is:

[0117]

[0118] Where Δh is the depth of the detection point, V w is the sea water speed.

[0119] like Figure 5 As shown in the figure, for data with non-zero offset, the angle θ between the correction wave and the normal line of the detection point is obtained. 2 When , the calculation formula for the difference between the travel time of the correction wave and the travel time of the reflected wave is approximated as:

[0120]

[0121] Among them, this is based on a certain degree of approximate assumption. When actually collecting data at sea, since the seabed depth is much greater than the cable depth, the time difference between the correction wave and the reflected wave can be approximately calculated.

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

[0123]

[0124] Where H is half offset, which is half of the offset X, T ref (H) is the travel time of the seismic reflection wave when the offset is H, which is the one-way travel time and the two-way reflection wave propagation time T 0 Half of (X), V rms is the stacking velocity of the strata above the reflection point.

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

[0126]

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

[0128] An embodiment of the present application provides a reference plane correction formula when the reflected wave is a downgoing wave. By analyzing the geometric relationship of the propagation ray path of the variable depth cable seismic downgoing wave, a reference plane correction formula for marine variable depth cable seismic data is derived to correct the source and receiving point to a unified reference plane.

[0129] Please refer to Figure 6 , which shows a flow chart of a variable depth cable data processing method provided by another exemplary embodiment of the present application. The method is executed by a computer device, such as Figure 6 As shown above Figure 3 Step 340 in the illustrated embodiment may be implemented as step 340a and step 340b.

[0130] Step 340a: Calculate the time difference correction value of the round-trip travel time corresponding to the first offset based on the reference plane correction formula; the first offset is any one of the multiple offsets.

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

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

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

[0134] In the embodiment of the present application, after the time difference correction value of the round-trip travel time corresponding to the first offset is calculated in step 340a, 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] Then, based on step 340a and step 340b, the round trip travel time corresponding to the second offset is corrected.

[0136] By analogy, the time difference correction value is calculated and applied one by one for the seismic wave data collected by the variable depth cable technology, so that the shot point and the detection point of the seismic wave data are located on the same reference plane, completing the reference plane correction of the variable depth cable data.

[0137] The embodiment of the present application provides a technical solution for correcting the two-way travel time of multiple seismic waves based on a datum correction formula, and performs time difference correction on multiple two-way travel times with the same offset parameters in turn to improve the efficiency of datum correction.

[0138] Please refer to Figure 7 , which shows a flow chart of a variable depth cable data processing method provided by another exemplary embodiment of the present 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 step 350 , step 360 , and step 370 .

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

[0140] In an embodiment of the present application, after step 340 corrects the two-way travel time of multiple seismic waves based on the reference plane correction formula, the computer device can generate corrected seismic wave data so that the corrected seismic wave data can be subjected to multiple wave suppression, velocity analysis, broadband processing, dynamic correction, and stacking quality control using theoretical assumptions in the conventional processing flow.

[0141] Step 360: performing correlation processing on the corrected seismic wave data to generate target seismic wave data; the target seismic wave data includes target two-way travel times of multiple seismic waves, target offset distances of multiple seismic waves, target paths of multiple seismic waves, and target directions of multiple seismic waves.

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

[0143] The target seismic wave data may be similar to the seismic wave data, including target two-way travel times of multiple seismic waves, target offsets of multiple seismic waves, target paths of multiple seismic waves, and target directions of multiple seismic waves.

[0144] Step 370: Based on the reference surface correction formula, inversely correct the target round-trip travel time corresponding to the target offset.

[0145] In the embodiment of the present application, according to the reference plane correction formula obtained in step 330, the computer device can perform inverse correction on the target seismic wave data.

[0146] Among them, since some processing modules can only be completed on data of a unified reference surface, such as a module for removing ocean multiple wave 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 wave interference, remove the interference of the corrected seismic wave data on a unified reference surface; then perform inverse correction. Although the seismic wave data is similar in form to the original seismic wave data, the seismic wave data that has been processed by the module and then inversely corrected has had interference removed, which is more conducive to improving the signal-to-noise ratio of subsequent processing modules.

[0147] The embodiment of the present application provides a technical solution for inverse correction based on the reference plane correction formula. Through the inverse reference plane correction, the reference plane correction result can be restored to the original variable depth cable seismic wave data, thereby realizing reversible correction processing and eliminating interference with the original variable depth cable seismic wave data.

[0148] Please refer to Figure 8 , which shows a flow chart of a variable depth cable data processing method provided by another exemplary embodiment of the present application. The method is executed by a computer device, such as Figure 8 As shown above Figure 7 Step 370 in the illustrated embodiment may be implemented as step 370a and step 370b.

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

[0150] In the embodiment of the present application, after step 360 performs correlation processing on the corrected seismic wave data to generate target seismic wave data, the computer device can use the reference plane correction formula obtained in step 330 to refer to the above Figure 6In the illustrated embodiment, the inverse of the time difference correction value for the target round trip travel time corresponding to the first target offset is calculated.

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

[0152] In the embodiment of the present application, after step 370a calculates the inverse 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 inverse 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] Then, based on step 370a and step 370b, the target round trip travel time corresponding to the second target offset is corrected.

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

[0155] The embodiment of the present application provides a supplementary technical solution for inverse correction based on the reference plane correction formula, and performs time difference inverse correction on multiple target round-trip travel times with the same target offset in sequence, so as to improve the efficiency of inverse reference plane correction.

[0156] In summary, this application aims at the problems of low velocity analysis accuracy and final imaging quality in conventional processing of variable depth cable data. By analyzing the geometric relationship of the propagation ray paths of variable depth cable seismic up / down waves, an accurate marine variable depth cable seismic data datum correction formula is derived, which can distinguish variable depth cable up / down waves and accurately correct the source and receiving point to a unified datum, thereby improving the resolution of the results of variable depth cable seismic data processing and the final imaging quality. Through the reverse datum correction, the datum correction result can be restored to the original input variable depth cable data, realizing basically reversible correction processing.

[0157] Based on the method shown in the above embodiment of the present application, the embodiment of the present application provides a schematic implementation flow of a variable depth cable data processing method, including the following steps.

[0158] S1: Fig. 9 As shown in the figure, the ocean variable depth cable up-going wave seismic data S(t i ,x j ).

[0159] Wherein, t represents the round-trip travel time of the reflected wave, x represents the offset distance, i represents the serial number of the sampling point along the time direction, j represents the serial 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.

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

[0161] S2: For the upgoing variable depth cable seismic data input in S1 above, the time difference correction value of the two-way travel time is obtained, and the two-way travel time of the above data is corrected for the time difference so that the shot point and the detection point are located on the same reference plane.

[0162] S21: Select the upgoing wave variable depth cable seismic data S(t i ,x j ), that is, the number of sampling points increases, while the offset parameter remains unchanged, for example, S(t i ,x 1 ), i = 1, 2…6000, j = 1, use the upgoing wave reference correction formula:

[0163]

[0164]

[0165] For S(t i ,x 1 ) The time difference correction value TX(t i )=Δt(θ), i=1, 2...6000.

[0166] In the embodiment of the present application, S(t i ,x 1 ) The depth of the variable depth cable is ΔH = 20 meters, and the seawater velocity V w =1500m / s, offset distance is X=25m, select the first sample point T of this track 0 (X) = 2ms, the superposition velocity V at this point rmsThe time difference correction value TX(t 1 ).

[0167] Similarly, the time difference correction values ​​TX(t 2 )…Txt 6000 ), that is, S(t i ,x 1 ) The time difference correction value TX(t i ), i=1, 2…6000.

[0168] S22: Use Sinc interpolation to apply the time difference correction value TX (t i ) to make time difference correction so that the gun check points are located on the same reference plane.

[0169] S23: Repeat steps S21 and S22 to select the upgoing wave variable depth cable seismic data S(t i ,x j ), such as S(t i ,x 2 ), i = 1, 2 ... 6000, j = 2, and the time difference correction value TX (t i ) and apply it so that the gun check point is located on the same reference plane; repeat the operation 120 times, so that S(t i ,x i ) All 120 records in the data have been corrected.

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

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

[0172] S3: Fig.12 As shown in the figure, the variable depth cable downlink seismic data S(T I ,X J ), I=1, 2...N, J=1, 2...M.

[0173] T represents the round-trip travel time of the reflected wave, X represents the offset, I represents the serial number of the sampling point along the time direction; J represents the serial number of the sampling point along the offset 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 direction.

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

[0175] S4: For the downlink variable depth cable seismic data input in S3 above, the time difference correction value of the two-way travel time is obtained, and the two-way travel time of the above data is corrected for the time difference so that the shot point and the detection point are located on the same reference plane.

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

[0177]

[0178]

[0179] For S(T I ,X 1 ) The time difference correction value TX(t I )=Δt(θ), I=1, 2...6000.

[0180] In the present embodiment, S(T I ,X 1 ) The depth of the variable depth cable is ΔH = 20 meters, and the seawater velocity V w =1500m / s, offset distance is X=25m, select the first sample point T of this track 0 (X) = 2ms, the superposition velocity V at this point rmsThe time difference correction value TX(t 1 ).

[0181] Similarly, the time difference correction values ​​TX(t 2 )…Txt 6000 ), that is, the time difference correction value TX(t I ), I = 1, 2…6000.

[0182] S42: Using Sinc interpolation, apply the time difference correction value TX (t I ) to make time difference correction so that the gun check points are located on the same reference plane.

[0183] S43: Repeat steps S41 and S42 to select the downlink wave variable depth cable seismic data S(T I ,X 1 ), such as S(T I ,X 2 ), I = 1, 2 ... 6000, J = 2, and the time difference correction value TX (t I ) and apply it so that the shot check points of this channel are located on the same reference plane. Repeat the operation M times, so that all M channel records of the data are corrected.

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

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

[0186] In the embodiment of the present application, in order to verify the effect of the upgoing wave reference plane correction and the downgoing wave reference plane correction, conventional dynamic correction is performed on each result data.

[0187] like Fig.15 As shown in Figure 1, after the upgoing wave reference plane correction is performed on the variable depth cable upgoing wave model data, the conventional dynamic correction result is performed; Fig.16 As shown in the figure, after the variable depth cable down-going wave model data is corrected for the down-going wave reference surface, the conventional dynamic correction result is performed; Fig.17As shown in Figure 1, the upgoing wave model data of the variable depth cable was not corrected for the upgoing wave reference surface, and the conventional dynamic correction result was directly performed; Fig.18 As shown in the figure, for the variable depth cable down-going wave model data, no down-going wave reference correction was performed, and the conventional dynamic correction result was directly performed.

[0188] contrast Fig.15 , Fig.16 , Fig.17 and Fig.18 It can be seen that due to the influence of the depth difference between the shot point and the detection point, directly performing conventional dynamic correction cannot level the upgoing wave channel set or the downgoing wave channel set of the original variable depth cable model data; however, after the upgoing wave reference plane correction processing and the downgoing wave reference plane correction processing of the embodiment of the present application, the shot point and the detection point of the variable depth cable data are located on the same reference plane, which effectively eliminates the influence of the depth difference between the shot and detection points. Therefore, conventional dynamic correction can also correctly level the upgoing wave or downgoing wave co-phase axis of the variable depth cable data.

[0189] For the reverse reference surface correction, compare the attached Fig. 9 With attached Fig.11 , Attachment Fig.12 With attached Fig.14 It can be seen that the variable depth cable data after the input upgoing wave reference plane correction and the variable depth cable data after the downgoing wave reference plane correction can be restored to the original input variable depth cable data after the upgoing wave anti-reference plane correction and the downgoing wave anti-reference plane correction respectively, so the reference plane correction technology is reversible.

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

[0191] In summary, the present application can accurately implement reversible reference correction of marine variable depth cable seismic data.

[0192] This application introduces a method for datum correction of marine variable depth cable seismic data to address the problem that conventional processing affects processing accuracy and imaging quality due to the non-coplanarity of variable depth cable acquisition seismic shots and checkpoints. This technical solution analyzes the ray path and direction information of variable depth cable acquisition, and derives the datum correction formulas for marine variable depth cable seismic data for up / down waves respectively. Using information such as velocity, offset, water velocity and cable depth, the variable depth cable data can be accurately corrected to sea level so that the shot checkpoints are located on the same datum, thereby improving the accuracy and imaging quality of conventional processing of variable depth cable seismic data. Through reverse datum correction, the datum correction result can be restored to the original input variable depth cable data, achieving basically reversible correction processing.

[0193] Fig.19 It shows a block diagram of a data processing device for a variable depth cable according to an exemplary embodiment of the present application. The device can be used to perform the following steps: 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 Fig.19 As shown, the device comprises:

[0194] The data acquisition module 1901 is used to acquire seismic wave data collected by the variable depth cable technology; the seismic wave data includes the round-trip travel time of multiple seismic waves, the offset distance of multiple seismic waves, the paths of multiple seismic waves, and the directions of multiple seismic waves; seismic waves include incident waves and reflected waves; the incident wave is the seismic wave between the shot point and the reflection point, the reflected wave is the seismic wave between the reflection point and the detection point, and the reflected wave is divided into an upgoing wave and a downgoing wave;

[0195] A setting module 1902 is used to set correction waves of the multiple seismic waves based on the offset distances of the multiple seismic waves, the paths of the multiple seismic waves, and the directions of the multiple seismic waves; the correction waves are axially symmetric with the incident waves about the normal line where the reflection point is located;

[0196] A determination module 1903 is used to determine the calculation formula of the difference between the travel time of the correction wave and the travel time of the reflection wave as the reference plane correction formula of 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 determination module 1903 is used to:

[0199] In response to the reflected wave being an upgoing wave, the reference plane correction formula is determined as:

[0200]

[0201] Where Δh is the depth of the detection point, V w is the seawater velocity, X is the offset, T 0 (X) is the two-way reflected wave propagation time, V rms is the stacking velocity of the strata above the reflection point.

[0202] In some embodiments, the determination module 1903 is used to:

[0203] In response to the reflected wave being a downlink wave, the reference plane correction formula is determined as:

[0204]

[0205] Where Δh is the depth of the detection point, V w is the seawater velocity, X is the offset, T 0 (X) is the two-way reflected wave propagation time, V rms is the stacking velocity of the strata above the reflection point.

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

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

[0208] The time difference correction value is applied to the two-way travel time corresponding to the first offset by means of Sinc interpolation to correct the two-way travel time corresponding to the first offset.

[0209] In some embodiments, the apparatus further comprises: an inverse correction module, configured to:

[0210] Generating corrected seismic wave data based on the corrected two-way travel times of the plurality of seismic waves;

[0211] Correlation processing is performed on the corrected seismic wave data to generate target seismic wave data; the target seismic wave data includes target two-way travel times of multiple seismic waves, target offset distances of multiple seismic waves, target paths of multiple seismic waves, and target directions of multiple seismic waves;

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

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

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

[0215] The target round-trip travel time corresponding to the first target offset is inversely corrected by applying an inverse number to the target round-trip travel time corresponding to the first target offset in a Sinc interpolation manner.

[0216] Fig. 20The structural block diagram of a computer device 2000 shown in an exemplary embodiment of the present application is shown. The computer device can be implemented as the server in the above-mentioned solution of the present application. The computer device 2000 includes a central processing unit (CPU) 2001, a system memory 2004 including a random access memory (RAM) 2002 and a read-only memory (ROM) 2003, and a system bus 2005 connecting the system memory 2004 and the central processing unit 2001. The computer device 2000 also includes a mass storage device 2006 for storing an operating system 2009, an application program 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 medium provide non-volatile storage for the computer device 2000. That is, the mass storage device 2006 may include a computer readable medium (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 by 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), Electronically Erasable Programmable Read-Only Memory (EEPROM) flash memory or other solid-state storage technology, CD-ROM, Digital Versatile Disc (DVD) or other optical storage, cassette, tape, disk storage or other magnetic storage devices. Of course, those skilled in the art will know that the computer storage medium is not limited to the above. The above-mentioned system memory 2004 and mass storage device 2006 can be collectively referred to as memory.

[0219] According to various embodiments of the present disclosure, the computer device 2000 can also be connected to a remote computer on the network through a network such as the Internet. That is, the computer device 2000 can be connected to the network 2008 through the 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, which is stored in the memory. The central processing unit 2001 implements all or part of the steps in the methods shown in the above embodiments by executing the at least one computer program.

[0221] In an exemplary embodiment, a chip is also provided. The chip includes a programmable logic circuit and / or program instructions. When the chip runs on a computer device, it is used to implement the variable depth cable data processing method in the above aspect.

[0222] In an exemplary embodiment, a computer program product is also provided, the computer program product includes computer instructions, the computer instructions are stored in a computer-readable storage medium. A processor of a computer device reads the computer instructions from the computer-readable storage medium, and the processor reads and executes the computer instructions from the computer-readable storage medium to implement the variable depth cable data processing method provided by the above-mentioned method embodiments.

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

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

[0225] Those skilled in the art should be aware that in one or more of the above examples, the functions described in the embodiments of the present application can be implemented with 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 codes on a computer-readable medium. Computer-readable media include computer storage media and communication media, wherein the communication media include any media that facilitates the transmission of a computer program from one place to another. The storage medium can be any available medium that a general or special-purpose computer can access.

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

Claims

1. A variable depth cable data processing method, It is characterized in that The method comprises: Acquire seismic wave data collected by variable depth cable technology; the seismic wave data includes round trip travel times of multiple seismic waves, offset distances of the multiple seismic waves, paths of the multiple seismic waves, and directions of the multiple seismic waves; the seismic waves include incident waves and reflected waves; the incident waves are seismic waves between a shot point and a reflection point, the reflected waves are seismic waves between the reflection point and a detection point, and the reflected waves are divided into upgoing waves and downgoing waves; Based on the offset distances of the multiple seismic waves, the paths of the multiple seismic waves, and the directions of the multiple seismic waves, setting correction waves of the multiple seismic waves; the correction waves are axially symmetric with the incident waves about the normal line where the reflection point is located; Determine the calculation formula of the difference between the travel time of the correction wave and the travel time of the reflected wave as the reference plane correction formula of the round-trip travel time; Based on the reference surface correction formula, the round-trip travel times of the plurality of seismic waves are corrected.

2. The method according to claim 1, It is characterized in that The calculation formula for the difference between the travel time of the correction wave and the travel time of the reflected wave is determined as the reference plane correction formula, including: In response to the reflected wave being an upgoing wave, the reference plane correction formula is determined as: Where Δh is the depth of the detection point, V w is the seawater velocity, X is the offset, T 0 (X) is the two-way reflected wave propagation time, V rms is the stacking velocity of the strata above the reflection point.

3. The method according to claim 1, It is characterized in that The calculation formula for the difference between the travel time of the correction wave and the travel time of the reflected wave is determined as the reference plane correction formula, including: In response to the reflected wave being a downlink wave, the reference plane correction formula is determined as: Where Δh is the depth of the detection point, V w is the seawater velocity, X is the offset, T 0 (X) is the two-way reflected wave propagation time, V rms is the stacking velocity of the strata above the reflection point.

4. The method according to any one of claims 1 to 3, It is characterized in that The step of correcting the two-way travel time of the plurality of seismic waves based on the reference correction formula comprises: Based on the reference surface correction formula, calculating a time difference correction value of a round trip travel time corresponding to a first offset; the first offset is any one of the multiple offsets; The time difference correction value is applied to the round-trip travel time corresponding to the first offset by means of Sinc interpolation to correct the round-trip travel time corresponding to the first offset.

5. The method according to any one of claims 1 to 3, It is characterized in that The method further comprises: generating corrected seismic wave data based on the corrected two-way travel times of the plurality of seismic waves; Correlation processing is performed on the corrected seismic wave data to generate target seismic wave data; the target seismic wave data includes target two-way travel times of multiple seismic waves, target offsets of the multiple seismic waves, target paths of the multiple seismic waves, and target directions of the multiple seismic waves; The target round-trip travel time corresponding to the target offset is inversely corrected based on the reference plane correction formula.

6. The method according to claim 5, It is characterized in that The inverse correction of the target round-trip travel time corresponding to the target offset based on the reference correction formula includes: Based on the reference plane correction formula, calculating the inverse of the time difference correction value of the target round trip travel time corresponding to the first target offset; the first target offset is any one of the multiple target offsets; The inverse number is applied to the target round trip travel time corresponding to the first target offset by means of Sinc interpolation to inversely correct the target round trip travel time corresponding to the first target offset.

7. A data processing device for a variable depth cable, It is characterized in that The device comprises: An acquisition module is used to acquire seismic wave data collected by the variable depth cable technology; the seismic wave data includes round-trip travel times of multiple seismic waves, offset distances of the multiple seismic waves, paths of the multiple seismic waves, and directions of the multiple seismic waves; the seismic waves include incident waves and reflected waves; the incident waves are seismic waves between a shot point and a reflection point, the reflected waves are seismic waves between the reflection point and a detection point, and the reflected waves are divided into upgoing waves and downgoing waves; A setting module, which sets correction waves of the multiple seismic waves based on the offset distances of the multiple seismic waves, the paths of the multiple seismic waves, and the directions of the multiple seismic waves; the correction waves are axially symmetric with the incident waves about the normal line where the reflection point is located; A determination module, used to determine the calculation formula of the difference between the travel time of the correction wave and the travel time of the reflected wave as the reference plane correction formula of the round-trip travel time; A correction module is used to correct the two-way travel time of the multiple seismic waves based on the reference surface correction formula.

8. A computer device, It is characterized in that The computer device comprises a processor and a memory, wherein the memory stores at least one computer instruction, and the at least one computer instruction is loaded and executed by the processor to implement the variable depth cable data processing method according to any one of claims 1 to 6.

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

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

Images