CN114993269B

CN114993269B - Method for extracting track water depth by using SEGY data

Info

Publication number: CN114993269B
Application number: CN202210386364.2A
Authority: CN
Inventors: 曾凡祥; 张旭东; 杨力; 王启
Original assignee: Guangzhou Marine Geological Survey; Southern Marine Science and Engineering Guangdong Laboratory Guangzhou
Current assignee: Guangzhou Marine Geological Survey; Southern Marine Science and Engineering Guangdong Laboratory Guangzhou
Priority date: 2022-04-13
Filing date: 2022-04-13
Publication date: 2023-08-18
Anticipated expiration: 2042-04-13
Also published as: CN114993269A

Abstract

The application relates to a method for extracting track water depth by utilizing SEGY data. The method comprises the following steps: performing data processing on SEGY data of each measuring point on the track to obtain the coordinates of the measuring point and the water depth of the measuring point of each measuring point, and forming the water depth of the track; the process for acquiring the coordinates of a measuring point and the water depth of the measuring point comprises the following steps: performing de-compiling on SEGY data of the measuring points to obtain GPS antenna head coordinates and sample point amplitude sets; acquiring the draught depth of a transducer, the sound velocity value of sea water and the starting moment of submarine searching; carrying out coordinate correction based on the GPS antenna head coordinates to obtain measuring point coordinates; screening out the maximum sample point amplitude from the sample point amplitude set; and calculating according to the seabed recording time, the seabed searching starting time, the transducer draft and the seawater sound velocity value corresponding to the maximum sampling point amplitude to obtain the water depth of the measuring point. The scheme provided by the application can be obtained based on the existing SEGY data processing without reconnaissance, so that the time is saved and the reconnaissance cost is greatly saved.

Description

Method for extracting track water depth by using SEGY data

Technical Field

The application relates to the technical field of marine survey, in particular to a method for extracting track water depth by using SEGY data.

Background

When marine survey is performed at present, the acquisition of the water depth value is mainly obtained by acquiring data of a depth finder, and the water depth data acquired by the depth finder is more accurate under the general condition, however, when the fault, aging or parameter setting of the observation equipment of the depth finder is wrong, the accuracy of the observed water depth value is poor and the requirements of corresponding projects cannot be met, and due to the high cost of marine measurement and retesting, the new observation equipment of the depth finder is utilized for retesting the water depth value, so that the time is consumed, and the additional expenditure of large survey cost is brought.

The related art provides a time-depth conversion method, which needs to manually pick up or manually interact with automatic pick-up parameters on a seismic system or a seismic image by each line, and specifically: and manually picking up the first arrival time by using a manual picking method on the seismic system or the seismic image by using each measuring line, and then performing time depth conversion to obtain the depth of the submarine horizon.

The manual picking method can accurately acquire the submarine recording delay value, but needs more time to process, and further processes are needed to be carried out on each measuring line jump point after the submarine delay is acquired, so that the submarine reflection signal recording time can be converted into the water depth value, and the efficiency is low.

Disclosure of Invention

In order to overcome the problems in the related art, the application provides a method for extracting the track water depth by using the SEGY data, which can be obtained based on the existing SEGY data processing, and greatly saves the survey cost while saving the time.

The first aspect of the application provides a method for extracting track water depth by using SEGY data, which comprises the following steps:

performing data processing on SEGY data of each measuring point on the track to obtain the coordinates of the measuring point and the water depth of the measuring point of each measuring point, thereby forming the water depth of the track; the track water depth is a mapping data set of the coordinates of the measuring points and the water depth of the measuring points of all the measuring points on the track;

the process for obtaining the coordinates of each measuring point and the water depth of the measuring point in the measuring point coordinates and the water depth of the measuring point by performing data processing on the SEGY data of each measuring point on the track comprises the following steps:

performing de-compiling on the SEGY data of the measuring points to obtain GPS antenna head coordinates and sample point amplitude sets;

acquiring the draught depth of a transducer, the sound velocity value of sea water and the starting moment of submarine searching;

carrying out coordinate correction based on the GPS antenna head coordinates to obtain the measuring point coordinates;

screening out the maximum sample point amplitude from the sample point amplitude set;

and calculating according to the seabed recording time, the seabed searching starting time, the transducer draft and the seawater sound velocity value corresponding to the maximum sample point amplitude to obtain the measuring point water depth.

In one embodiment, the screening the maximum sample amplitude from the set of sample amplitudes includes:

and searching the sampling point amplitude set based on the principle that the sum of the sampling point amplitudes of the window is maximum according to the preset time window size, and obtaining the maximum sampling point amplitude and the corresponding submarine recording moment.

In one embodiment, according to the preset time window size, searching the set of sample amplitudes based on the principle that the sum of the sample amplitudes of the window is maximum, to obtain the maximum sample amplitude and the corresponding submarine recording time, includes:

searching the sample point amplitude set by taking the sample point corresponding to the seabed searching starting moment as a starting point, taking the sample point at the latest recording moment in the sample point amplitude set as an end point, gradually moving backward by a preset sliding step length, and recording the sample point amplitude sum in each time window to obtain a window sample point amplitude sum data set;

and selecting and obtaining the maximum window sample point amplitude sum from the window sample point amplitude sum data set, taking the maximum amplitude in the maximum window sample point amplitude sum as the maximum sample point amplitude, and taking the moment corresponding to the maximum sample point amplitude as the submarine recording moment.

In one embodiment, the calculating the measuring point water depth according to the seabed recording time corresponding to the maximum sample amplitude, the seabed searching starting time, the transducer draft and the seawater sound velocity value includes:

calculating according to the following calculation formula to obtain the water depth of the measuring point;

wherein D represents the water depth of the measuring point; v represents the sound velocity value of the sea water; d represents the transducer draft; and T represents the sound wave propagation time, wherein the sound wave propagation time is the interval duration between the submarine recording moment and the submarine searching starting moment.

In one embodiment, after obtaining the coordinates and the water depth of each measuring point, the method further includes:

and filtering the water depth value of the measuring point of each measuring point to obtain the track water depth.

In one embodiment, the filtering the water depth value of the measuring point of each measuring point includes:

acquiring a water depth value filtering threshold value;

and removing the water depths of the measuring points, which are larger than the water depth value filtering threshold value, in the water depths of the measuring points of each measuring point.

acquiring a preset seabed amplitude value;

and removing the water depth of the measuring point of which the water depth is smaller than 1/4 value of the preset seabed amplitude value.

In one embodiment, after the maximum sample amplitude is selected from the set of sample amplitudes, the method includes:

obtaining a maximum amplitude sampling point ordinal number of a sampling point corresponding to the maximum sampling point amplitude;

determining the submarine recording moment according to the maximum amplitude sampling point ordinal number, the recording delay value and the data sampling interval; the record delay value and the data sampling interval are obtained by the SEGY data de-encoding.

In one embodiment, the determining the submarine recording moment according to the maximum amplitude sample number, the recording delay value and the data sampling interval comprises:

calculating according to the following calculation formula to obtain the recording interval duration;

Δt＝N*a+b；

wherein Δt represents the recording interval duration; n represents the maximum amplitude sample point ordinal number; a represents a data sampling interval; b represents a recording delay value;

and combining the seabed searching starting time and the recording interval time to obtain the seabed recording time.

In one embodiment, before the coordinate correction based on the GPS antenna head coordinates, the method includes:

respectively carrying out position detection on the GPS antenna head and the transducer under the same state to obtain the actual measurement coordinates of the GPS antenna head and the actual measurement coordinates of the transducer under the same coordinate system;

calculating to obtain a coordinate offset according to the GPS antenna head measured coordinate and the transducer measured coordinate; the coordinate offset includes a horizontal axis offset and a vertical axis offset;

and correcting the GPS antenna head coordinate according to the coordinate offset to obtain the measuring point coordinate.

The technical scheme provided by the application can comprise the following beneficial effects:

according to the method for extracting the track water depth by using the SEGY data, the GPS antenna head coordinates and the sample point amplitude set are obtained based on SEGY data of the measuring points, and the deviation between the positions of the GPS antenna head and the transducer, which is caused by the large volume of a survey ship, is corrected through coordinate correction, so that the accurate measuring point coordinates are obtained; the method comprises the steps that a maximum sampling point amplitude corresponding sampling point is searched in a sampling point amplitude set, namely, a sampling point of submarine data is surveyed, the data acquisition time corresponding to the sampling point, namely, the submarine recording time can be determined, so that the time for transmitting sound waves from a transducer back to the transducer through the submarine can be known according to the submarine recording time and the submarine searching starting time, the depth between the submarine position of the measuring point and the transducer can be calculated by combining a seawater sound velocity value, and the measuring point water depth of the measuring point can be obtained by adding the draught depth of the transducer on the basis; because the acquisition process of the water depth of the measuring point is based on the existing SEGY data processing, reconnaissance is not needed, and the time is saved and the reconnaissance cost is greatly saved; and the maximum and continuous amplitude of the submarine signal reflection in the SEGY data is considered, so that the recording time of the maximum submarine amplitude is obtained by using the SEGY data, and the water depths of all measuring points on the track are calculated, so that the data are more accurate and the processing efficiency is higher compared with the manual pick-up process.

It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the application as claimed.

Drawings

The foregoing and other objects, features and advantages of the application will be apparent from the following more particular descriptions of exemplary embodiments of the application as illustrated in the accompanying drawings wherein like reference numbers generally represent like parts throughout the exemplary embodiments of the application.

FIG. 1 is a schematic flow chart of a process for acquiring coordinates and water depth of a measuring point according to an embodiment of the present application;

FIG. 2 is a flow chart of a method for screening for maximum sample amplitude according to an embodiment of the present application;

FIG. 3 is a schematic flow chart of a track water depth extraction method for removing water depth jump points according to an embodiment of the application;

fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present application.

Detailed Description

Preferred embodiments of the present application will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present application are shown in the drawings, it should be understood that the present application may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the application to those skilled in the art.

The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used in this specification and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any or all possible combinations of one or more of the associated listed items.

It should be understood that although the terms "first," "second," "third," etc. may be used herein to describe various information, these information should not be limited by these terms. These terms are only used to distinguish one type of information from another. For example, first information may also be referred to as second information, and similarly, second information may also be referred to as first information, without departing from the scope of the application. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the present application, the meaning of "a plurality" is two or more, unless explicitly defined otherwise.

Example 1

When the fault, aging or parameter setting of the sounding device is wrong, the observed water depth value precision is poor, the requirements of corresponding projects cannot be met, and due to the high cost of the marine survey re-survey, the new sounding device is used for re-survey of the water depth value, so that time is consumed, and additional expenditure of large survey cost is brought.

The time-depth conversion method provided by the related technology needs to manually pick up or manually interact with automatic pick-up parameters on the seismic system or the seismic image line by line, and has low efficiency.

Aiming at the problems, the embodiment of the application provides a method for extracting the track water depth by using the SEGY data, which can be obtained based on the existing SEGY data processing without reconnaissance, and greatly saves the reconnaissance cost while saving time.

The following describes the technical scheme of the embodiment of the present application in detail with reference to the accompanying drawings.

FIG. 1 is a flow chart of a process for acquiring coordinates and water depth of a measurement point according to an embodiment of the present application.

A method for extracting track water depth by using SEGY data, comprising: and performing data processing on the SEGY data of each measuring point on the track to obtain the measuring point coordinates and the measuring point water depth of each measuring point, thereby forming the track water depth.

The track water depth is a mapping data set of the coordinates of the measuring points and the water depth of the measuring points of all the measuring points on the track; the water depth of the measuring point on each measuring point coordinate position on the track can be obtained according to the water depth of the track, and then a water depth curve on the track is formed.

Referring to fig. 1, in an embodiment of the present application, a process for acquiring coordinates of a measurement point and a water depth of the measurement point includes:

101. performing de-compiling on SEGY data of the measuring points to obtain GPS antenna head coordinates and sample point amplitude sets;

the SEGY data is seismic data organized by taking a seismic channel as a unit, the storage mode of each data body is to store the amplitude value of one sampling point corresponding to each time point, the SEGY file is a binary data file, and the data is acquired from the SEGY, so that the corresponding seismic channel record information and the corresponding sampling point data can be read by means of de-compiling. Therefore, after the SEGY data is decompressed, a group of sample amplitude data of the current measuring point, namely a sample amplitude set, can be obtained.

102. Acquiring the draught depth of a transducer, the sound velocity value of sea water and the starting moment of submarine searching;

in the embodiment of the application, the draught depth of the transducer, the sound velocity value of the sea water and the seabed searching starting time can be a parameter set pre-stored in a memory, and can also be obtained by inputting in real time by an operator.

It should be noted that, the execution timing of step 102 is not strictly limited in the embodiment of the present application, and in the practical application process, step 102 may be executed before step 105, that is, step 102 may also be executed before step 101 or between step 103 and step 105.

103. Carrying out coordinate correction based on the GPS antenna head coordinates to obtain measuring point coordinates;

in the process of acquiring the SEGY data, a survey ship has a certain volume, a certain deviation exists between the position of the GPS antenna head and the position of the transducer, and the coordinates recorded in the SEGY data are generally the position of the GPS antenna head and are not the real position of the transducer, so that in order to ensure the data accuracy of the track water depth, the track is often required to be fitted according to the coordinate data of the measuring line, the course of the measuring line is calculated, and the position translation is carried out according to the relative position and the course angle of the transducer and the GPS antenna head, so that the coordinates of the real shallow-profile transducer are obtained.

Therefore, the GPS antenna head coordinate obtained by the SEGY data de-encoding in step 101 cannot be directly used as the measurement point coordinate of the current measurement point, and the GPS antenna head coordinate needs to be corrected by a coordinate translation method or the like, so that an accurate measurement point coordinate is obtained.

In the embodiment of the application, because the survey ship is fixed in the process of surveying a track, the deviation value between the position of the GPS antenna head and the position of the transducer is relatively fixed; on the basis, after the actual positions of the GPS antenna head and the transducer in the same state are detected, the offset is calculated, and the GPS antenna head coordinates under each measuring point can be corrected by using the offset.

Specific:

correcting the coordinates of the GPS antenna head according to the coordinate offset to obtain coordinates of the measuring point; the method comprises the steps of adding a transverse axis offset on the basis of a GPS antenna head transverse coordinate to obtain a transverse coordinate of a measuring point coordinate; and adding the vertical axis offset on the basis of the vertical coordinate of the GPS antenna head to obtain the vertical coordinate of the measuring point coordinate.

104. Screening out the maximum sample point amplitude from the sample point amplitude set;

in the surveying process, the current measuring point depth is the depth of the corresponding submarine position of the measuring point coordinate, namely the point with the longest acoustic wave propagation path, and the maximum sample point amplitude is obtained from the sample point amplitude data.

In the traditional time-depth conversion method, the maximum sample point amplitude is extracted root by mainly adopting a manual picking method and a method for automatically picking up first arrival wave time delay through man-machine interaction, so that the efficiency is low and the error of extracted data is high.

In the embodiment of the application, a time window with a preset size is utilized, and a sampling point amplitude set is searched based on the principle that the sum of sampling point amplitudes of the window is maximum, so that the maximum sampling point amplitude is obtained, and the submarine recording moment corresponding to the maximum sampling point amplitude is extracted to be used for calculating the water depth of the measuring point of the current measuring point.

105. And calculating according to the seabed recording time, the seabed searching starting time, the transducer draft and the seawater sound velocity value corresponding to the maximum sampling point amplitude to obtain the water depth of the measuring point.

The process of surveying the water depth by utilizing sound waves is that the sound waves are sent out by the transducer position and then reflected back to the transducer position through the sea bottom, the moment of sending out the sound waves is taken as the sea bottom searching starting moment, the sound waves reflected back by the sea bottom are received by the transducer as the sea bottom recording moment corresponding to the maximum sample point amplitude, the interval duration before the sea bottom searching starting moment and the sea bottom recording moment is the duration required by the sound waves to make a round trip between the transducer position and the sea bottom position, and the distance between the transducer position and the sea bottom position can be obtained by combining the propagation speed of the sound waves in sea water, namely the sea water sound velocity value, and the distance from the current measuring point coordinate water surface to the current measuring point coordinate sea bottom, namely the measuring point water depth, can be obtained by adding the draught of the transducer.

Specifically, the water depth of the measuring point can be calculated according to the following calculation formula:

It should be noted that, the propagation speed of the sound wave in the sea will be affected by the sea condition, so the sound speed value of the sea adopted in the embodiment of the application is the actually measured sound speed value.

Example two

Step 104 in the first embodiment of the present application is designed, in which the sample amplitude set is searched by using time window sliding, the maximum sum of sample amplitudes in the time window is used as a screening principle, a time window area with the maximum sum of sample amplitudes is found, the maximum value in the time window area is used as the maximum sample amplitude in the seismic data of the current measuring point, and the screening speed is improved by using a time window sliding screening mode, so as to achieve the purpose of efficient data extraction.

FIG. 2 is a flow chart of a method for screening for maximum sample amplitude according to an embodiment of the present application.

Referring to fig. 2, the screening method of the maximum sample amplitude includes:

201. searching the sample point amplitude set by taking the sample point corresponding to the seabed searching starting moment as a starting point and taking the sample point at the latest recording moment in the sample point amplitude set as an end point, gradually moving backward by a preset sliding step length, and recording the sample point amplitude sum in each time window to obtain a window sample point amplitude sum data set;

in the embodiment of the present application, the time window size is set according to the actual requirements of searching fineness and searching speed, and in the embodiment of the present application, the time window size is set as the pulse length of the SEGY profile data, for example: when one pulse length is 1ms, the time window size is set to 1ms.

The searching process in the embodiment of the application is as follows:

taking a sample point corresponding to the seabed searching starting moment as a starting point of time window movement, and calculating the sum of all sample point amplitudes in the current time window as a window sample point amplitude sum S1 when the time window is positioned at the starting point; then, the time window is moved backwards by a preset sliding step length, and the sum S2 of all sample amplitudes in the time window at the moment is calculated; and similarly, after each time window is moved backwards by a preset sliding step length, calculating to obtain the sum of all sample point amplitudes in the current time window until the time window moves to a sample point at the latest recording moment in the sample point amplitude set, and obtaining the sum Sn of all sample point amplitudes in the current time window; s1 to Sn form the window sample amplitude sum dataset described above.

It should be noted that, the preset sliding step length is a numerical value set by the operator according to the actual searching requirement, and may be adjusted according to the actual situation, which is not limited herein.

In the practical application process, the moving direction of the time window is not unique, or the sample point corresponding to the starting moment of the submarine search is taken as the end point, and the sample point at the latest recording moment in the sample point amplitude set is taken as the starting point to move the time window.

It should be noted that, when calculating the water depth by using the SEGY data, the seafloor searching start time of the profile may be set manually, for example, the seafloor searching start time is set to be 2ms, that is, signals within 2ms are directly ignored, and the amplitude is not calculated, because the data acquisition process is affected by the environment and the equipment during the initial searching, and the interference information is more.

202. And selecting and obtaining the maximum window sample point amplitude sum from the window sample point amplitude sum data set, taking the maximum amplitude in the maximum window sample point amplitude sum as the maximum sample point amplitude, and taking the moment corresponding to the maximum sample point amplitude as the seabed recording moment.

And screening out the maximum value in S1 to Sn to obtain the maximum window sample point amplitude sum Smax, wherein a time window area corresponding to the maximum window sample point amplitude sum Smax is provided with a plurality of sample point amplitudes, the maximum value is the maximum sample point amplitude, and the maximum value corresponds to a specific recording moment, namely the submarine recording moment in the application.

In the process of searching the maximum window sample point amplitude sum, the time window position after each movement and the sample point amplitude sum thereof can be recorded, and after the time window sliding is finished, the maximum window sample point amplitude sum is uniformly screened; or after S1 and S2 are obtained, firstly performing a numerical comparison, and performing the numerical comparison with S3 by taking the larger value in S1 and S2 as the reference of the next numerical comparison, and the like, and after Sn is obtained, completing the final numerical comparison, thus obtaining the maximum window sample point amplitude sum.

It should be noted that the above description of the process of searching for the maximum window sample amplitude sum is only an example, and is not a unique limitation of the present application.

In the embodiment of the application, the submarine recording moment is calculated based on the maximum amplitude sample point ordinal number:

after screening out the maximum sample point amplitude from the sample point amplitude set, obtaining the maximum amplitude sample point ordinal number of the sample point corresponding to the maximum sample point amplitude;

determining the submarine recording moment according to the maximum amplitude sampling point ordinal number, the recording delay value and the data sampling interval; wherein the record delay value and the data sampling interval are obtained by the SEGY data de-encoding.

Specifically, after the recording interval duration is calculated according to the following calculation formula, the submarine recording moment is obtained by combining the submarine searching starting moment and the recording interval duration:

Δt＝N*a+b；

wherein Δt represents the recording interval duration; n represents the maximum amplitude sample point ordinal number; a represents a data sampling interval; b denotes a recording delay value.

In practical applications, the data sampling interval is the time interval between the current sample and the next sample, which is typically small, only 0.010 ms-0.018 ms.

The embodiment of the application provides a screening method of maximum sample amplitude, which searches a sample amplitude set of a current measuring point by moving a preset time window, and extracts the maximum sample amplitude and the corresponding seabed recording moment thereof in a time window area corresponding to the maximum value of the sample amplitude sum based on the principle that the window sample amplitude sum is maximum; and (3) carrying out coarse screening on a large amount of sample amplitude data by using a time window, after finding out the maximum value of the window sample amplitude sum, carrying out secondary screening on the sample amplitude in the time window corresponding to the window sample amplitude sum maximum value, and finally finding out the maximum sample amplitude, wherein compared with the line-by-line point-by-point extraction, the method for screening the maximum sample amplitude by using the sliding time window can save a large amount of data extraction time, and compared with a manual picking method, the accuracy of data extraction is higher.

Example III

In shallow profile measurement, there may be a condition that ping is lost or reflection on the seabed is weak, so that the maximum amplitude searched is not located on the seabed but above or below the seabed, so in order to screen out data jump points caused by ping loss or weak reflection on the seabed, the embodiment of the application provides a track water depth extraction method for removing water depth jump points, which filters the water depth values of the measurement points of each measurement point, deletes the water depth jump points located outside a preset threshold range, stores correct measurement point water depth data, and forms the track water depth with high accuracy.

Fig. 3 is a flow chart of a track water depth extraction method for removing water depth jump points according to an embodiment of the application.

Referring to fig. 3, the track water depth extraction method for removing the water depth jump points comprises the following steps:

301. performing data processing on SEGY data of each measuring point on the track to obtain the coordinates and the water depth of the measuring point of each measuring point;

in the embodiment of the present application, the specific content of step 301 is described in detail in the above-mentioned first embodiment and second embodiment, and will not be described here again.

302. Acquiring a water depth jump point reference range;

in the embodiment of the application, the water depth jump point reference range can be determined by any one parameter of a water depth value filtering threshold value and a preset seabed amplitude magnitude value or determined by the two parameters together.

Exemplary:

when the water depth jump point reference range is determined by the water depth value filtering threshold, firstly acquiring the water depth value filtering threshold, taking a numerical region larger than the water depth value filtering threshold as the water depth jump point reference range, and correspondingly locating data points in the range as position points below the seabed;

when the water depth jump point reference range is determined by a preset seabed amplitude value, firstly acquiring the preset seabed amplitude value, calculating to obtain a 1/4 value of the preset seabed amplitude value, taking a numerical range which is lower than the 1/4 value of the preset seabed amplitude value as the water depth jump point reference range, and correspondingly locating data points in the range as position points above the seabed;

when the water depth jump point reference range is determined by the water depth value filtering threshold value and the preset seabed amplitude value, taking a numerical range which is smaller than 1/4 value of the preset seabed amplitude value and larger than the water depth value filtering threshold value as the water depth jump point reference range.

303. And removing the water depth of the measuring point in the water depth jump point reference range of the measuring point of each measuring point.

Specific:

and when the water depth jump point reference range is determined by the water depth value filtering threshold value, removing the water depths of the measuring points, which are larger than the water depth value filtering threshold value, from the water depths of the measuring points of each measuring point.

When the water depth jump point reference range is determined by a preset seabed amplitude value, eliminating the water depths of the measuring points of which the water depths are smaller than 1/4 value of the preset seabed amplitude value;

when the water depth jump point reference range is jointly determined by the water depth value filtering threshold value and the preset seabed amplitude value, removing the water depths of the measuring points of which the water depths are smaller than 1/4 value of the preset seabed amplitude value and larger than the water depth value filtering threshold value.

In the actual surveying process, under the general condition, when shallow water is measured, the time interval of each data is short, so that the data points are large in number and densely distributed, the occurrence probability of ping loss or weak submarine reflection and the like is low, so that the water depth jump points cannot be too much, the water depth jump points are directly deleted, the surveying project cannot be greatly influenced, interpolation supplement is not needed after the water depth jump points are removed, and the surveying precision cannot be influenced.

The embodiment of the application provides a track water depth extraction method for removing water depth jump points, which is used for removing the water depth jump points in a water depth value filtering mode after the measurement point coordinates and the measurement point water depths of all measurement points are obtained, wherein abnormal points with overlarge amplitudes are removed through a water depth value filtering threshold value and/or abnormal points with overlarge amplitudes are removed through preset seabed amplitude values, so that the situation that the water depth jump points do not exist in the obtained track water depth is ensured, and the accuracy of the track water depth is influenced.

Example IV

The application also provides an electronic device and corresponding embodiments corresponding to the method embodiment of the first embodiment.

Referring to fig. 4, an electronic device 1000 includes a memory 1010 and a processor 1020; the memory 1010 has stored thereon executable code that, when processed by the processor 1020, can cause the processor 1020 to perform some or all of the methods described above.

The processor 1020 may be a central processing unit (Central Processing Unit, CPU), but may also be other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), field programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.

Memory 1010 may include various types of storage units, such as system memory, read Only Memory (ROM), and persistent storage. Where the ROM may store static data or instructions that are required by the processor 1020 or other modules of the computer. The persistent storage may be a readable and writable storage. The persistent storage may be a non-volatile memory device that does not lose stored instructions and data even after the computer is powered down. In some embodiments, the persistent storage device employs a mass storage device (e.g., magnetic or optical disk, flash memory) as the persistent storage device. In other embodiments, the persistent storage may be a removable storage device (e.g., diskette, optical drive). The system memory may be a read-write memory device or a volatile read-write memory device, such as dynamic random access memory. The system memory may store instructions and data that are required by some or all of the processors at runtime. Furthermore, memory 1010 may comprise any combination of computer-readable storage media including various types of semiconductor memory chips (DRAM, SRAM, SDRAM, flash memory, programmable read-only memory), magnetic disks, and/or optical disks may also be employed. In some implementations, memory 1010 may include readable and/or writable removable storage devices such as Compact Discs (CDs), digital versatile discs (e.g., DVD-ROMs, dual-layer DVD-ROMs), blu-ray discs read only, super-density discs, flash memory cards (e.g., SD cards, min SD cards, micro-SD cards, etc.), magnetic floppy disks, and the like. The computer readable storage medium does not contain a carrier wave or an instantaneous electronic signal transmitted by wireless or wired transmission.

The specific manner in which the respective modules perform the operations in the electronic apparatus in the above embodiments has been described in detail in the embodiments related to the method, and will not be described in detail here.

The aspects of the present application have been described in detail hereinabove with reference to the accompanying drawings. In the foregoing embodiments, the descriptions of the embodiments are focused on, and for those portions of one embodiment that are not described in detail, reference may be made to the related descriptions of other embodiments. Those skilled in the art will also appreciate that the acts and modules referred to in the specification are not necessarily required for the present application. In addition, it can be understood that the steps in the method of the embodiment of the present application may be sequentially adjusted, combined and pruned according to actual needs, and the modules in the device of the embodiment of the present application may be combined, divided and pruned according to actual needs.

Furthermore, the method according to the application may also be implemented as a computer program or computer program product comprising computer program code instructions for performing part or all of the steps of the above-described method of the application.

Alternatively, the application may also be embodied as a non-transitory machine-readable storage medium (or computer-readable storage medium, or machine-readable storage medium) having stored thereon executable code (or a computer program, or computer instruction code) which, when executed by a processor of an electronic device (or electronic device, server, etc.), causes the processor to perform part or all of the steps of the above-described method according to the application.

Those of skill would further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the application herein may be implemented as electronic hardware, computer software, or combinations of both.

The flowcharts and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems and methods according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.

The foregoing description of embodiments of the application has been presented for purposes of illustration and description, and is not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the various embodiments described. The terminology used herein was chosen in order to best explain the principles of the embodiments, the practical application, or the improvement of technology in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.

Claims

1. A method for extracting track water depth by using SEGY data, comprising the steps of:

2. The method for extracting track water depth using SEGY data as claimed in claim 1, wherein said screening out the maximum sample amplitude from the set of sample amplitudes comprises:

3. The method for extracting the track water depth by using the SEGY data according to claim 2, wherein the searching the sample amplitude set based on the principle of maximum sum of window sample amplitudes according to the preset time window size to obtain the maximum sample amplitude and the corresponding submarine recording moment comprises the following steps:

4. The method for extracting track water depth by using SEGY data according to claim 1, wherein the calculating the measuring point water depth according to the submarine recording time corresponding to the maximum sample amplitude, the submarine search starting time, the transducer draft and the seawater sound velocity value comprises:

5. The method for extracting the track water depth by using the SEGY data according to claim 1, wherein after obtaining the coordinates and the water depth of each measuring point, the method further comprises:

6. The method for extracting the track water depth by using the SEGY data according to claim 5, wherein the filtering the water depth value of the measuring point of each measuring point comprises:

acquiring a water depth value filtering threshold value;

7. The method for extracting the track water depth by using the SEGY data according to claim 5 or 6, wherein the filtering the water depth value of the measuring point of each measuring point comprises:

acquiring a preset seabed amplitude value;

8. The method for extracting track water depth using SEGY data as claimed in claim 1, wherein after screening out the maximum sample amplitude from the set of sample amplitudes, the method comprises:

9. The method of extracting track water depth using SEGY data as claimed in claim 8, wherein said determining said sub-sea recording time based on said maximum amplitude sample number, recording delay value and data sampling interval comprises:

Δt＝N*a+b；

10. The method for extracting track water depth using SEGY data according to claim 1, comprising, before the coordinate correction based on the GPS antenna head coordinates: