The invention concerns correction of clock drift in seismic nodes, and in particular concerns a method for correction of clock drift in seismic nodes, a seismic node and a seismic node handling system.
Seismic surveying may be performed using a number of different solutions. Example solutions are e.g. seismic streamer cables towed behind a vessel, ocean bottom seismic cables, or autonomous seismic recorders/nodes arranged on the ocean bottom. The autonomous seismic recorders may be individually placed on the ocean bottom by e.g. remotely operated vehicles (ROV), by dropping the autonomous seismic recorders into the sea from a vessel, or alternatively the autonomous seismic nodes may be arranged on the ocean bottom attached to a cable deployed from the vessel. The cable may be a rope or a wire. The autonomous seismic nodes may be recovered attached to the cable, picked up by an ROV or picked up from the sea after floating to the surface.
Autonomous seismic nodes are independent seismic recorders that are able to operate on their own while on the sea floor. Autonomous seismic nodes typically comprises one or more compartments containing sensors (geophones or accelerometers and hydrophones), a data recording device and a power source. Autonomous seismic nodes are thus self-contained and do not need to be connected to anything nor communicate with anything when positioned on the sea floor. When hereafter referring to seismic nodes, these are autonomous seismic nodes as explained above.
Each seismic node contains a slave clock providing a reference time for the various sensor data recorded by the node. Before placement on the seabed the slave clock in each seismic node is calibrated and synchronized with a master clock on the vessel for the same absolute time reference. The slave clocks in the nodes are subject to drift during the time the seismic nodes spend on the seabed. This drift is often referred to as a combination of an ageing effect and other effects coming from temperature change etc. A typical drift for a slave clock in a seismic node may be around 0.5-5 milliseconds per 100 days, but may also be higher than 5 milliseconds.
The clock drift of the slave clock while on the seabed is typically measured when the nodes are onboard the vessel again, but clock drift may also be measured while on seabed from an ROV using some kind of connector, optical communication or other forms of wireless communication. After retrieval of the seismic nodes from the seabed, the slave clock in each seismic node is therefore again compared with the vessel master clock and a clock drift, Δt, of each slave clock is measured and recorded. The clock drift is the deviation from the master clock.
Before processing, the time reference of recorded data in each node is corrected for measured clock drift (Δt). The correction for this assumed clock drift can be done and applied to the seismic data in the node, onboard the vessel or later during processing to ensure that all the seismic data recorded by all the seismic nodes are assigned as near as possible to the same correct time of reference which by definition is the vessel master clock. The correction for clock drift is also used for timing reference in other integrated systems such as navigation computers etc. The systems today synchronize the slave clock(s) in each node before deployment and measures the clock drift in the slave clock(s) in each node after retrieval. These two measurements are then used for correction of the time of reference of the recorded sensor data in that node. The two measurements typically assumes and provide a linear correction function for the clock drift.
The present invention provides a method for correction of clock drift in at least one slave clock in a seismic node. The method comprises performing a number of clock drift measurements of the at least one slave clock in the seismic node. A clock drift correction function of the at least one seismic node slave clock is calculated by curve fitting the number of clock drift measurements to a 2nd order polynomial.
The clock drift measurements are performed both before deployment and after retrieval of the seismic node.
The invention provides a method for correcting for clock drift of at least one slave clock in at least one seismic node after a seismic recording by the at least one seismic node resulting in data sets of recorded seismic sensor data as a function of time. The method comprising obtaining a number of clock drift measurements of the at least one seismic node slave clock, and calculating a clock drift correction function of the at least one seismic node slave clock by curve fitting the number of clock drift measurements to a 2nd order polynomial.
A number of clock drift measurements may include obtaining at least one clock drift measurement of the at least one slave clock before recording of seismic sensor data by the seismic node. The at least one clock drift measurement may comprise synchronizing the at least one slave clock with a master clock. The at least one clock drift measurement may include obtaining a first set of clock drift measurements of the at least one slave clock during a first time-period. Obtaining a number of clock drift measurements may further comprise obtaining a second set of clock drift measurements of the at least one slave clock during a time-period after recording of seismic sensor data by the seismic node. The number of clock drift measurements may be obtained when the seismic node is in a storage position in a seismic node storage system. The number of clock drift measurements may be obtained when the seismic node is in an intermediate measurement position in a seismic node handling system. The method may further comprise correcting a time of reference of the recorded seismic sensor data by the 2nd order polynomial clock drift correction function. Correcting the time of reference of the recorded seismic sensor data in the seismic node may be done before or during downloading of the recorded seismic sensor data from the seismic node and to a central computer. The method may be performed in the at least one seismic node. The method may be performed by a clock drift measurement device.
The invention also provides a seismic node comprising a slave clock, and where the seismic node is configured for performing the method above. The seismic node may further be configured for receiving a master clock signal and obtaining a clock drift between the slave clock and the master clock signal.
The invention also provides a handling system for seismic nodes, the handling system comprising a master clock, wherein the handling system is configured for performing the method above. A clock drift between a slave clock in a seismic node and the master clock may be obtained by a clock drift measurement internally in each seismic node.
The invention also provides a computer program for performing the method above.
The synchronization and clock drift measurements are performed while the seismic nodes are connected to a master clock or measuring device before deployment of the seismic nodes and again connected to the master clock or the measuring device after retrieving the seismic node from the seabed. The seismic nodes may be connected to the master clock or the measuring device in a separate location before docking the seismic nodes in the seismic node storage system. The seismic nodes may also be connected to the master clock or the measuring device when the seismic nodes are in their docking position in a seismic node storage system. In the seismic node storage system, the seismic nodes may be docked and connected to the master clock or the measuring device, a central computer or a data management system over a long time, and even as long as the nodes are in their docking position. The number of clock drift measurements performed before and after the seismic nodes have been on the seabed are thus not limited to shortly before deployment and shortly after retrieving of the seismic nodes and nor limited as to the number of measurements performed or the period over which the measurements are done.
An improved fit of the clock drift measurements to the 2nd order polynomial, and thus improved accuracy for the clock drift correction function, is obtained by increasing the number of clock drift measurements performed before recording and after recording by the seismic nodes. Measuring the clock drift before disconnecting the seismic nodes from the master clock or measuring device and after connecting again the seismic nodes to the master clock or measuring device, and fitting these clock drift measurements to a 2nd order polynomial, reduces the possible maximum timing error of the timing reference of the recorded data. This results in increased accuracy and quality of the seismic survey data obtained with the seismic node and thus an improved resolution of the underground. Such an improved clock drift correcting function may also allow for use of less expensive slave clocks in the seismic nodes with better timing accuracy compared to use of prior art linear correction methods.
Example embodiments of the invention will now be described with reference to the followings drawings, where:
Example embodiments are described with reference to the drawings. The same reference numerals are used for the same or similar features in all the drawings and throughout the description.
At least one clock drift measurement of the seismic node slave clock may be obtained before recording of seismic sensor data by the seismic node. The slave clocks in the seismic nodes are synchronized with a master clock before recording of seismic sensor data by the seismic nodes. The first clock drift measurement may thus be provided by the synchronization of the seismic node slave clock with the master clock. However, two or more clock drift measurements of the seismic node slave clock may be obtained before recording by the seismic node. These two or more clock drift measurements may be provided by the synchronization measurement and additional clock drift measurements. These clock drift measurements may be performed during a first specified time-period. Before recording of seismic sensor data, the seismic node is disconnected from the master clock. After recording of seismic sensor data by the seismic node, the seismic node is again connected to the master clock. At least two clock drift measurements of the seismic node slave clock may be obtaining during a second specified time-period after connecting the seismic node to the master clock.
Clock drift correction functions for a slave clock in a seismic node are illustrated in
The 2nd order polynomial function represents a correction function for the clock drift of the at least one slave clock in the at least one seismic node from the time of the first clock drift measurement and to the time of the final clock drift measurement. The time of reference of the recorded seismic sensor data by that node, may then be corrected by the 2nd order polynomial correction function 32 for the clock drift of the slave clock in that node for better data quality in seismic processing. The real clock drift is shown as curve 36. As can be seen from the example in
The time drift d at any time t of a slave clock are expressed as a 2nd order polynomial:
d(t)=at1+bt+c (eq. 1)
where a, b and c are coefficients.
The coefficients a, b and c may be determined based on three clock drift measurements. The first clock drift measurement may be the synchronization of the seismic node slave clock with the master clock at t=0. This implies that d(t)=0=c for t=0. The 2nd order polynomial 32 may thus be calculated with only one measurement (t0) before deployment where t0 may be the time when synchronization or a separate measurement takes place.
In the example below the coefficients a, b and c are determined based on four clock drift measurements in addition to the synchronization the seismic node slave clock with the master clock. The calculation of the coefficients a and b are based on eq.3 and eq.4 below as for the example with four clock drift measurements, however with only three measurements d(t0) and d(t1) are zero.
Determination of the coefficients a, b and c based on four clock drift measurements. These four clock drift measurements are:
d(t0),d(t1),d(t2),d(t3)
The first clock drift measurement d(t0), and the second clock drift measurement d(t1), are performed before disconnecting the seismic node slave clock from the master clock before recording of seismic sensor data.
The first clock drift measurement, d(t0), may be performed after synchronization of the slave clock. The first clock drift measurement is performed for t0=0. The synchronization of the slave clock takes place at t<0. Synchronization may be performed just before the first clock drift measurement, but may also take place well before t=0, e.g. at t=−1000 s or earlier. Earlier may be minutes, hours or days. Tests show that good accuracy of the 2nd order correction function are achieved with synchronization 2-3 days and 20 days before the first clock drift measurement. The first clock drift measurement, d(t0) may also be the synchronization of the seismic node slave clock with the master clock.
The second clock drift measurement d(t1) is performed after the first clock drift measurement at t>0. After the second clock drift measurement, the seismic node may be disconnected from the master clock. The seismic node may then be deployed to the seabed for seismic surveying. The seismic nodes in a seismic survey may be deployed for a various amount of time before retrieval from the seabed. In the example in
The coefficient c in the 2nd order polynomial is calculated as follows:
C=d(t0) (eq.2)
The coefficient a, is found using the double derivate of the 2nd order polynomial function. (eq.3)
The coefficient b is found by substitution (eq.4)
The interval between the first clock drift measurement and the second clock drift measurement may typically be 10 minutes. However, the time periods may also be shorter or longer. An interval of about 10 minutes has been demonstrated to provide a stable solution for a Chip Scale Atomic Clock (CSAC) and a good estimate for the 2nd order correction function for the clock drift as compared with the true clock drift for a slave clock. The accuracy of the clock drift measurement itself is about 11 ns with today's technical solutions. By stable solution means a measurement situation with less chances of calculating an inaccurate slope of the 2nd order polynomial curve. An inaccurate slope will lead to larger uncertainties in the determination of the coefficient a. A longer time interval between the measurements than 10 minutes may further improve the accuracy of the estimated clock drift as compared to the true clock drift for the slave clock.
The accuracy of the correction function is improved for large time differences between the last clock drift measurement d(t3) and the first clock drift measurement d(t0). The improved accuracy also depends on the aging of the crystal in the slave clock. The accuracy as compared with a prior art linear solution, may be improved by 2 to about 1000 times or more. A maximum error over a 100 day period using the Chip Scale Atomic Clock may be reduced to +/−50 to 100 micro second compared to +/−1000 micro second with a linear correction.
More than four time-drift measurements may also be performed. Increasing the number of measurements will increase the accuracy of the 2nd order curve fit and hence the clock drift correction function. The difference between consecutive clock drift measurements in the first time period may be averaged to improve the accuracy of the difference d(t1)−d(t0). Likewise, the difference between consecutive clock drift measurements in the second time period may be averaged to improve the accuracy of the difference d(t4)−d(t3).
Using an average of the clock drift measurements to improve the result, will have the same effect as increasing the time between t1 and t0 and also between t3 and t2. As an example, 10 clockdrift measurements may be performed during the first time interval (before disconnecting from the master clock). Each of these 10 measurements may also be executed at an interval of about 10 minutes. The first clock drift measurement is then made at t=0, the second at t=10 minutes, the third at t=20 minutes etc. Two clock drift measurements may be performed during the second time interval (after again connecting to the master clock), also with an interval of about 10 minutes between each measurement. This may further improve the accuracy of the estimated clock drift function. Other alternatives further improving the accuracy may be possible such as two clock drift measurements during the first time interval, and then a number of clock drift measurements (e.g. 5, 10 or more) during the second time interval. Also e.g. 5, 10 or more clock drift measurements during the first time interval and the second time interval may also be possible. It may also be envisaged that clock drift may be measured as long as the node is connected to the master clock, both during the first time interval (before disconnecting from the master clock) or during the second time interval (after again connecting to the master clock). Clock drift of the slave clock may be measured intermittently or nearly continuously depending on the master clock signal. The interval and number of clock drift measurements performed to ensure a clock drift correction function with acceptable accuracy may also depend on type of clock oscillator.
A possible example procedure for measuring the clock drift of the slave clock in the seismic node may be as follows:
The seismic nodes are connected to the master clock or measurement device. At t<0 the slave clocks are synchronized to the master clock. At t=0 the measurements of clock drift are initiated. The clock drift measurements may be performed a number of times or as long as the node is connected to the master clock or measurement device. It is also possible to use the synchronization to the master clock as the only clock drift measurement before deployment of the seismic nodes onto the seabed. The seismic nodes are disconnected from the master clock or measurement device and deployed onto the seabed. After recording on the seabed, the seismic nodes are retrieved. Onboard the vessel, the seismic nodes are again connected to the master clock or measurement device and clock drift measurements performed. The clock drift measurements may be transmitted to the data management system (DMS) or the central computer onboard the vessel or stored in a memory of the seismic node or seismic node slave clock.
The clock drift measurements may be obtained when the seismic node is in a storage position or docking position in a seismic node storage system. The final storage positions may be in a node storage system having a node storage area(s), a node storage rack(s) or node storage shelf/shelves. The clock drift measurements may be obtained when the seismic node is not in a storage position in the seismic node storage system. The clock drift measurements may be obtained when the node is e.g. in synchronization position or preparation station before deployment, in an intermediate measurement position before or after deployment, a readout station after deployment or other positions in the seismic node handling system.
The connection to the master clock or measurement device may be a cable based connection or a wireless connection.
The clock drift measurements may be obtained in the at least one seismic node. The master clock signal is received in each seismic node, and clock drift measurements between a slave clock signal and the received master clock signal obtained. The clock drift correction function may also be calculated in the seismic node. The clock drift measurements and or the calculation of the clock drift correction function may be performed by the at least one slave clock in the seismic node, in a processor in the seismic node or in a hardware and/or software component in the seismic node. The clock drift measurement may also be obtained by a clock drift measuring device outside the node. The clock drift correction function may be calculated by the clock drift measuring device or in e.g. a central computer of the seismic node handling system.
The time of reference of the recorded seismic sensor data are corrected by the 2nd order polynomial clock drift correction function. The time of reference of the recorded seismic sensor data in the seismic node may be corrected before or during downloading of the recorded seismic sensor data from the seismic node and to a central computer. The recorded seismic sensor data may be transmitted from the at least one seismic node and to a central computer over a cable, a fibre optic cable or a wireless transmission, where the wireless transmission comprising at least one of optical transmission, radio broadcasting transmission, or RF transmission. The recorded seismic sensor data may be corrected for clock drift by use of the 2nd order clock drift function before transmittal or transmitted in uncorrected form. The recorded seismic sensor data may be corrected for clock drift in the central computer or later corrected for clock drift in other external computing facilities.
When the clock drift measurements are obtained when the nodes are in the storage position, this may allow the nodes to be connected to a clock drift measurement system as long as the nodes are docked. Such a solution opens up for the possibility of performing a large number of clock drift measurements further improving the accuracy of the clock drift correction function achieved through the curve fit of these measurements to the 2nd order polynomial.
The seismic nodes may typically also include sensors, internal memory, processors and batteries as required for operation as a seismic node on the seabed. The sensors may e.g. include at least one of geophones, hydrophones, recorders, compasses, tilt sensors and transponders. The batteries may be non-rechargeable or rechargeable. If rechargeable batteries are used, a built in battery charger may also be included. The seismic node and seismic node handling system also include software and/or firmware for executing the procedures performed by and in the nodes and by and in the seismic node handling system.
A communication line for the master clock signal, a communication line for the clock drift signal and any communication lines between the seismic nodes and a seismic node handling system may be a cable based communication line, a fibre optic based communication line or a wireless transmission. The cable connection may be a single line, an electrical cable, or a fibre optic cable. The wireless transmission may e.g. be optical transmission, radio broadcasting transmission, or RF transmission. The charging of batteries and any power transfer may be a cable based solution or a wireless solution. Wireless power transfer includes e.g. induction.
Some example embodiments implementing the method above are provided below.
Example embodiments of procedures for measuring clock drift Δti are illustrated in
The seismic nodes 1 are connected to the vessel master clock 4 before deployment onto the seabed and after retrieval when the seismic nodes are onboard the vessel. The seismic node's internal slave clock 2 sends its internal slave clock synchronization signal 3 to the vessel master clock 4. The vessel master clock 4 compares the slave clock synchronization signal 3 with its own master clock synchronization signal and provides a measurement of the time difference, Δti called clock drift. The internal slave clock synchronization signal 3 and the master clock synchronization signal 5 may also both be sent to a separate clock drift measuring device 11 (
The example embodiment in
In the embodiments in
The input for receiving the signal from the master clock may be provided in a seismic node connector 13 used for connecting the seismic node to the master clock. The seismic node may be connected to the master clock in a storage position. The storage position may be a docking station. The storage position/docking stations comprises a corresponding docking connector 14 for connecting to the seismic node connector 13. The seismic node connector and the docking connector may be provided with plugs or connecting devices for at least one of electrical connection, connectors for optical fibres or wireless communication. The seismic node connector may be designed as a common single node interface 29 with the docking connector 14. This common single node interface 29 may thus be the common interface with the seismic node handling system 20 onboard the vessel. The seismic node interface 29 may thus also provide an interface for e.g. the synchronization and calibration signals received from the master clock before deployment of the seismic nodes, for the signals received from the master clock after recovery of the seismic node, for data communication between a data management system and the seismic node, for transmission of the measured clock drift, and for power and charging of the chargeable batteries in the nodes.
The seismic node may include a transceiver/transmitter for transmitting the measured clock drift (Δti) to the data management system 25 or a central computer 26. The clock drift may be also be transmitted by the separate comparator device 15 (
An example embodiment of a seismic node handling system 20 for clock drift measurement internally in the seismic nodes, is illustrated in
As each node receives the master clock synchronization signal 5 when docked or connected, the clock drift can be measured as soon as the node is connected or docked upon retrieval. The master clock signal is the same signal as used for calibration and synchronization of the slave clocks in the nodes before deployment.
The measurement of the slave clock drift between the slave clock signal and the master clock signal is performed in the nodes themselves by software and/or hardware in each node. Each node 1 transmits the measured time-drift via the data communication line or a separate clock drift communication line 10 to the data management system (DMS) 25 or the central computer 26 onboard the vessel.
The recorded data in each node are corrected for clock drift. As the measurement of the clock drift is performed in the seismic node, the correction of recorded data may be performed in each node before or during downloading of the recorded data from the node to a central computer. Alternatively, the correction of recorded data may also be performed after the recorded data has been downloaded from the seismic node to the central computer. The recorded data are stored in data storage mediums 27, e.g. tapes, or discs or other machine readable medium.
Having described embodiments of the invention it will be apparent to those skilled in the art that other embodiments incorporating the concepts may be used. These and other examples of the invention illustrated above are intended by way of example only and the actual scope of the invention is to be determined from the following claims.
Number | Date | Country | Kind |
---|---|---|---|
20180319 | Mar 2018 | NO | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/NO2019/050047 | 3/1/2019 | WO |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2019/168411 | 9/6/2019 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
7388541 | Yang | Jun 2008 | B1 |
20090143990 | Kamata | Jun 2009 | A1 |
Number | Date | Country |
---|---|---|
2 071 361 | Jun 2009 | EP |
2 525 240 | Nov 2012 | EP |
20160971 | Oct 2017 | NO |
2008070754 | Jun 2008 | WO |
2017213512 | Dec 2017 | WO |
Entry |
---|
Olofsson, Bjorn, et al, “Ensuring Correct Clock Timing in Ocean Bottom Node Acquisition.” Paper presented at the 2010 SEG Annual Meeting, Denver, Colorado, Oct. 2010 (Year: 2010). |
International Search Report dated Jun. 7, 2019 in corresponding International Application No. PCT/NO2019/050047. |
International Preliminary Report on Patentability dated May 12, 2020 in corresponding International Application No. PCT/NO2019/050047. |
Norwegian Search Report dated Sep. 12, 2018 in corresponding Norwegian Application No. 20180319. |
F. Vernotte et al., “Uncertainties of drift coefficients and extrapolation errors: Application to clock error prediction”, Metrologia, 2001, vol. 38, No. 4, pp. 325-342. |
Number | Date | Country | |
---|---|---|---|
20210041586 A1 | Feb 2021 | US |