CN111435171A - Clock drift correction method and device for ocean bottom node seismic data and storage medium - Google Patents
Clock drift correction method and device for ocean bottom node seismic data and storage medium Download PDFInfo
- Publication number
- CN111435171A CN111435171A CN201910034554.6A CN201910034554A CN111435171A CN 111435171 A CN111435171 A CN 111435171A CN 201910034554 A CN201910034554 A CN 201910034554A CN 111435171 A CN111435171 A CN 111435171A
- Authority
- CN
- China
- Prior art keywords
- clock drift
- shot point
- drift correction
- node
- point data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012937 correction Methods 0.000 title claims abstract description 123
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004590 computer program Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 13
- 230000005284 excitation Effects 0.000 claims description 5
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 13
- 230000008859 change Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the specification provides a clock drift correction method, a clock drift correction device and a storage medium for ocean bottom node seismic data, wherein the method comprises the following steps: extracting an annular relation cannon taking a target seabed node as a circle center from a seismic observation system of a target work area; each shot point in the annular relation gun has the same offset distance; determining shot point data of the annular relation gun acquired by the target seabed node; carrying out first arrival homophase axis leveling on the shot point data of the annular relation gun; determining the clock drift category of the submarine node according to the overall trend of shot point data after the first arrival homophase axis is leveled; according to the clock drift type of the submarine node, matching a corresponding clock drift correction model from a preset clock drift correction model set; and performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model. The embodiment of the specification can improve the accuracy of the clock drift correction result of the seismic data of the ocean bottom nodes.
Description
Technical Field
The present disclosure relates to the field of seismic data processing technologies for ocean bottom nodes, and in particular, to a clock drift correction method and apparatus for ocean bottom node seismic data, and a storage medium.
Background
In ocean bottom node seismic exploration, ocean bottom nodes are used as receiving equipment of seismic signals and are generally provided with independent clock systems. The subsea node may be synchronized with the onboard main GPS time on the survey vessel power system prior to release. After the power supply of the subsea node is disconnected, the clock of the subsea node drifts according to a certain trend, so that the clock deviates from the GPS time. After release, the drift amount of the clock of the subsea node gradually increases as the subsea node increases in subsea acquisition time. And after the submarine nodes are recovered to the exploration ship and the power supply is switched on, the submarine nodes and the main GPS time are synchronized again, and finally, the total clock drift amount of the submarine nodes in the acquisition process is obtained.
Therefore, in order to improve the quality of seismic data acquired by the ocean bottom nodes (for convenience of description, seismic data acquired by the ocean bottom nodes is hereinafter simply referred to as ocean bottom node seismic data), it is generally necessary to perform clock drift correction on the ocean bottom node seismic data. And moreover, clock drift correction of the seismic data of the submarine nodes is taken as a primary step of quality control of the seismic data of the submarine nodes, and the clock drift correction has important significance on subsequent quality control of positioning, multi-component rotation and the like of the seismic data of the submarine nodes.
However, the current clock drift correction scheme for the ocean bottom node seismic data has low correction accuracy, so that the time correctness of the ocean bottom node seismic data is difficult to guarantee.
Disclosure of Invention
An object of the embodiments of the present disclosure is to provide a method and an apparatus for correcting clock drift of ocean bottom node seismic data, and a storage medium, so as to improve accuracy of a clock drift correction result of the ocean bottom node seismic data.
In order to achieve the above object, in one aspect, the present specification provides a method for correcting clock drift of seismic data of an ocean bottom node, including:
extracting an annular relation cannon taking a target seabed node as a circle center from a seismic observation system of a target work area; each shot point in the annular relation gun has the same offset distance;
determining shot point data of the annular relation gun acquired by the target seabed node;
carrying out first arrival homophase axis leveling on the shot point data of the annular relation gun;
determining the clock drift category of the submarine node according to the overall trend of shot point data after the first arrival homophase axis is leveled;
according to the clock drift type of the submarine node, matching a corresponding clock drift correction model from a preset clock drift correction model set;
and performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model.
In an embodiment of the present specification, the leveling the shot point data of the circular relation gun by the first arrival homophase axis includes:
and performing dynamic correction processing on the shot point data of the annular relation cannon so as to flatten the first arrival of the shot point data of the annular relation cannon to a homophase axis.
In an embodiment of this specification, before leveling the first arrival in-phase axis of the shot point data of the circular relation gun, the method further includes:
and sequencing the annular relation cannons according to the shot point excitation time sequence.
In an embodiment of the present specification, the set of clock drift correction models includes the following clock drift correction models:
a linear type correction model;
a quadratic curve type correction model;
and (4) a jump type correction model.
In another aspect, an embodiment of the present specification further provides a clock drift correction apparatus for ocean bottom node seismic data, including:
the shot point range determining module is used for extracting annular relation shots with the target seabed nodes as the circle centers from the earthquake observation system of the target work area; each shot point in the annular relation gun has the same offset distance;
the shot point data determining module is used for determining shot point data of the annular relation shots acquired by the target seabed node;
the first arrival homophase axis flattening module is used for flattening the shot point data of the annular relation gun;
the integral trend determining module is used for determining the clock drift category of the submarine node according to the integral trend of shot point data after the first arrival homophase axis is leveled;
the drift type determining module is used for matching a corresponding clock drift correction model from a preset clock drift correction model set according to the clock drift type of the submarine node;
and the clock drift correction module is used for performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model.
In an embodiment of the present specification, the leveling the shot point data of the circular relation gun by the first arrival homophase axis includes:
and performing dynamic correction processing on the shot point data of the annular relation cannon so as to flatten the first arrival of the shot point data of the annular relation cannon to a homophase axis.
In an embodiment of the present specification, the apparatus for correcting clock drift of seismic data at a subsea node further includes:
and the shot point sorting module is used for sorting the circular relational guns according to the shot point excitation time sequence before the first arrival in-phase axis leveling module levels the shot point data of the circular relational guns.
In an embodiment of the present specification, the set of clock drift correction models includes the following clock drift correction models:
a linear type correction model;
a quadratic curve type correction model;
and (4) a jump type correction model.
In another aspect, the present specification further provides a clock drift correction apparatus for ocean bottom node seismic data, including a memory, a processor, and a computer program stored on the memory, where the computer program is executed by the processor to perform the following steps:
extracting an annular relation cannon taking a target seabed node as a circle center from a seismic observation system of a target work area; each shot point in the annular relation gun has the same offset distance;
determining shot point data of the annular relation gun acquired by the target seabed node;
carrying out first arrival homophase axis leveling on the shot point data of the annular relation gun;
determining the clock drift category of the submarine node according to the overall trend of shot point data after the first arrival homophase axis is leveled;
according to the clock drift type of the submarine node, matching a corresponding clock drift correction model from a preset clock drift correction model set;
and performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model.
In another aspect, the embodiments of the present specification further provide a computer storage medium, on which a computer program is stored, and the computer program, when executed by a processor, implements the following steps:
extracting an annular relation cannon taking a target seabed node as a circle center from a seismic observation system of a target work area; each shot point in the annular relation gun has the same offset distance;
determining shot point data of the annular relation gun acquired by the target seabed node;
carrying out first arrival homophase axis leveling on the shot point data of the annular relation gun;
determining the clock drift category of the submarine node according to the overall trend of shot point data after the first arrival homophase axis is leveled;
according to the clock drift type of the submarine node, matching a corresponding clock drift correction model from a preset clock drift correction model set;
and performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model.
According to the technical scheme provided by the embodiment of the specification, in the embodiment of the specification, the annular relation cannons taking the target seabed nodes as the circle centers are extracted from the earthquake observation system of the target work area; and partial shot points in all relation shot lines can be contained in the range of the annular relation shot, namely the shot time of the shot points in the range is distributed in the whole acquisition process of the target seabed node, so that the whole variation trend of clock drift of the seabed node seismic data can be effectively described. In addition, the change of the underground medium can cause the sound wave velocity to change correspondingly with the change of the offset distance, and because the shots in the annular relation shot range have the same offset distance, namely have the same sound wave velocity, the first arrival coaxial axes of the signals can be conveniently leveled to the same time reference based on the same sound wave velocity. Therefore, the embodiments of the present description can obtain a more accurate clock drift correction result of the ocean bottom node seismic data.
Drawings
In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present specification, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort. In the drawings:
FIG. 1 is a flow chart of clock drift correction of ocean bottom node seismic data in some embodiments of the present description;
FIG. 2 is a schematic diagram of a circular relationship cannon in an exemplary embodiment of the present description;
FIG. 3 is a diagram illustrating the flare data of a circular relationship gun after first arrival homophase axis flattening in an exemplary embodiment of the present description;
FIG. 4 is a diagram illustrating a result of clock drift correction performed on the shot data shown in FIG. 3 in an exemplary embodiment of the present disclosure;
FIG. 5 is a block diagram of a clock drift correction apparatus for ocean bottom node seismic data in some embodiments of the present description;
fig. 6 is a block diagram of a clock drift correction apparatus for seismic data at a subsea node according to further embodiments of the present disclosure.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present specification, the technical solutions in the embodiments of the present specification will be clearly and completely described below with reference to the drawings in the embodiments of the present specification, and it is obvious that the described embodiments are only a part of the embodiments of the present specification, and not all of the embodiments. All other embodiments obtained by a person skilled in the art based on the embodiments in the present specification without any inventive step should fall within the scope of protection of the present specification.
Referring to fig. 1, a method for clock drift correction of ocean bottom node seismic data of some embodiments of the present description may include the steps of:
s101, extracting an annular relation cannon taking a target seabed node as a circle center from a seismic observation system of a target work area; and each shot point in the annular relation gun has the same offset distance.
In some embodiments provided herein, the seismic is an artificial exploration seismic. The corresponding target work area is the area range of the appointed exploration earthquake. The seismic observation system is preset with an arrangement mode of relative spatial position relation between shot points and receiving point arrangement in seismic exploration. Seismic surveys of different survey areas may employ different seismic observation systems.
In some embodiments provided herein, a circular relationship gun refers to a set of multiple guns distributed in a circular shape in a target work area, that is, a circular relationship gun is a part of guns selected from the set of guns in the target work area.
In some embodiments provided herein, the purpose of extracting the circular relation cannon with the target seabed node as the center is two points: firstly, partial shot points in all relation shot lines can be contained in the range of the annular relation shot, namely the shot time of the shot points in the range is distributed in the whole acquisition process of the target seabed node, so that the whole variation trend of clock drift of the seabed node seismic data can be effectively described. And secondly, the change of the underground medium can cause the corresponding change of the acoustic velocity along with the change of the offset, and because the shots in the annular relation shot range have the same offset, namely the same acoustic velocity, the first arrival homodromous axis of the signal can be leveled to the same time reference based on the same acoustic velocity when the dynamic correction processing is carried out on the shot seismic data.
In some embodiments provided herein, the inner circle radius and the outer circle radius of the circular relationship cannon can be set according to specific requirements. For example, in one exemplary embodiment, the outer radius of the circular relationship cannon may be determined in accordance with the maximum offset of the target seafloor node. On the basis of the determination of the outer circle radius, the inner circle radius of the annular relation gun can be determined according to the requirement that each gun point in the annular relation gun has the same offset distance. In an exemplary implementation scenario, a circular relationship cannon centered on a target seafloor node may be as shown in fig. 2. In fig. 2, O is the position of the target seabed node, and R1 is the inner circle radius of the ring-shaped relation cannon, and the length thereof is 4400 m; r2 is the outer radius of a circular relationship cannon, and the length of the cannon is 4500 m.
S102, determining shot point data of the annular relation cannons collected by the target seabed nodes.
In some embodiments provided herein, the target subsea node may receive seismic signals from various shots within the target work area. Therefore, after the shot range is determined in the manner shown in step S101, for the target seabed node, shot data of each shot in the circular relation gun is also determined.
S103, carrying out first arrival homophase axis leveling on the shot point data of the annular relation gun.
In some embodiments provided in this specification, the first-arrival in-phase axis flattening may be performed by any suitable technique, and is not limited in this specification. For example, in an exemplary embodiment, the first arrival of the shot point data of the circular relational guns can be leveled to the in-phase axis by performing dynamic correction processing on the shot point data of the circular relational guns, so that the influence of the offset on the shot point data can be eliminated.
In some embodiments provided in this specification, before performing first arrival in-phase axis flattening on the shot point data of the circular relational guns, the circular relational guns may be further sorted according to the shot point excitation time sequence, so as to perform subsequent first arrival in-phase axis flattening processing.
And S104, determining the clock drift type of the submarine node according to the integral trend of the shot point data after the first arrival homophase axis is leveled.
In some embodiments provided herein, the seismic data after the first arrival event is leveled off may more significantly reflect the overall trend of clock drift of the ocean bottom node seismic data. Therefore, after the first arrival event is leveled, the clock drift category of the subsea node can be determined based on the overall trend of the shot data after the first arrival event is leveled. The inventors of the present application found through experiments and studies that: clock drift of ocean bottom node seismic data generally presents one of three clock drift categories of linear trend, quadratic curve trend and jump trend.
In an exemplary embodiment, the overall trend of the shot data after the first arrival event is leveled is shown in FIG. 3. In fig. 3, the abscissa represents the sea-bottom node line number and the geophone number, and the ordinate represents the sea-bottom node seismic trace time, as is clear from fig. 3: the integral trend of shot point data after the first arrival homophase axis is leveled is a jump trend. In that
And S105, matching a corresponding clock drift correction model from a preset clock drift correction model set according to the clock drift category of the submarine node.
In some embodiments provided herein, clock drift in view of ocean bottom node seismic data may generally exhibit three clock drift categories, a linear trend, a quadratic trend, and a transition trend. Therefore, a corresponding clock drift correction model can be set in advance for each clock drift category through research experiments, that is, it is possible to set in advance: a linear-type correction model, a quadratic curve-type correction model, and a jump-type correction model.
The linear type correction model expression may be expressed as: t ═ at
Wherein T represents a linear-type drift correction amount, a represents a linear-type model coefficient, and T represents a ratio of a time length from the start of the seabed node to the current seismic recording time to the total recording time length.
The quadratic curve type correction model expression can be expressed as: t ═ bt2
Wherein T represents a quadratic curve type drift correction amount, b represents a quadratic curve type model coefficient, and T represents a ratio of a time length from the start of the seabed node to the current seismic recording time to a total recording time length.
The hopping correction model expression can be expressed as: t ═ at + c
Wherein T represents a jump type drift correction amount, a represents a linear model coefficient, c represents a clock drift jump variable, and T represents the ratio of the time length from the start of the seabed node to the current seismic recording time to the total recording time.
And S106, performing clock drift correction on shot point data acquired by the target seabed node according to the matched clock drift correction model.
In some embodiments provided in this specification, the performing, according to the matched clock drift correction model, clock drift correction on the shot data collected by the target subsea node is that: and substituting shot point data acquired by the target seabed node into the matched clock drift correction model, and calculating according to the matched clock drift correction model to obtain corrected seismic data.
In an exemplary embodiment, the ocean bottom node seismic data shown in FIG. 3 may be corrected to the ocean bottom node seismic data shown in FIG. 4 based on the above-described embodiments of the present specification. In fig. 4, the abscissa represents the sea-bottom node line number and the wave-detecting point number, and the ordinate represents the sea-bottom node seismic trace time, and it can be seen from fig. 4 that the first arrival co-axial of the corrected seismic data is on the uniform time reference.
Therefore, in the above embodiments of the present specification, since the circular relation cannons using the target seabed nodes as the circle centers are extracted from the seismic observation system of the target work area; and partial shot points in all relation shot lines can be contained in the range of the annular relation shot, namely the shot time of the shot points in the range is distributed in the whole acquisition process of the target seabed node, so that the whole variation trend of clock drift of the seabed node seismic data can be effectively described. In addition, the change of the underground medium can cause the sound wave velocity to change correspondingly with the change of the offset distance, and because the shots in the annular relation shot range have the same offset distance, namely have the same sound wave velocity, the first arrival coaxial axes of the signals can be conveniently leveled to the same time reference based on the same sound wave velocity. Therefore, the embodiments of the present description can obtain a more accurate clock drift correction result of the ocean bottom node seismic data.
Referring to fig. 5, a clock drift correction apparatus for ocean bottom node seismic data according to some embodiments of the present disclosure may include:
the shot point range determining module 51 may be configured to extract an annular relation shot using a target seabed node as a circle center from a seismic observation system of a target work area; each shot point in the annular relation gun has the same offset distance;
a shot point data determining module 52, configured to determine shot point data of the circular relation shots acquired by the target seabed node;
a first arrival homophase axis flattening module 53, configured to flatten the first arrival homophase axis of the shot point data of the circular relation gun;
the overall trend determining module 54 may be configured to determine the clock drift category of the subsea node according to an overall trend of the shot data after the first arrival homophase axis is leveled;
a drift category determining module 55, configured to match a corresponding clock drift correction model from a preset clock drift correction model set according to the clock drift category of the subsea node;
the clock drift correction module 56 may be configured to perform clock drift correction on shot data acquired by the target subsea node according to the matched clock drift correction model.
In some embodiments of the present disclosure, the apparatus for correcting clock drift of seismic data of ocean bottom nodes may further include a shot ordering module (not shown in fig. 5), which may be configured to order the circular relational shots according to a shot firing time order before the first arrival in-phase axis leveling module levels the shot data of the circular relational shots.
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functions of the various elements may be implemented in the same one or more software and/or hardware implementations of the present description.
Referring to fig. 6, a clock drift correction apparatus for ocean bottom node seismic data according to further embodiments of the present disclosure may include a memory, a processor, and a computer program stored on the memory, the computer program when executed by the processor performing the steps of:
extracting an annular relation cannon taking a target seabed node as a circle center from a seismic observation system of a target work area; each shot point in the annular relation gun has the same offset distance;
determining shot point data of the annular relation gun acquired by the target seabed node;
carrying out first arrival homophase axis leveling on the shot point data of the annular relation gun;
determining the clock drift category of the submarine node according to the overall trend of shot point data after the first arrival homophase axis is leveled;
according to the clock drift type of the submarine node, matching a corresponding clock drift correction model from a preset clock drift correction model set;
and performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model.
While the process flows described above include operations that occur in a particular order, it should be appreciated that the processes may include more or less operations that are performed sequentially or in parallel (e.g., using parallel processors or a multi-threaded environment).
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method or apparatus that comprises the element.
As will be appreciated by one skilled in the art, embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, the description may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the description may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
This description may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
The above description is only an example of the present specification, and is not intended to limit the present specification. Various modifications and alterations to this description will become apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present specification should be included in the scope of the claims of the present specification.
Claims (10)
1. A clock drift correction method for ocean bottom node seismic data is characterized by comprising the following steps:
extracting an annular relation cannon taking a target seabed node as a circle center from a seismic observation system of a target work area; each shot point in the annular relation gun has the same offset distance;
determining shot point data of the annular relation gun acquired by the target seabed node;
carrying out first arrival homophase axis leveling on the shot point data of the annular relation gun;
determining the clock drift category of the submarine node according to the overall trend of shot point data after the first arrival homophase axis is leveled;
according to the clock drift type of the submarine node, matching a corresponding clock drift correction model from a preset clock drift correction model set;
and performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model.
2. The method for clock drift correction of ocean bottom node seismic data according to claim 1, wherein said flattening of the shot data of said circular relationship guns by first arrival event comprises:
and performing dynamic correction processing on the shot point data of the annular relation cannon so as to flatten the first arrival of the shot point data of the annular relation cannon to a homophase axis.
3. The method for clock drift correction of ocean bottom node seismic data according to claim 1, wherein prior to the flattening of the first arrival in-phase axis of the shot data of said circular relationship cannons, further comprising:
and sequencing the annular relation cannons according to the shot point excitation time sequence.
4. The method for clock drift correction of ocean bottom node seismic data of claim 1, wherein said set of clock drift correction models comprises the following clock drift correction models:
a linear type correction model;
a quadratic curve type correction model;
and (4) a jump type correction model.
5. A clock drift correction apparatus for ocean bottom nodal seismic data, comprising:
the shot point range determining module is used for extracting annular relation shots with the target seabed nodes as the circle centers from the earthquake observation system of the target work area; each shot point in the annular relation gun has the same offset distance;
the shot point data determining module is used for determining shot point data of the annular relation shots acquired by the target seabed node;
the first arrival homophase axis flattening module is used for flattening the shot point data of the annular relation gun;
the integral trend determining module is used for determining the clock drift category of the submarine node according to the integral trend of shot point data after the first arrival homophase axis is leveled;
the drift type determining module is used for matching a corresponding clock drift correction model from a preset clock drift correction model set according to the clock drift type of the submarine node;
and the clock drift correction module is used for performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model.
6. The apparatus for clock drift correction of ocean bottom nodal seismic data as described in claim 5, wherein said flattening of the shot point data of said circular relationship cannons on a first arrival event comprises:
and performing dynamic correction processing on the shot point data of the annular relation cannon so as to flatten the first arrival of the shot point data of the annular relation cannon to a homophase axis.
7. The apparatus for clock drift correction of ocean bottom node seismic data of claim 5, further comprising:
and the shot point sorting module is used for sorting the circular relational guns according to the shot point excitation time sequence before the first arrival in-phase axis leveling module levels the shot point data of the circular relational guns.
8. The apparatus for clock drift correction of ocean bottom node seismic data of claim 5, wherein said set of clock drift correction models comprises the following clock drift correction models:
a linear type correction model;
a quadratic curve type correction model;
and (4) a jump type correction model.
9. An apparatus for clock drift correction of seismic data at a subsea node, comprising a memory, a processor, and a computer program stored on the memory, wherein the computer program when executed by the processor performs the steps of:
extracting an annular relation cannon taking a target seabed node as a circle center from a seismic observation system of a target work area; each shot point in the annular relation gun has the same offset distance;
determining shot point data of the annular relation gun acquired by the target seabed node;
carrying out first arrival homophase axis leveling on the shot point data of the annular relation gun;
determining the clock drift category of the submarine node according to the overall trend of shot point data after the first arrival homophase axis is leveled;
according to the clock drift type of the submarine node, matching a corresponding clock drift correction model from a preset clock drift correction model set;
and performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model.
10. A computer storage medium having a computer program stored thereon, the computer program, when executed by a processor, performing the steps of:
extracting an annular relation cannon taking a target seabed node as a circle center from a seismic observation system of a target work area; each shot point in the annular relation gun has the same offset distance;
determining shot point data of the annular relation gun acquired by the target seabed node;
carrying out first arrival homophase axis leveling on the shot point data of the annular relation gun;
determining the clock drift category of the submarine node according to the overall trend of shot point data after the first arrival homophase axis is leveled;
according to the clock drift type of the submarine node, matching a corresponding clock drift correction model from a preset clock drift correction model set;
and performing clock drift correction on the shot point data acquired by the target seabed node according to the matched clock drift correction model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910034554.6A CN111435171B (en) | 2019-01-15 | 2019-01-15 | Clock drift correction method and device for ocean bottom node seismic data and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910034554.6A CN111435171B (en) | 2019-01-15 | 2019-01-15 | Clock drift correction method and device for ocean bottom node seismic data and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111435171A true CN111435171A (en) | 2020-07-21 |
| CN111435171B CN111435171B (en) | 2023-02-28 |
Family
ID=71580588
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910034554.6A Active CN111435171B (en) | 2019-01-15 | 2019-01-15 | Clock drift correction method and device for ocean bottom node seismic data and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111435171B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111856581A (en) * | 2020-07-27 | 2020-10-30 | 广州海洋地质调查局 | OBS clock drift correction method and processing terminal |
| CN112462425A (en) * | 2020-10-28 | 2021-03-09 | 中国石油天然气集团有限公司 | Method and device for identifying cross interference source in mixed data of submarine nodes |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1712992A (en) * | 2004-06-25 | 2005-12-28 | 中国石油化工股份有限公司 | Correcting method and system for gun point or demodulator probe in earthquake prospection |
| CN104032309A (en) * | 2013-03-06 | 2014-09-10 | 中国石油天然气股份有限公司 | Clock drift correction method in cathodic protection potential synchronous detection |
| CN104483705A (en) * | 2014-12-31 | 2015-04-01 | 中国石油天然气集团公司 | Three-dimensional residual static correction method |
| US20150168576A1 (en) * | 2012-06-04 | 2015-06-18 | Fairfield Industries Incorporated D/B/A Fairfieldnodal | Seismic data analysis using ocean bottom node data collection |
| CN105259570A (en) * | 2015-11-12 | 2016-01-20 | 中国海洋石油总公司 | Seaborne time shifting earthquake receiving point displacement correction method |
| US20160245945A1 (en) * | 2015-02-24 | 2016-08-25 | Seabed Geosolutions B.V. | Single vessel range navigation and positioning of an ocean bottom seismic node |
| CN106291699A (en) * | 2016-08-26 | 2017-01-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | In seismic prospecting, the method for position is correctly blown out in the inclined detection method of big gun and location |
| US20170017006A1 (en) * | 2014-03-19 | 2017-01-19 | Cgg Services Sa | Node position correction in presence of time distortions |
| CN106855640A (en) * | 2017-02-03 | 2017-06-16 | 中国石油天然气集团公司 | A kind of submarine cable geological data first arrival computational methods and device |
| CN106950601A (en) * | 2017-03-03 | 2017-07-14 | 中国石油天然气集团公司 | Static correcting method and device |
| CN107765298A (en) * | 2017-07-14 | 2018-03-06 | 中国石油天然气集团公司 | A kind of method and device for determining near-surface velocity model |
| CN107783186A (en) * | 2017-09-29 | 2018-03-09 | 中国石油天然气集团公司 | A kind of near-surface velocity modeling method and device |
| AR107076A1 (en) * | 2015-12-18 | 2018-03-21 | Bp Corp North America Inc | SYSTEM AND METHOD FOR CORRECTION OF WATCH DERIVA IN RECEIVERS |
| CN107918156A (en) * | 2017-10-30 | 2018-04-17 | 中国石油天然气集团公司 | Detect the method and device of subsea node collection earthquake data polarity |
| CN108445533A (en) * | 2018-03-08 | 2018-08-24 | 青岛海洋地质研究所 | The method of long offset structure phytal zone rate pattern based on OBS |
| US20180335537A1 (en) * | 2015-11-23 | 2018-11-22 | Westerngeco Llc | Gradient-Based 4D Seabed Acquisition Positioning |
-
2019
- 2019-01-15 CN CN201910034554.6A patent/CN111435171B/en active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1712992A (en) * | 2004-06-25 | 2005-12-28 | 中国石油化工股份有限公司 | Correcting method and system for gun point or demodulator probe in earthquake prospection |
| US20150168576A1 (en) * | 2012-06-04 | 2015-06-18 | Fairfield Industries Incorporated D/B/A Fairfieldnodal | Seismic data analysis using ocean bottom node data collection |
| CN104032309A (en) * | 2013-03-06 | 2014-09-10 | 中国石油天然气股份有限公司 | Clock drift correction method in cathodic protection potential synchronous detection |
| US20170017006A1 (en) * | 2014-03-19 | 2017-01-19 | Cgg Services Sa | Node position correction in presence of time distortions |
| CN104483705A (en) * | 2014-12-31 | 2015-04-01 | 中国石油天然气集团公司 | Three-dimensional residual static correction method |
| US20160245945A1 (en) * | 2015-02-24 | 2016-08-25 | Seabed Geosolutions B.V. | Single vessel range navigation and positioning of an ocean bottom seismic node |
| CN105259570A (en) * | 2015-11-12 | 2016-01-20 | 中国海洋石油总公司 | Seaborne time shifting earthquake receiving point displacement correction method |
| US20180335537A1 (en) * | 2015-11-23 | 2018-11-22 | Westerngeco Llc | Gradient-Based 4D Seabed Acquisition Positioning |
| AR107076A1 (en) * | 2015-12-18 | 2018-03-21 | Bp Corp North America Inc | SYSTEM AND METHOD FOR CORRECTION OF WATCH DERIVA IN RECEIVERS |
| CN106291699A (en) * | 2016-08-26 | 2017-01-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | In seismic prospecting, the method for position is correctly blown out in the inclined detection method of big gun and location |
| CN106855640A (en) * | 2017-02-03 | 2017-06-16 | 中国石油天然气集团公司 | A kind of submarine cable geological data first arrival computational methods and device |
| CN106950601A (en) * | 2017-03-03 | 2017-07-14 | 中国石油天然气集团公司 | Static correcting method and device |
| CN107765298A (en) * | 2017-07-14 | 2018-03-06 | 中国石油天然气集团公司 | A kind of method and device for determining near-surface velocity model |
| CN107783186A (en) * | 2017-09-29 | 2018-03-09 | 中国石油天然气集团公司 | A kind of near-surface velocity modeling method and device |
| CN107918156A (en) * | 2017-10-30 | 2018-04-17 | 中国石油天然气集团公司 | Detect the method and device of subsea node collection earthquake data polarity |
| CN108445533A (en) * | 2018-03-08 | 2018-08-24 | 青岛海洋地质研究所 | The method of long offset structure phytal zone rate pattern based on OBS |
Non-Patent Citations (2)
| Title |
|---|
| 侯昆鹏: "OBN地震采集质控技术", 《CPS/SEG北京2018国际地球物理会议暨展览电子论文集》 * |
| 刘乃儒: "陆地节点仪器在滩浅海地震勘探中的应用研究", 《物探装备》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111856581A (en) * | 2020-07-27 | 2020-10-30 | 广州海洋地质调查局 | OBS clock drift correction method and processing terminal |
| CN111856581B (en) * | 2020-07-27 | 2022-02-22 | 广州海洋地质调查局 | OBS clock drift correction method and processing terminal |
| CN112462425A (en) * | 2020-10-28 | 2021-03-09 | 中国石油天然气集团有限公司 | Method and device for identifying cross interference source in mixed data of submarine nodes |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111435171B (en) | 2023-02-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3055715B1 (en) | System and method for regularizing seismic data | |
| CN111435171B (en) | Clock drift correction method and device for ocean bottom node seismic data and storage medium | |
| US9285494B2 (en) | Method for managing shots in a multi-vessel seismic system | |
| CN111352152B (en) | Rapid quality control method and device for seismic data observation system | |
| CN104199103A (en) | Static correction method and static correction device for converted wave | |
| US20150355362A1 (en) | Magnetic field measurement via streamer cables | |
| CN109696705A (en) | Deep-sea OBS uplink and downlink Seismic imaging data correcting method and device | |
| CN112764102A (en) | Seismic data acquisition method and device | |
| CN107544092B (en) | Data processing method and device in vertical seismic profiling (VSP) acquisition | |
| CN110376644A (en) | Earthquake data before superposition joining method and device | |
| CN103501535A (en) | Method and device for locating node | |
| CN111965729B (en) | Real-time monitoring method, system and device for vibroseis combination center | |
| CN108801249A (en) | A kind of inertial navigation system | |
| CN109581494A (en) | Prestack migration method and device | |
| CN112649875A (en) | Method and device for positioning submarine node geophone | |
| CN111624656B (en) | Node clock drift correction method and device based on near offset seismic signals | |
| CN112241020B (en) | Method and device for determining undersampling rate in sparse seismic data acquisition | |
| CN108761523B (en) | Method and device for optimizing combined excitation among seismic sources | |
| CN109143281B (en) | Reducing frequency search space for global navigation satellite system acquisition | |
| CN112444876A (en) | Seismic data processing method and device | |
| CN105301635B (en) | A kind of construction method and device of marine streamer time-lapse seismic data | |
| RU2684491C2 (en) | Method for automatically assigning wireless seismic acquisition units to topographic positions | |
| CN111913215B (en) | Backscatter amplitude-preserving offset imaging method, apparatus and computer storage medium | |
| CN111624655B (en) | Method and device for determining residual static correction value of first-motion wave | |
| CN114063165B (en) | Three-dimensional seismic data splicing method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |