TEMPERATURE COMPENSATED PARABOLIC SKEW ESTIMATION FOR CLOCKS ON AN AUTONOMOUS SEISMIC NODE

Abstract

Disclosed is a method for determining skew measurements for clock errors in an autonomous seismic node. A parabolic fit may be used to estimate the clock drift of an ocean bottom seismic node during node deployment. A temperature and/or frequency trend and a real-time temperature measurement may be used to compute a temperature corrected parabolic trend. The temperature and/or frequency trend may be measured in a laboratory on a node by node basis or it may be a single trend that is suitable for all nodes. The method may include measuring clock skew prior to node deployment and after node recovery, correcting the pre and post deployment skew measurements based on a temperature and/or frequency trend and/or to a constant reference temperature, and/or computing a parabolic trend of the skew measurements of the clock based on the temperature corrected pre and post deployment skew measurements.

Claims

1. A method for modeling clock drift of a seismic node, comprising measuring clock skew prior to node deployment to determine pre-deployment skew measurements; measuring clock skew after node recovery to determine post-recovery skew measurements; correcting the pre-deployment skew measurements and the post-recovery skew measurements to a constant reference temperature; computing a parabolic trend of clock skew measurements based on the temperature corrected pre-deployment and post-recovery skew measurements; and correcting the parabolic trend to a variable experiment temperature.

2. The method of claim 1, further comprising measuring a plurality of temperatures of the node during node deployment.

3. The method of claim 2, wherein the measuring step is performed continuously.

4. The method of claim 1, wherein the variable experiment temperature comprises a plurality of recorded node temperatures during deployment.

5. The method of claim 1, wherein correcting the pre-deployment and post-recovery skew measurements comprises using the following formula:
Temp Adjusted Skew Measurement=Skewmeasured.sub.to.sup.t[TFCF*(TmeasuredTreference)dt].

6. The method of claim 1, wherein correcting the parabolic trend step comprises using the following formula:
TempAdjustedSkewEstimate(t)=ParabolicSkewEstimate(t)+.sub.to.sup.t[TFCF*(TmeasuredTreference)dt].

7. The method of claim 1, further comprising performing each of these steps on a plurality of deployed nodes.

8. The method of claim 1, further comprising performing each of these steps on a single node of a plurality of deployed nodes and using the computed parabolic trend for the plurality of deployed nodes.

9. A method for modeling clock drift of a seismic node, comprising measuring clock skew prior to node deployment to determine pre-deployment skew measurements; measuring clock skew after node recovery to determine post-recovery skew measurements; measuring a plurality of temperatures of the node during node deployment; and computing a parabolic trend of the skew measurements of the clock based on the measured temperatures.

10. The method of claim 9, further comprising correcting the pre-deployment skew measurements and the post-recovery skew measurements to a constant reference temperature.

11. The method of claim 10, further comprising correcting the parabolic trend from the temperature corrected skew measurements.

12. A method for modeling clock drift of a seismic node, comprising determining a parabolic fit of clock skew measurements of a plurality of ocean bottom seismic nodes based on temperature measurements during operation of the ocean bottom seismic nodes.

13. The method of claim 12, further comprising correcting the parabolic fit based on the temperature measurements.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0012] The following drawings form part of the present specification and are included to further demonstrate certain aspects of the present invention. The invention may be better understood by reference to one or more of these drawings in combination with the detailed description of specific embodiments presented herein.

[0013] FIG. 1 illustrates a chart comparing measured skew against a linear estimation of the skew using pre-deployment and post recovery skew measurements.

[0014] FIG. 2 illustrates a chart detailing the impact of temperature dependent frequency changes on the shape of the skew trend.

[0015] FIG. 3 illustrates one embodiment of a temperature corrected parabolic estimation method according to the present disclosure.

DETAILED DESCRIPTION

[0016] Various features and advantageous details are explained more fully with reference to the non limiting embodiments that are illustrated in the accompanying drawings and detailed in the following description. Descriptions of well-known starting materials, processing techniques, components, and equipment are omitted so as not to unnecessarily obscure the invention in detail. It should be understood, however, that the detailed description and the specific examples, while indicating embodiments of the invention, are given by way of illustration only, and not by way of limitation. Various substitutions, modifications, additions, and/or rearrangements within the spirit and/or scope of the underlying inventive concept will become apparent to those skilled in the art from this disclosure. The following detailed description does not limit the invention.

[0017] Reference throughout the specification to one embodiment or an embodiment means that a particular feature, structure, or characteristic described in connection with an embodiment is included in at least one embodiment of the subject matter disclosed. Thus, the appearance of the phrases in one embodiment or in an embodiment in various places throughout the specification is not necessarily referring to the same embodiment. Further, the particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.

[0018] The present disclosure provides a method for determining a parabolic fit for a clock skew of an autonomous seismic node at a variety of temperatures. As mentioned above, prior modeling techniques for skew adjustment do not adjust for temperature and do not use parabolic models. In one embodiment, the disclosure utilizes a temperature and/or frequency trend and a real-time temperature measurement to compute a temperature corrected parabolic trend. The temperature and/or frequency trend may be measured in a laboratory on a node by node basis or it may be a single trend that is suitable for all nodes of the particular survey.

Parabolic Estimations of Clock Drift

[0019] One application for the present disclosure is for autonomous seismic nodes. As is known in the art, such autonomous seismic nodes are self-contained independent seismometers, and each seismic node generally includes one or more seismic sensors (e.g., geophone and hydrophone sensors), a data recording unit, a reference clock for time synchronization, and a power source, such as those nodes and internal components described in more detail at U.S. Pat. Nos. 9,494,700 and 9,523,780, incorporated herein by reference. As described herein, one type of clock used in ocean bottom seismic nodes (such as Applicant's MANTA node) is a regulated crystal oscillator. The oven controlled crystal oscillator (OCXO) design of the clock helps maintain the oscillator at a constant temperature and minimize the impact of temperature changes on the clock frequency.

[0020] As with all seismic acquisition systems, ocean bottom seismic nodes must accurately know the time of each individual recorded amplitude measurement. This is typically done through an atomic clock or crystal oscillator within or coupled to the seismic node. For example, a crystal oscillator is designed to oscillate at a pre-described frequency and provides a relative reference of the time of each sample during the deployment period of the ocean bottom node. Even though efforts are made to minimize the impact of temperature, the clock is still very temperature dependent. For example, the oscillation frequency of the crystal clock is sensitive to temperature changes. Temperature changes affect the clock screw, which makes it difficult to fit accurate trends to the clock for skew estimation purposes (such as a parabolic trend).

[0021] FIG. 1 illustrates a chart comparing measured skew against a linear estimation of the skew using pre-deployment and post recovery skew measurements. As seen from FIG. 1, at a constant temperature and assuming that the aging for a clock is linear, the frequency deviation is linear and the resulting timing skew is parabolic. Empirical testing conducted by the Applicant confirms that the skew is parabolic. Thus, an assumption of a linear fit, as is typically done for conventional skew estimation, does not optimally model the skew. Based on the present disclosure, however, it is possible to accurately estimate the parabolic trend from the pre-deployment and post recovery skew measurements. If the node is maintained at a constant temperature for the duration of the seismic experiment the parabolic estimation using pre-deployment and post recovery measurements may be used to determine skew corrections. Unfortunately, crystal oscillator frequency is temperature dependent and the normal ocean bottom node operational process includes large changes in temperature.

[0022] In one embodiment, the clock drift/skew (i.e., the change in clock error over time) can be modeled by the following formula:

[00001]$\mathrm{Skew}.Math..Math.\left(\mathrm{ms}\right)=\left[\frac{{\mathrm{df}}_0}{f_{\mathrm{nom}}}+\left(T-T_{\mathrm{ref}}\right).Math.\right].Math.t+\frac{1}{2}.Math.A.Math.t^2+C$

[0023] where

[0024] f=frequency

[0025] T=Temperature

[0026] Tref=Reference Temperature (25 C.)

[0027] t=time since skew calibration (e.g., since the start of the test)

[0028] =Temperature Coefficient (d(df/F))/dT)

[0029] A=Aging factor (assumes a constant for simplicity, i.e., the frequency changes linearly with time).

[0030] In deployment of ocean bottom nodes, the nodes and related operational process encounters large changes in the external temperature of the node. For example, as part of the offshore process the pre-deployment and post-recovery skew measurements are performed on the vessel at a very different temperature than the temperature the node is exposed to while deployed on the seafloor and/or in a water column for the duration of the experiment. This temperature/frequency dependence causes a distortion in the skew trend which makes it non-parabolic. Further, the change in oscillation frequency (d/t) temperature makes it impossible to accurately fit a parabolic trend to the pre-deployment and post recovery skew measurements without first compensating the measurements for temperature. The present disclosure provides methods for such temperature compensation.

[0031] FIG. 2 is a chart illustrating the impact of temperature dependent frequency changes on the shape of the skew trend, assuming that the clock is linear (i.e., there is no aging component). FIG. 2 illustrates the skew for a deployment time of a seismic node (i.e., a skew/time curve).First line A shows skew at a constant temperature and second line B shows skew at a variable temperature. At constant temperature (see line A), the parabolic trend computation will estimate a straight line and thus frequency deviation will be constant for the full experiment. At variable temperatures (see line B), the change in temperature changes the gradient of the skew curve and thus makes it impossible to reliably fit a parabolic function to the data. Testing conducted by the Applicant for different clocks confirms this effect of variable temperatures on the skew/time curve. For clocks which have a very stable frequency with temperature, the parabolic trend estimation is very accurate (errors of less than 1 ms). However, when the temperature frequency dependence is higher, the parabolic estimation does not represent the true skew and results in errors which are larger in magnitude than that achieved with a linear correction.

[0032] Thus, a need exists for an improved method and system for modeling clock drift using a parabolic fit during applications where temperature varies during use of the clock.

Temperature Corrected Parabolic Estimation

[0033] In one embodiment, the disclosure utilizes a temperature/frequency trend and a real-time temperature measurement to compute a temperature corrected parabolic trend. The temperature/frequency trend may be measured in a laboratory on a node by node basis or it may be a single trend that is suitable for all nodes for a particular survey.

[0034] In one embodiment, the disclosed temperature corrected parabolic estimation method 300 (see FIG. 3) includes the following steps. First, as illustrated in block 302, method 300 comprises measuring skew immediately prior to node deployment and immediately after node recovery. For example, most nodes are deployed and recovered from a marine surface vessel, and the skew may be measured while the nodes are on board the vessel. Further, step 302 may be performed while the node is connected to an external reference clock. Still further, to accurately estimate the quadratic component of the skew it may be necessary to make multiple skew measurements over multiple minutes. Second, as illustrated in block 304, method 300 comprises measuring the temperature of the node throughout the seismic recording and/or other experiment. Regarding step 304, the temperature measurements may be performed continuously, real time, at various predetermined intervals, and/or estimated based on various temperature models. Third, as illustrated in block 306, method 300 comprises correcting the pre-deployment and post-recovery skew measurements to a constant reference temperature. Regarding step 306, the correction may be performed as if the node had been at a single reference temperature for the duration of the experiment. In one embodiment, the following formula is used for correction in step 306: Temp Adjusted Skew Measurement=Skewmeasured.sub.to.sup.t[TFCF*(TmeasuredTreference)dt], where t=time, T=temperature, and TFCF=temperature-frequency calibration factor. Fourth, as illustrated in block 308, method 300 comprises computing a parabolic trend from the temperature corrected pre-deployment and post recovery skew measurements (from step 306). In one embodiment, the computation of the parabolic trend is performed using well established techniques. Fifth, as illustrated in block 310, method 300 comprises correcting the parabolic trend to the skew at the recording temperature. In one embodiment, the following formula may be used for correcting the parabolic trend: TempAdjustedSkewEstimate(t)=ParabolicSkewEstimate(t)+.sub.to.sup.t[TFCF*(TmeasuredTreference)dt], where t=time, T=temperature and TFCF=temperature-frequency calibration factor. One or more of these steps may be omitted.

[0035] In one embodiment, the temperature/frequency trend could be measured in a laboratory on a node by node basis or it could be a single trend that is suitable for all nodes for the particular seismic survey. In one embodiment, a factory or laboratory may calibrate the temperature dependence of the oscillation frequency of every node resulting in a temperature-frequency calibration factor (TFCF) for each node. This calibration may assume a linear relationship between frequency and temperature or it may assume a more complex relationship. In one embodiment, the TFCF may be archived in a database that is accessible during the node recovery, harvesting, and/or data-download process.

[0036] All of the methods disclosed and claimed herein can be made and executed without undue experimentation in light of the present disclosure. While the apparatus and methods of this invention have been described in terms of preferred embodiments, it will be apparent to those of skill in the art that variations may be applied to the methods and in the steps or in the sequence of steps of the method described herein without departing from the concept, spirit and scope of the invention. In addition, modifications may be made to the disclosed apparatus and components may be eliminated or substituted for the components described herein where the same or similar results would be achieved. All such similar substitutes and modifications apparent to those skilled in the art are deemed to be within the spirit, scope, and concept of the invention.

[0037] Many other variations in the disclosed method are within the scope of the invention. For example, the seismic node may be one deployed on the seabed or in other parts of the ocean or body of water. The clock may be located on an autonomous seismic node or other device. The described process may be utilized on any clock that is operated at a variety of temperatures and/or when the clock skew may change based on a variety of operational and/or temperature environments. The operation and/or environment temperature of the clock may be measured continuously, measured at predetermined intervals, and/or estimated based on various temperature models. It is emphasized that the foregoing embodiments are only examples of the very many different methods and uses that are possible within the scope of the present invention.

[0038] Although the invention(s) is/are described herein with reference to specific embodiments, various modifications and changes can be made without departing from the scope of the present invention(s), as presently set forth in the claims below. Accordingly, the specification and figures are to be regarded in an illustrative rather than a restrictive sense, and all such modifications are intended to be included within the scope of the present invention(s). Any benefits, advantages, or solutions to problems that are described herein with regard to specific embodiments are not intended to be construed as a critical, required, or essential feature or element of any or all the claims.

[0039] Unless stated otherwise, terms such as first and second are used to arbitrarily distinguish between the elements such terms describe. Thus, these terms are not necessarily intended to indicate temporal or other prioritization of such elements. The terms coupled or operably coupled are defined as connected, although not necessarily directly, and not necessarily mechanically. The terms a and an are defined as one or more unless stated otherwise. The terms comprise (and any form of comprise, such as comprises and comprising), have (and any form of have, such as has and having), include (and any form of include, such as includes and including) and contain (and any form of contain, such as contains and containing) are open-ended linking verbs. As a result, a system, device, or apparatus that comprises, has, includes or contains one or more elements possesses those one or more elements but is not limited to possessing only those one or more elements. Similarly, a method or process that comprises, has, includes or contains one or more operations possesses those one or more operations but is not limited to possessing only those one or more operations.