A marine seismic survey line leveling method, processing system and storage medium

By grouping marine seismic survey lines and calculating intersection differences, and using a semi-systematic difference adjustment method, the problem of inconsistent marine seismic survey line data was solved, rapid leveling was achieved, the accuracy of the survey lines and the applicability of the data were improved, and more accurate imaging of underground structures and resource exploration were supported.

CN117706611BActive Publication Date: 2026-06-02GUANGZHOU MARINE GEOLOGICAL SURVEY

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
GUANGZHOU MARINE GEOLOGICAL SURVEY
Filing Date
2023-11-14
Publication Date
2026-06-02

AI Technical Summary

Technical Problem

In existing technologies, marine seismic survey line data suffers from inconsistent reflection times at the intersections of different survey lines due to the influence of tides and waves, resulting in discontinuous stratigraphy, which affects the accuracy of the survey lines and the applicability of the data. Furthermore, traditional leveling methods are inefficient and cannot calculate the optimal adjustment value.

Method used

By grouping the seismic survey lines, extracting the first arrival wave time value from the seabed, calculating the intersection difference, and calculating the adjustment value based on the intersection difference, the survey lines are gradually adjusted using a semi-systematic difference adjustment method until the intersection difference approaches 0, thus achieving rapid leveling.

Benefits of technology

It improves the processing efficiency of seismic line closure error, ensures consistent seismic wave arrival time, enhances the accuracy of underground geological structure imaging and inversion, and supports geological and resource exploration activities.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN117706611B_ABST
    Figure CN117706611B_ABST
Patent Text Reader

Abstract

The application discloses a marine seismic survey line leveling method, a processing system and a storage medium, and comprises the following steps: acquiring a seismic operation survey line; extracting a seabed first arrival wave time value according to the seismic operation survey line; grouping the seismic operation survey line to obtain a main survey line group and a liaison survey line group; obtaining an intersection difference according to the seabed first arrival wave time value, the main survey line group and the liaison survey line group; calculating a difference adjustment value of the seismic operation survey line according to the intersection difference; and completing seismic survey line leveling according to the difference adjustment value. The application can quickly level the seismic survey line, reduce the closure error of the seismic survey line and improve the processing efficiency of the closure error of the seismic survey line by grouping the seismic operation survey line first, then calculating the intersection difference (i.e. the difference value of the seabed first arrival wave time value) of two groups of seismic operation survey lines in different directions, and finally calculating the corresponding difference adjustment value and adjusting all the seismic operation survey lines.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention relates to the field of seismic exploration data processing technology, and in particular to a method, processing system and storage medium for leveling marine seismic survey lines. Background Technology

[0002] Currently, marine seismic survey data obtained through shallow seismic profiling and single-channel seismic measurements often suffer from inconsistent reflection times at the seabed and stratigraphic boundaries at the intersections of different survey lines due to the influence of operational environments such as tides and waves. This inconsistency in the timeframes of the profiles leads to incomplete stratigraphic closure in subsequent seismic interpretations, resulting in stratigraphic discontinuities. Consequently, researchers are unable to accurately determine strata using the survey lines, ultimately affecting the applicability of the data and reducing the accuracy of the survey lines. Therefore, it is necessary to level the seismic survey lines to reduce closure error.

[0003] Current methods for eliminating seismic line closure errors mainly utilize seismic processing systems such as Paradigm and OMEGA. The process involves grouping two intersecting seismic lines together, selecting one group of intersecting lines, observing the discrepancies in the first arrival wave or sea-layer elevation at the intersection points, measuring the adjustment value for the selected line, and then adjusting the corresponding line. This process is repeated for the next group, until all intersecting seismic line groups have been selected and adjusted. This method is inefficient because measuring the discrepancies at the intersection points on the seismic lines is time-consuming; furthermore, since a single seismic line may intersect with multiple lines but only has one adjustment value, this method cannot calculate the optimal adjustment value. Summary of the Invention

[0004] The purpose of this invention is to overcome the problems existing in related technologies. One of the purposes of this invention is to provide a method, processing system and storage medium for leveling marine seismic survey lines, so as to quickly level seismic survey lines, reduce seismic survey line closure error and improve the processing efficiency of seismic survey line closure error.

[0005] The technical solution adopted in this invention is:

[0006] On one hand, embodiments of the present invention provide a method for leveling marine seismic survey lines, comprising the following steps:

[0007] Seismic survey lines are acquired; each seismic survey line is generated based on continuous shot point records.

[0008] Based on the aforementioned seismic survey line, the first arrival wave time value of the seafloor was extracted;

[0009] The seismic survey lines are grouped into main survey line groups and connecting survey line groups, wherein the seismic survey lines in the main survey line groups and the seismic survey lines in the connecting survey line groups have different directions.

[0010] The intersection difference is obtained based on the first arrival wave time value, the main survey line group, and the connecting survey line group; the intersection difference is the difference between the first arrival wave time values ​​at the intersection of the seismic operation survey line in the main survey line group and the seismic operation survey line in the connecting survey line group.

[0011] The adjustment value of the seismic survey line is calculated based on the intersection difference.

[0012] Based on the aforementioned adjustment value, the seismic survey line leveling is completed.

[0013] Further, the step of extracting the first arrival wave time value from the seafloor based on the seismic survey line includes:

[0014] Obtain the first arrival wave from the seabed along the aforementioned seismic survey line;

[0015] Extract the first arrival time value of the seabed and the corresponding geodetic coordinate value from the first arrival time value of the seabed;

[0016] The time values ​​of the first arrival wave on the seabed are subjected to jump point deletion and filtering.

[0017] Furthermore, the grouping of the seismic survey lines to obtain the main survey line group and the connecting survey line group includes:

[0018] All seismic operation survey lines that meet the first condition are assigned to the main survey line group. The first condition is: the direction is the first direction, and there are no intersections between the seismic operation survey lines.

[0019] All seismic operation survey lines that meet the second condition are assigned to the connecting survey line group. The second condition is: the direction is the second direction, and there are no intersections between the seismic operation survey lines; the first direction is different from the second direction.

[0020] Further, obtaining the intersection difference based on the first arrival wave time value, the main survey line group, and the connecting survey line group includes:

[0021] Obtain the first and second work survey lines; the first work survey line comes from the main survey line group; the second work survey line comes from the connecting survey line group;

[0022] Obtain the geodetic coordinates of each shot point on the first and second survey lines, and convert the geodetic coordinates into planar coordinates.

[0023] Obtain the intersection point of the first and second survey lines;

[0024] Based on the plane coordinates, calculate the distance between the shot points on the first working survey line and the shot points on the second working survey line;

[0025] Identify the two shot points with the smallest distance values ​​between the first and second survey lines;

[0026] Based on the two shot points with the smallest distance values ​​and the surrounding shot points, the first and second seabed first arrival wave times are interpolated and calculated; the surrounding shot points are the shot points located around the two shot points with the smallest distance values.

[0027] The first seabed arrival time value is the interpolated first seabed arrival time value of the intersection point on the first working survey line; the second seabed arrival time value is the interpolated first seabed arrival time value of the intersection point on the second working survey line.

[0028] The difference between the first first arrival wave time value and the second first arrival wave time value is calculated as the intersection difference.

[0029] Further, the step of calculating the adjustment value of the seismic survey line based on the intersection difference includes:

[0030] Calculate the differences between all the aforementioned intersection points;

[0031] The arithmetic mean of the differences between all the intersection points is calculated as the adjustment value.

[0032] Furthermore, the step of completing the seismic survey line leveling based on the adjustment value includes:

[0033] The amount of time required to adjust for each of the aforementioned seismic survey lines is set as the current calculation parameter;

[0034] Based on the magnitude of the intersection difference, the current calculation parameter is set to half the magnitude of the intersection difference;

[0035] Adjust the first arrival wave time value based on the current calculation parameters;

[0036] Based on the new first arrival wave time value at the seabed, the intersection difference is recalculated and statistically analyzed;

[0037] Based on the newly calculated intersection difference, reset the current calculation parameters, and then return to adjust the first arrival wave time value on the seabed based on the current calculation parameters until the systematic difference of all intersection differences approaches 0; output the current calculation parameters as the survey line adjustment value;

[0038] Based on the adjustment values ​​of the seismic survey lines, the seismic survey line leveling is completed.

[0039] Further, the step of leveling the seismic survey line based on the survey line adjustment value includes:

[0040] If the seismic survey line belongs to the main survey line group, then the current calculation parameter is set to the survey line adjustment value of the main survey line group;

[0041] If the seismic survey line belongs to the connecting survey line group, then the current calculation parameter plus the adjustment value is set as the survey line adjustment value of the connecting survey line group.

[0042] Furthermore, the method also includes:

[0043] Determine whether the intersection difference exceeds a preset threshold;

[0044] If the intersection difference does not exceed the threshold, then seismic line leveling is not required.

[0045] If the intersection difference exceeds the threshold, then seismic survey line leveling is performed to obtain the survey line adjustment value.

[0046] On the other hand, embodiments of the present invention also provide a marine seismic survey line leveling system, comprising:

[0047] The first module is used to acquire seismic operation survey lines; each seismic operation survey line is generated based on continuous shot point records.

[0048] The second module is used to extract the first arrival wave time value from the seabed based on the seismic survey line.

[0049] The third module is used to group the seismic operation lines into main line groups and connecting line groups, wherein the seismic operation lines in the main line groups and the seismic operation lines in the connecting line groups have different directions.

[0050] The fourth module is used to obtain the intersection difference based on the first arrival wave time value, the main survey line group, and the connecting survey line group; the intersection difference is the difference between the first arrival wave time values ​​at the intersection of the seismic operation survey line in the main survey line group and the seismic operation survey line in the connecting survey line group.

[0051] The fifth module is used to calculate the adjustment value of the seismic survey line based on the intersection difference;

[0052] The sixth module is used to complete the leveling of the seismic survey line based on the aforementioned adjustment value.

[0053] On the other hand, embodiments of the present invention also provide a computer-readable storage medium storing computer-executable instructions for causing a computer to execute the aforementioned marine seismic survey line leveling method.

[0054] The beneficial effects of this invention are: it provides a method, processing system, and storage medium for leveling marine seismic survey lines, which involves acquiring seismic survey lines; each seismic survey line is generated based on continuous shot point records; the first arrival wave time value is extracted from the seismic survey line; the seismic survey lines are grouped into a main survey line group and a connecting survey line group, wherein the seismic survey lines in the main survey line group and the seismic survey lines in the connecting survey line group have different directions; the intersection difference is obtained based on the first arrival wave time value, the main survey line group, and the connecting survey line group; the intersection difference is the difference between the first arrival wave time values ​​at the intersection points of the seismic survey lines in the main survey line group and the seismic survey lines in the connecting survey line group; the adjustment value of the seismic survey line is calculated based on the intersection difference; and the seismic survey line leveling is completed based on the adjustment value. This invention first groups the seismic survey lines, then calculates the intersection difference (i.e., the difference in the first arrival time value of the seafloor wave) between two groups of seismic survey lines in different directions, then calculates the corresponding adjustment value, and finally adjusts all the seismic survey lines. This can quickly level the seismic survey lines, reduce the seismic survey line closure error, and improve the processing efficiency of the seismic survey line closure error. Attached Figure Description

[0055] Figure 1 This is a flowchart of the marine seismic survey line leveling method provided in the embodiments of the present invention;

[0056] Figure 2 This is a flowchart of the marine seismic survey line leveling process provided in an embodiment of the present invention;

[0057] Figure 3 This is a seismic survey line profile provided in an embodiment of the present invention;

[0058] Figure 4 This is a schematic diagram of the first arrival wave time curve provided in an embodiment of the present invention;

[0059] Figure 5 This is a diagram showing the distribution of intersection differences before leveling the seismic profile, provided in an embodiment of the present invention.

[0060] Figure 6 This is a diagram showing the distribution of intersection differences before leveling of the seismic profile, provided in an embodiment of the present invention.

[0061] Figure 7 This is a diagram showing the distribution of intersection differences after leveling the seismic profile, provided in an embodiment of the present invention.

[0062] Figure 8 This is a diagram showing the distribution of intersection differences after leveling the seismic profile, provided in an embodiment of the present invention.

[0063] Figure 9 This is a flight path diagram formed by two seismic survey lines and shot points provided in an embodiment of the present invention;

[0064] Figure 10 This is a schematic diagram of the seismic profile before leveling provided in an embodiment of the present invention;

[0065] Figure 11 This is a schematic diagram of the leveled seismic profile provided in an embodiment of the present invention. Detailed Implementation

[0066] Preferred embodiments of the invention will now be described in more detail with reference to the accompanying drawings. While preferred embodiments of the invention are shown in the drawings, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that the invention will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.

[0067] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The singular forms “a,” “the,” and “the” used in this invention 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 or all possible combinations of one or more of the associated listed items.

[0068] It should be understood that although the terms "first," "second," "third," etc., may be used in this invention to describe various information, this 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, without departing from the scope of this invention, first information may also be referred to as second information, and similarly, second information may also be referred to as first information. Thus, features defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this invention, "a plurality of" means two or more, unless otherwise explicitly specified.

[0069] The embodiments of the present invention will be further described below with reference to the accompanying drawings.

[0070] On one hand, embodiments of the present invention provide a method for leveling marine seismic survey lines. The method of the present invention can ensure that the arrival time of seismic waves in seismic data is consistent across different locations, thereby obtaining more accurate imaging and inversion results of subsurface geological structures, and thus supporting various geological and resource exploration activities.

[0071] This embodiment discloses a method for leveling marine seismic survey lines. Specifically, refer to... Figure 1 The method includes:

[0072] S100. Obtain seismic survey lines; each seismic survey line is generated based on continuous shot point records.

[0073] S200. Based on the seismic survey lines, the time value of the first arrival wave on the seabed is extracted;

[0074] S300. Seismic operation survey lines are grouped into main survey line group and connecting survey line group. The seismic operation survey lines in the main survey line group and the seismic operation survey lines in the connecting survey line group have different directions.

[0075] S400. Based on the first arrival wave time value, the main survey line group, and the connecting survey line group, the intersection difference is obtained. The intersection difference is the difference in the first arrival wave time value at the intersection of the seismic operation survey line in the main survey line group and the seismic operation survey line in the connecting survey line group.

[0076] S500. Based on the intersection difference, the adjustment value of the seismic operation survey line is calculated;

[0077] S600. Based on the adjustment value, complete the leveling of the seismic survey line.

[0078] By executing the process shown in steps S100-S600 above, the purpose of this embodiment of the invention is to: achieve intersection leveling of seismic profiles in a small-scale marine seismic survey network project, reduce seismic profile closure error; ensure that the arrival time of seismic waves in the seismic data is consistent at different locations, so as to obtain more accurate imaging and inversion results of underground geological structures, thereby supporting various geological and resource exploration activities.

[0079] As an optional implementation, in this embodiment of the invention, the seismic survey line is generated by curve fitting based on multiple recorded shot points.

[0080] As an optional implementation, embodiments of the present invention acquire seismic operation survey lines. Marine seismic operation survey lines are seismic exploration activities conducted in the ocean, aiming to obtain subsurface geological information by recording and analyzing the propagation of seismic waves beneath the seabed. These survey lines typically involve generating seismic waves on the ocean surface using equipment such as seismic air guns, and then recording the propagation and reflection of the seismic waves using receiving equipment such as hydrophones. The method of acquiring marine seismic operation survey lines generally includes the following steps:

[0081] Project planning: Determine the research objectives, area, and required seismic survey line density. This requires consideration of factors such as geological characteristics and resource exploration objectives.

[0082] Selecting vessels and instruments: Choose an appropriate research vessel and equip it with the necessary seismic instruments, such as seismic air guns (to generate seismic waves), hydrophones (to receive seismic waves), and seismic recorders (to record data).

[0083] Vessel deployment: Deploy the research vessel to the research area and ensure that the instruments are properly installed and calibrated.

[0084] Use of seismic emission sources: Seismic emission sources, such as seismic air guns, are used at the starting point of the survey line to generate seismic waves. These waves propagate through the water and interact with the seabed strata.

[0085] Hydrophone deployment: Hydrophones are deployed along the survey line, typically suspended in the water or directly installed on the seabed using underwater equipment. Hydrophones are used to receive the propagation and reflection of seismic waves.

[0086] Operation of the seismic recorder: The seismic recorder is used to record seismic data received through the hydrophone, including the arrival time and amplitude of seismic waves.

[0087] Data acquisition and storage: Real-time or post-event acquisition and storage of seismic data. This may involve transmitting data to a seismic recorder on board and performing quality control.

[0088] In marine seismic operations, work is typically conducted along a survey line. The vessel tows a seismic cable equipped with a geophone, and the ship navigates along the pre-designed survey line while simultaneously acquiring data. For shallow exploration, shallow profiling techniques (such as shallow profiling transducers) are often used to transmit and receive signals for data acquisition. In this case, the vessel will also operate along the survey line.

[0089] The purpose of this invention is to level seismic survey lines, a step in geophysical exploration. Its main objective is to correct acquired seismic data so that the arrival times of seismic waves at different locations are consistent. This step is crucial for subsequent seismic data processing and interpretation, and its specific uses include:

[0090] Improving data quality: Leveling seismic survey lines allows for more accurate determination of seismic wave arrival times. This helps improve data quality and reduce errors caused by differences in start times at different locations.

[0091] Enhancing imaging accuracy: In subsurface structure imaging, the arrival time of seismic waves is one of the key factors in generating geological models. Leveling can better align seismic records from different locations in time, contributing to the generation of more accurate geological structure images.

[0092] Reducing uncertainty in stratigraphic interpretation: The arrival time of seismic waves is crucial for interpreting the properties of subsurface strata. By leveling seismic survey lines, the uncertainty in seismic wave arrival times can be reduced, allowing geologists to interpret strata at different depths with greater confidence.

[0093] Reducing data processing complexity: Leveling simplifies subsequent data processing. Aligning the arrival times of seismic waves reduces data processing complexity, making subsequent inversion and model building smoother.

[0094] Improving seismic inversion results: In seismic inversion, leveling helps reduce the time offset of data at different locations, improving the accuracy of the inversion results. This is crucial for identifying underground structures and resource distribution.

[0095] Overall, leveling seismic survey lines is a crucial step in data processing during geophysical exploration. It helps improve data quality and reduce errors, thereby providing a more reliable basis for the accurate interpretation of subsurface structures and resource exploration.

[0096] This embodiment specifically discloses step S200, which involves extracting the first arrival wave time value from the seabed based on the seismic survey line, including:

[0097] S210. Obtain the first arrival wave from the seabed along the seismic survey line;

[0098] S220. Extract the time value of the first arrival wave from the seabed and the corresponding geodetic coordinate value.

[0099] S230. Perform jump point deletion and filtering on the time value of the first arrival wave on the seabed.

[0100] As an optional implementation, step S210 obtains the first arrival wave from the seabed in the seismic survey line: Marine seismic measurement is an acoustic measurement method, and the data is mostly recorded in SEGY format. For SEGY data from shallow seismic profiles, single-channel seismic surveys, etc., the geodetic coordinates of the shot points and the shot time are generally included. Since the acoustic reflection from the seabed is much stronger than that from seawater, the arrival time of the acoustic signal can be extracted based on the amplitude intensity of the signal received by the equipment. The shot number, longitude, latitude, and first arrival wave time value are stored in a text file or database.

[0101] As an optional implementation, the data filtering principle in this embodiment of the invention is as follows: Due to the limited sampling rate of seismic data, such as a sampling rate of 41 μs and a sound speed of 1500 m / s, it is equivalent to one sample point every 41 / 1000 / 1000*1500 = 0.0615 m. If the seabed is relatively flat, there will be obvious stepped patterns, meaning the data accuracy is insufficient and the graphic display is not smooth enough. Therefore, data filtering is required to fit the real seabed topography. For example, IIR Butterworth low-pass filtering is generally used, with the filtering window set to 10 s (i.e., a frequency of 0.1 Hz), which retains the portion of seabed topography changes with a frequency below 0.1 Hz.

[0102] As an optional implementation, the principle of skip point deletion in this embodiment of the invention is as follows: Due to data acquisition reasons, the data acquired by some shot numbers may have strong reflection signals above the seabed, resulting in the extracted seabed first arrival wave time values ​​not being from the seabed. Skip points will appear on the point / curve diagram connecting the seabed first arrival wave time values ​​of all shot numbers, and these need to be deleted. Based on the one-dimensional curve diagram of the extracted seabed first arrival wave time values ​​of the survey line, the seabed first arrival wave time values ​​of the shot points corresponding to these skip points can be deleted from a text document or from self-developed software. After deleting the skip points, the seabed first arrival wave time values ​​of all shot points will be the seismic travel time of the corresponding shot number on the seabed.

[0103] It should be noted that earthquake travel time refers to the time it takes for seismic waves (usually P-waves and S-waves) to travel from the point of earthquake occurrence to different locations on Earth. Earthquake travel time is an important concept in seismology, providing information about the Earth's internal structure.

[0104] As an optional implementation, step S230 will perform jump point deletion and filtering on the seabed first arrival time values: Due to the influence of the working environment and equipment factors, the noise of some shot points in the survey profile may be relatively large, and the extracted seabed first arrival time values ​​may be incorrect, so these jump points need to be removed; Since the sampling interval of the seismic profile record is a finite decimal, which cannot be infinitely small, the extracted seabed first arrival time values ​​may have a sawtooth phenomenon, so it is necessary to select appropriate parameters to filter the data. The filtered seabed first arrival time value curve will be more consistent with the actual seabed topography.

[0105] The first arrival time value is an important parameter for assessing the geological structure of the seabed. Therefore, by analyzing the first arrival time value, this invention can obtain information about the velocity distribution, interfaces, and reflectors of seabed geological bodies. Based on the changes in the first arrival time value, the interface conditions between different strata can be inferred, such as sedimentary layers and lithological variations.

[0106] This embodiment specifically discloses step S300, which involves grouping seismic survey lines to obtain a main survey line group and a connecting survey line group, including:

[0107] S310. All seismic operation lines that meet the first condition are assigned to the main line group. The first condition is: the direction is the first direction and there are no intersections between the seismic operation lines.

[0108] S320. All seismic operation lines that meet the second condition are assigned to the connecting line group. The second condition is: the direction is the second direction, and there are no intersections between the seismic operation lines; the first direction is different from the second direction.

[0109] As an optional implementation, in this embodiment of the invention, step S300 groups the seismic survey lines to obtain main survey line groups and connecting survey line groups. The survey line network is often crisscrossed, with main survey lines in one direction and connecting survey lines in another. There are no intersections between main survey lines or connecting survey lines; intersections generally only exist between main survey lines and connecting lines. To better analyze the differences in intersection points, it is necessary to group all main survey lines into one group and the connecting survey lines into another.

[0110] In this invention, the main survey line groups are typically deployed in key exploration areas. Through a dense network of main survey lines, high-resolution information on subsurface structures can be obtained. This helps to more accurately interpret important features such as the boundaries of geological bodies, reflectors, faults, and lithological variations. Connecting survey line groups can provide data that intersects with or connects to the main survey line groups. By comparing and synthesizing the data with the main survey line group data, the subsurface geological structure can be better interpreted, and noise and uncertainties in the data can be eliminated.

[0111] This embodiment specifically discloses step S400, which obtains the intersection difference based on the first arrival wave time value, the main survey line group, and the connecting survey line group, including:

[0112] S410, Obtain the first and second work survey lines; the first work survey line comes from the main survey line group; the second work survey line comes from the connecting survey line group;

[0113] S420. Obtain the geodetic coordinates of each shot point on the first and second work survey lines, and convert the geodetic coordinates into planar coordinates.

[0114] S430, Obtain the intersection point of the first and second work survey lines;

[0115] S440. Based on the plane coordinates, calculate the distance between the shot points on the first working survey line and the shot points on the second working survey line.

[0116] S450: Obtain the two shot points with the smallest distance between the first and second working survey lines;

[0117] S460. Based on the two shot points with the smallest distance values ​​and the surrounding shot points, interpolate to calculate the first and second seabed first arrival wave times; the surrounding shot points are the shot points located around the two shot points with the smallest distance values.

[0118] The first seabed arrival time value is the interpolated first seabed arrival time value of the intersection point on the first working survey line; the second seabed arrival time value is the interpolated first seabed arrival time value of the intersection point on the second working survey line.

[0119] S470. Calculate the difference between the first and second first arrival wave times at the seabed as the intersection difference.

[0120] As an optional implementation, step S400 of this embodiment of the invention obtains the intersection difference based on the first arrival wave time value, the main survey line group, and the connecting survey line group: if two seismic survey lines intersect, then these two seismic survey lines have an intersection point. Based on the geodetic coordinates of each shot point of the seismic survey line, the geodetic coordinates are converted into plane coordinates to facilitate the calculation of the distance between the two points. Based on the distance values ​​between each shot point of the two seismic survey lines, the pair of shot points with the smallest distance is searched. Then, based on the situation of this pair of shot points and their surrounding shot points, the difference in the first arrival wave time values ​​of the seabed at these two intersection points is interpolated and calculated; this is the intersection difference between the two seismic survey lines.

[0121] This invention, through analysis of intersection differences, can obtain information on velocity variations, interface dip angles, and underground faults between different underground structures. Measuring intersection differences provides a means to verify and correct underground velocity models, thereby improving the accuracy of underground structure models. The analysis of intersection differences can be used to constrain the establishment of underground velocity models. By comparing the first arrival time values ​​at the intersection of the main survey line group and the connecting survey line group, inconsistencies and irregularities in the velocity model can be identified. Further interpretation and correction of the intersection differences can optimize the underground velocity model and improve its interpretability of underground structures.

[0122] This embodiment specifically discloses step S500, which calculates the adjustment value of the seismic survey line based on the intersection difference, including:

[0123] S510, Calculate the difference between all intersection points;

[0124] S520. Calculate the arithmetic mean of the differences between all intersection points as the adjustment value.

[0125] As an optional implementation, in this embodiment of the invention, step S500 calculates the adjustment value of the seismic operation survey line based on the intersection difference: all intersection differences are statistically analyzed, where the intersection difference is the difference between the first arrival wave time values ​​at the intersection of the connecting survey line and the main survey line, and the arithmetic mean of all intersection differences is calculated. This arithmetic mean is the amount of time required to make overall adjustments to the connecting survey line, which is the adjustment value of the seismic operation survey line.

[0126] This invention, through the calculation of adjustment values, can correct systematic errors and biases in the data, thereby improving seismic imaging and positioning results. By correcting intersection differences during the calculation process, the quality and clarity of seismic profiles can be improved, enabling explorers to more accurately locate underground targets, identify geological interfaces, and reveal structural details.

[0127] This invention specifically discloses step S600, which involves leveling the seismic survey line based on the adjustment value, including:

[0128] S610. Set the amount of time that needs to be adjusted for each seismic survey line as the current calculation parameter;

[0129] S620. Based on the magnitude of the intersection difference, set the current calculation parameter to half the magnitude of the intersection difference;

[0130] S630. Adjust the seabed first arrival wave time value according to the current calculation parameters;

[0131] S640. Based on the new first arrival wave time value, recalculate and statistically analyze the intersection difference;

[0132] S650. Based on the newly calculated intersection difference, reset the current calculation parameters, and then return to adjust the seabed first arrival wave time value according to the current calculation parameters until the systematic difference of all intersection differences approaches 0; output the current calculation parameters as the survey line adjustment value.

[0133] S660. Based on the adjustment values ​​of the seismic survey lines, complete the leveling of the seismic survey lines.

[0134] As an optional implementation, in this embodiment of the invention, step S600 completes the seismic line leveling based on the adjustment value, using a semi-systematic difference adjustment method: After adding an adjustment value to the connecting seismic lines, the intersection difference of the specific seismic operation lines is still large, and the seismic profile non-compliance value still has a large value, which is not conducive to subsequent seismic interpretation work such as stratigraphic division. Therefore, it is still necessary to further calculate the adjustment amount of the first arrival wave of the seafloor for each seismic operation line (i.e., the line adjustment value). First, set the time parameter that needs to be adjusted for each seismic operation line (i.e., the current calculation parameter). Based on the size of the intersection difference, set half the adjustment amount for each seismic operation line (i.e., set the current calculation parameter to half the size of the intersection difference), and iterate continuously until the systematic difference of all intersection differences approaches 0.

[0135] According to the measurement line adjustment value, the present invention can determine how to adjust the height of each measuring point, thereby realizing the leveling of the seismic measurement line.

[0136] This invention specifically discloses step S660, which involves leveling the seismic survey line based on the survey line adjustment value, including:

[0137] S670. If the seismic survey line belongs to the main survey line group, then set the current calculation parameters to the survey line adjustment value of the main survey line group.

[0138] S680. If the seismic survey line belongs to the connecting survey line group, then add the adjustment value to the current calculation parameters and set it as the survey line adjustment value of the connecting survey line group.

[0139] As an optional implementation, in this embodiment of the invention, step S660 completes the seismic survey line leveling based on the survey line adjustment value: if the seismic operation survey line is the main survey line, then the adjustment amount calculated this time (i.e., the current calculation parameter) is the final value that the seismic operation survey line needs to be adjusted; if the seismic operation survey line is a connecting survey line, then the adjustment value of the seismic operation survey line is the value calculated this time (i.e., the current calculation parameter) plus the adjustment difference value.

[0140] As an optional implementation, this embodiment of the invention performs seismic profile SEGY adjustment: the seismic SEGY data is continuously recorded shot-by-shot, where each shot-point records acoustic reflection signal data for a fixed time length, and each shot-point records data from multiple sampling points. Adjusting the SEGY data for a seismic profile involves shifting the sample data from each shot-point of the SEGY data upwards by a fixed time value, thereby achieving an overall vertical shift of the entire seismic profile.

[0141] As described in the background section, the disadvantages of the related technology include:

[0142] 1) Low operational efficiency. Because it is necessary to measure the magnitude of discrepancies at the intersection of seismic profiles, a lot of time is required.

[0143] 2) The adjustment value of the survey line profile is estimated manually, and it is impossible to calculate the optimal adjustment value of each survey line profile when there are many intersections of the survey line network.

[0144] To address the aforementioned problems in related technologies, this application proposes a method for leveling seismic profiles using a semi-system difference adjustment, thereby reducing profile closure error and improving work efficiency. The optimal adjustment value for each seismic operation line is calculated based on the principle of minimizing the sum of intersection differences. Furthermore, a method is proposed to calculate the size of seismic profile intersections using the extracted seafloor first arrival wave time values ​​from all shot points of the seismic operation line. This method yields the non-compliance values ​​(i.e., intersection differences) at the intersections of all seismic operation line profiles, which can be used to calculate the profile adjustment value of the seismic operation line.

[0145] Optionally, the marine seismic survey line leveling method specifically disclosed in the embodiments of the present invention further includes:

[0146] S700, Determine whether the intersection difference exceeds the preset threshold;

[0147] S800. If the intersection difference does not exceed the threshold, seismic line leveling is not required.

[0148] S900. If the intersection difference exceeds the threshold, the seismic survey line is leveled to obtain the survey line adjustment value.

[0149] To explain in detail the principle of the technical solution of the present invention, the overall process of the marine seismic survey line leveling method of the present invention will be described below with reference to some specific embodiments. It is easy to understand that the following is an explanation of the technical principle of the present invention and should not be regarded as a limitation of the present invention.

[0150] refer to Figure 2 The detailed steps of the entire marine seismic survey line leveling method are as follows:

[0151] 1) Extraction of first arrival wave from the seabed: Extract the time value and corresponding geodetic coordinate value of the first arrival wave from the SEGY data of the seismic survey line;

[0152] 2) Skip point removal and data filtering: Skip points are removed from the first arrival time values ​​of the seafloor waves extracted from each seismic survey line, and appropriate filtering parameters are selected for filtering, etc.

[0153] 3) Grouping of seismic operation survey lines: All seismic operation survey lines are grouped into main survey lines and connecting survey lines. Generally, there are no intersections between main survey lines and connecting survey lines, but there may be intersections between main survey lines and connecting survey lines.

[0154] 4) Intersection difference calculation: Calculate the non-discrepancy value (referred to as intersection difference) of the seabed first arrival time values ​​for all the seismic survey lines with extracted seabed first arrival time values, and calculate the information of all intersection points;

[0155] 5) Adjustment value calculation: Calculate the average value 'a' of the intersection differences between the connecting survey lines and the main survey line (i.e., the adjustment value), which is the amount of time 'a' required for all connecting survey lines to be moved up or down as a whole, so that the average value of the intersection differences between the connecting survey lines and the main survey line tends to 0; Assuming that each seismic operation survey line has a systematic error (i.e., the current calculation parameter), refer to the semi-systematic error adjustment method in gravity measurement, and follow the principle that the sum of all intersection differences tends to 0, iteratively calculate the adjustment value of the seafloor first arrival wave time value (i.e., the survey line adjustment value) for each seismic operation survey line;

[0156] 6) Calculation of seismic line adjustment value: Add the seismic line adjustment value of each seismic operation line to the adjustment value a to obtain the amount of time required to translate each seismic operation line profile (i.e., the seismic line adjustment value).

[0157] 7) Seismic SEGY profile adjustment: The segy profile data of each seismic operation line is shifted upward or downward to obtain the result data after profile leveling.

[0158] To better compare the advantages and disadvantages of the method of this invention, shallow seismic profile SEGY data from a well site survey project were used for testing. (Reference) Figure 3 Seismic survey line profile, seafloor first arrival time extracted from seismic profile data, reference. Figure 4 The diagram shows the first arrival time curve of the seabed. The green line represents the extracted seabed time, and the red line represents the filtered data. It can be seen that the first arrival time of the seabed before processing fluctuates significantly, up to 0.05ms. The curve after filtering is more consistent with the actual seabed topography.

[0159] After extracting the first arrival times of the seabed for each seismic survey line, the time differences of the first arrival times at the intersections of the seismic survey lines are calculated based on the positional relationships of each shot point. The intersection differences are as follows: Figure 5 Distribution map of intersection differences before seismic profile leveling and Figure 6 The distribution of intersection differences before seismic profile leveling is shown in the diagram.

[0160] from Figure 5 and Figure 6 It can be seen that the intersection points of the seismic survey lines in this embodiment of the invention are relatively large, with the maximum value of the intersection point reaching 2.448 ms (for a two-way trip). If calculated using a sound speed of 1500 m / s, the closure difference between the two survey line profiles at this intersection point is:

[0161]

[0162] The seismic profile closure error is too large and will be unusable in later profile layer interpretations; it needs to be reduced or eliminated.

[0163] The amount of translation required for each seismic profile along the survey line is calculated using this method, as shown in Table 1 below.

[0164] Table 1 Adjustment Amounts for Seismic Profiles on Operational Survey Lines

[0165]

[0166]

[0167] After each seismic survey line is shifted and adjusted, the intersection difference is recalculated to obtain the magnitude of the seismic profile closure error.

[0168] refer to Figure 7 Distribution map of intersection differences after seismic profile leveling and Figure 8 The distribution map of intersection differences after seismic profile leveling, from Figure 7 , Figure 8 The intersection difference distribution data shows that after leveling each seismic survey line, the maximum value, standard error, systematic error, and average error of the difference in seafloor first arrival time values ​​at the intersection points of the seismic survey lines all significantly decreased. The maximum difference in seafloor first arrival time at the intersection point of the largest seismic profile after leveling was only 0.124 ms, which translates to a depth difference of only 0.093 m, greatly improving accuracy and no longer affecting subsequent seismic interpretation. The specific formula for converting intersection difference to height difference is as follows:

[0169]

[0170] To better compare the seismic profile data after profile leveling, survey lines 21ZH1L05 and 21ZH1L10 were selected for comparison.

[0171] refer to Figure 9 The embodiment of this invention provides a flight path diagram formed by two seismic survey lines and shot points, wherein survey lines 21ZH1L05 and 21ZH1L10 have different directions; survey line 21ZH1L05 is oriented vertically downwards, while survey line 21ZH1L10 is oriented horizontally to the left; and then refer to... Figure 10 The seismic profile diagram before leveling provided in the embodiments of the present invention and Figure 11 As can be seen from the schematic diagram of the leveled seismic profile provided in this embodiment of the invention, before the profile was leveled, the closure error of the profile at the intersection point was large and there was obvious misalignment. After the seismic profile was leveled by this method, the closure error of the seismic profile was basically eliminated.

[0172] By comparing the difference in the first arrival wave at the intersection of the seismic survey lines before and after leveling using the method of this embodiment of the invention, and the profile diagram of the seismic survey lines, it can be seen that using the method of this embodiment of the invention to level the seismic survey line profile can achieve efficient operation, greatly reduce the closure error of each seismic survey line profile, and meet the needs of subsequent seismic data interpretation.

[0173] Compared with the prior art, the technical solution of the present invention has the following advantages:

[0174] 1. Using the difference in the first arrival time of the seabed at the intersection of each seismic operation line (referred to as the intersection difference), the overall adjustment value of the connecting lines is first calculated. Then, based on the overall adjusted intersection difference, and taking the principle of minimizing the arithmetic mean of all intersection differences, the amount that each seismic operation line still needs to be adjusted is iteratively calculated. Finally, the time amount that each seismic operation line ultimately needs to be adjusted is calculated, thereby achieving the effect of reducing the seismic profile closure error.

[0175] 2. Utilizing the strong acoustic reflection signal from the seabed, calculate the first arrival time value of the seabed wave for all seismic operation lines, that is, obtain the geodetic coordinates and first arrival time value of the seabed wave for each shot point of each seismic operation line.

[0176] 3. Divide the seismic operation survey lines into main survey lines and connecting survey lines. According to the distance relationship between each point, compare and calculate point by point to find the intersection of the seismic operation survey lines. Then, interpolate and calculate the difference in the first arrival time value of the seabed of the seismic operation survey lines at the theoretical intersection point (intersection difference).

[0177] 4. The extracted first arrival wave data of the seabed along the seismic survey line is processed by removing jump points and filtering the data to make the first arrival wave curve more consistent with the actual seabed topography.

[0178] On the other hand, embodiments of the present invention provide a marine seismic survey line leveling system, comprising:

[0179] The first module is used to acquire seismic operation lines; each seismic operation line is generated based on continuous shot point records.

[0180] The second module is used to extract the first arrival wave time value from the seismic survey line.

[0181] The third module is used to group the seismic operation lines into main line groups and connecting line groups. The seismic operation lines in the main line groups and the seismic operation lines in the connecting line groups have different directions.

[0182] The fourth module is used to obtain the intersection difference based on the first arrival wave time value, the main survey line group, and the connecting survey line group; the intersection difference is the difference in the first arrival wave time value at the intersection of the seismic operation survey line in the main survey line group and the seismic operation survey line in the connecting survey line group.

[0183] The fifth module is used to calculate the adjustment value of the seismic survey line based on the intersection difference;

[0184] The sixth module is used to level the seismic survey lines based on the adjustment values.

[0185] It is understood that the content of the above method embodiments is applicable to the system embodiments. The specific functions implemented in the system embodiments are the same as those in the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above methods.

[0186] The marine seismic survey line leveling processing system of this invention includes a memory and a processor.

[0187] The processor can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor can be a microprocessor or any conventional processor.

[0188] Memory can include various types of storage units, such as system memory, read-only memory (ROM), and permanent storage devices. ROM can store static data or instructions required by the processor or other modules of the computer. Permanent storage devices can be read-write storage devices. Permanent storage devices can be non-volatile storage devices that retain stored instructions and data even when the computer is powered off. In some embodiments, permanent storage devices use high-capacity storage devices (e.g., magnetic or optical disks, flash memory) as permanent storage devices. In other embodiments, permanent storage devices can be removable storage devices (e.g., floppy disks, optical drives). System memory can be a read-write storage device or a volatile read-write storage device, such as dynamic random access memory. System memory can store some or all of the instructions and data required by the processor during operation. Furthermore, memory can include any combination of computer-readable storage media, including various types of semiconductor memory chips (DRAM, SRAM, SDRAM, flash memory, programmable read-only memory), and disks and / or optical disks can also be used. In some implementations, the memory may include removable storage devices that are readable and / or writable, such as laser discs (CDs), read-only digital versatile optical discs (e.g., DVD-ROMs, dual-layer DVD-ROMs), read-only Blu-ray discs, ultra-high density optical discs, flash memory cards (e.g., SD cards, mini SD cards, Micro-SD cards, etc.), magnetic floppy disks, etc. Computer-readable storage media do not contain carrier waves or transient electronic signals transmitted wirelessly or via wired connections.

[0189] The memory stores executable code, which, when processed by the processor, can cause the processor to execute some or all of the methods described above.

[0190] On the other hand, embodiments of the present invention provide a computer-readable storage medium storing computer-executable instructions for causing a computer to perform the marine seismic survey line leveling method as described above.

[0191] In summary, the technical solution of this invention provides a method, processing system, and storage medium for leveling marine seismic survey lines. The method calculates the first arrival wave time values ​​at all seafloor points of the seismic survey lines, along with the magnitude and statistical information of the differences at each intersection point. Based on the distribution and statistical information of the intersection point differences, the overall adjustment value of the connecting survey lines is calculated. Then, based on the overall adjusted intersection point difference values, the adjustment value of the first arrival wave time value of each seismic survey line is iteratively calculated using a semi-systematic difference adjustment method, i.e., the amount of adjustment still needed for each seismic survey line. This allows for the calculation of the final adjustment time required for each seismic survey line, enabling the adjustment of all seismic survey lines. This method has good processing efficiency, first adjusting the connecting lines as a whole, then adjusting each line individually. The adjusted survey lines have small closure errors, facilitating researchers to accurately determine the bottom layer using the survey lines.

[0192] It will be understood by those skilled in the art that all or some of the steps and systems in the methods disclosed above can be implemented as software, firmware, hardware, and suitable combinations thereof. Some or all of the physical components can be implemented as software executed by a processor, such as a central processing unit, digital signal processor, or microprocessor, or as hardware, or as an integrated circuit, such as an application-specific integrated circuit. Such software can be distributed on a computer-readable medium, which can include computer storage media (or non-transitory media) and communication media (or transient media). As is known to those skilled in the art, the term computer storage media includes volatile and non-volatile, removable and non-removable media implemented in any method or technology for storing information (such as computer-readable instructions, data structures, program modules, or other data). Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technologies, CD-ROM, digital versatile disc (DVD) or other optical disc storage, magnetic cartridges, magnetic tape, disk storage or other magnetic storage devices, or any other medium that can be used to store desired information and is accessible to a computer. Furthermore, as is known to those skilled in the art, communication media typically contain computer-readable instructions, data structures, program modules, or other data in modulated data signals such as carrier waves or other transmission mechanisms, and may include any information delivery medium.

[0193] Unless otherwise specifically stated, the relative arrangement, numerical expressions, and values ​​of the components and steps described in these embodiments do not limit the scope of this application. It should also be understood that, for ease of description, the dimensions of the various parts shown in the drawings are not drawn to actual scale. Techniques, methods, and devices known to those skilled in the art may not be discussed in detail, but where appropriate, such techniques, methods, and devices should be considered part of the specification. In all examples shown and discussed herein, any specific values ​​should be interpreted as merely exemplary and not as limitations. Therefore, other examples of exemplary embodiments may have different values. It should be noted that similar reference numerals and letters in the following drawings denote similar items; therefore, once an item is defined in one drawing, it need not be further discussed in subsequent drawings.

[0194] In the description of this application, it should be understood that the orientation or positional relationship indicated by directional terms such as "front, back, up, down, left, right", "horizontal, vertical, horizontal" and "top, bottom" is usually based on the orientation or positional relationship shown in the accompanying drawings, and is only for the convenience of describing this application and simplifying the description. Unless otherwise stated, these directional terms do not indicate or imply that the device or element referred to must have a specific orientation or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation on the scope of protection of this application; the directional terms "inner" and "outer" refer to the inner and outer contours relative to the outline of each component itself.

[0195] For ease of description, spatial relative terms such as "above," "on top of," "on the upper surface of," "above," etc., are used herein to describe the spatial positional relationship of a device or feature as shown in the figures to other devices or features. It should be understood that spatial relative terms are intended to encompass different orientations in use or operation beyond the orientation of the device as described in the figures. For example, if the device in the figures were inverted, a device described as "above" or "on top of" other devices or structures would subsequently be positioned as "below" or "under" other devices or structures. Thus, the exemplary term "above" can include both "above" and "below." The device may also be positioned in other different ways (rotated 90 degrees or in other orientations), and the spatial relative descriptions used herein will be interpreted accordingly.

[0196] The above is a detailed description of the preferred embodiments of the present invention. However, the present invention is not limited to the above embodiments. Those skilled in the art can make various equivalent modifications or substitutions without departing from the spirit of the present invention. All such equivalent modifications or substitutions are included within the scope defined by the claims of the present invention.

Claims

1. A method for leveling marine seismic survey lines, characterized in that, Includes the following steps: Seismic survey lines are acquired; each seismic survey line is generated based on continuous shot point records. Based on the aforementioned seismic survey line, the first arrival wave time value of the seafloor was extracted; The seismic survey lines are grouped into main survey line groups and connecting survey line groups, wherein the seismic survey lines in the main survey line groups and the seismic survey lines in the connecting survey line groups have different directions. The intersection difference is obtained based on the first arrival wave time value on the seabed, the main survey line group, and the connecting survey line group; The intersection difference is the difference in the first arrival time value of the seafloor wave at the intersection of the seismic operation survey line in the main survey line group and the seismic operation survey line in the connecting survey line group. The adjustment value of the seismic survey line is calculated based on the intersection difference. Based on the aforementioned adjustment value, the seismic survey line leveling is completed.

2. The method for leveling marine seismic survey lines according to claim 1, characterized in that, The extraction of the first arrival time value of the seafloor wave based on the seismic survey line includes: Obtain the first arrival wave from the seabed along the aforementioned seismic survey line; Extract the first arrival time value of the seabed and the corresponding geodetic coordinate value from the first arrival time value of the seabed; The time values ​​of the first arrival wave on the seabed are subjected to jump point deletion and filtering.

3. The method for leveling marine seismic survey lines according to claim 1, characterized in that, The process of grouping the seismic survey lines to obtain main survey line groups and connecting survey line groups includes: All seismic operation survey lines that meet the first condition are assigned to the main survey line group. The first condition is: the direction is the first direction, and there are no intersections between the seismic operation survey lines. All seismic operation survey lines that meet the second condition are assigned to the connecting survey line group. The second condition is: the direction is the second direction, and there are no intersections between the seismic operation survey lines; the first direction is different from the second direction.

4. The method for leveling marine seismic survey lines according to claim 3, characterized in that, The step of obtaining the intersection difference based on the first arrival wave time value, the main survey line group, and the connecting survey line group includes: Obtain the first and second work survey lines; the first work survey line comes from the main survey line group; the second work survey line comes from the connecting survey line group; Obtain the geodetic coordinates of each shot point on the first and second survey lines, and convert the geodetic coordinates into planar coordinates. Obtain the intersection point of the first and second survey lines; Based on the plane coordinates, calculate the distance between the shot points on the first working survey line and the shot points on the second working survey line; Identify the two shot points with the smallest distance values ​​between the first and second survey lines; Based on the two shot points with the smallest distance values ​​and the surrounding shot points, the first and second seabed first arrival wave times are interpolated and calculated; the surrounding shot points are the shot points located around the two shot points with the smallest distance values. The first seabed arrival time value is the interpolated first seabed arrival time value of the intersection point on the first working survey line; the second seabed arrival time value is the interpolated first seabed arrival time value of the intersection point on the second working survey line. The difference between the first first arrival wave time value and the second first arrival wave time value is calculated as the intersection difference.

5. The method for leveling marine seismic survey lines according to claim 1, characterized in that, The step of calculating the adjustment value of the seismic survey line based on the intersection difference includes: Calculate the differences between all the aforementioned intersection points; The arithmetic mean of the differences between all the intersection points is calculated as the adjustment value.

6. The method for leveling marine seismic survey lines according to claim 5, characterized in that, The step of leveling the seismic survey line based on the adjustment value includes: The amount of time required to adjust for each of the aforementioned seismic survey lines is set as the current calculation parameter; Based on the magnitude of the intersection difference, the current calculation parameter is set to half the magnitude of the intersection difference; Adjust the first arrival wave time value based on the current calculation parameters; Based on the new first arrival wave time value at the seabed, the intersection difference is recalculated and statistically analyzed; Based on the newly calculated intersection difference, reset the current calculation parameters, and then return to adjust the first arrival wave time value on the seabed based on the current calculation parameters until the systematic difference of all intersection differences approaches 0; output the current calculation parameters as the survey line adjustment value; Based on the adjustment values ​​of the seismic survey lines, the seismic survey line leveling is completed.

7. The method for leveling marine seismic survey lines according to claim 6, characterized in that, The adjustment value is based on the measured line. Complete the leveling of seismic survey lines, including: If the seismic survey line belongs to the main survey line group, then the current calculation parameter is set to the survey line adjustment value of the main survey line group; If the seismic survey line belongs to the connecting survey line group, then the current calculation parameter plus the adjustment value is set as the survey line adjustment value of the connecting survey line group.

8. The method for leveling marine seismic survey lines according to claim 6, characterized in that, The method further includes: Determine whether the intersection difference exceeds a preset threshold; If the intersection difference does not exceed the threshold, then seismic line leveling is not required. If the intersection difference exceeds the threshold, then seismic survey line leveling is performed to obtain the survey line adjustment value.

9. A marine seismic survey line leveling and processing system, characterized in that, include: The first module is used to acquire seismic operation survey lines; each seismic operation survey line is generated based on continuous shot point records. The second module is used to extract the first arrival wave time value from the seabed based on the seismic survey line. The third module is used to group the seismic operation lines into main line groups and connecting line groups, wherein the seismic operation lines in the main line groups and the seismic operation lines in the connecting line groups have different directions. The fourth module is used to obtain the intersection difference based on the first arrival wave time value on the seabed, the main survey line group, and the connecting survey line group; The intersection difference is the difference in the first arrival time value of the seafloor wave at the intersection of the seismic operation survey line in the main survey line group and the seismic operation survey line in the connecting survey line group. The fifth module is used to calculate the adjustment value of the seismic survey line based on the intersection difference; The sixth module is used to complete the leveling of the seismic survey line based on the aforementioned adjustment value.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer-executable instructions for causing a computer to perform the method as described in any one of claims 1 to 8.