The invention relates generally to calibrating an accelerometer, such as an accelerometer used in a seismic sensor for a marine surveying application.
Seismic surveying is used for identifying subterranean elements, such as hydrocarbon reservoirs, fresh water aquifers, gas injection reservoirs, and so forth. In performing seismic surveying, seismic sources can be placed above a surface under which is located a subterranean structure. One type of seismic surveying is marine seismic surveying, in which seismic sensors can be towed in a body of water or placed on a sea floor above the subterranean structure.
Seismic sensors are typically calibrated by a manufacturer of the seismic sensors prior to delivering such seismic sensors to customers. However, after some amount of use in the field, the sensitivity of seismic sensors can change over time. Some conventional techniques exist to check for sensitivity of seismic sensors in the field for determining whether such seismic sensors exhibit seriously degraded performance. However, such conventional techniques suffer from lack of accuracy. Moreover, conventional sensitivity checking techniques provide relative, not absolute, sensitivity values of seismic sensors, in which a sensitivity of one seismic sensor is relative to the sensitivity of another seismic sensor. To perform proper recalibration of seismic sensors, operators typically have to send the seismic sensors back to the manufacturer, which is costly and time-consuming.
In general, according to an embodiment, a method of calibrating an accelerometer includes rotating a carrier structure that carries the accelerometer, and receiving signals measured by the accelerometer as the carrier structure is rotated. At least one calibration parameter is computed according to the received signals, where the at least one calibration parameter is for use in calibrating the accelerometer.
Other or alternative features will become apparent from the following description, from the drawings, and from the claims.
In the following description, numerous details are set forth to provide an understanding of the present invention. However, it will be understood by those skilled in the art that the present invention may be practiced without these details and that numerous variations or modifications from the described embodiments are possible.
The electronic devices 103 can include sensors, steering or navigation devices, air gun controllers (or other signal source controllers), positioning devices, and so forth. Also depicted in
Although the sources 105 are depicted as being separate from the streamer 102, the sources 105 can also be part of the streamer 102 in a different implementation.
The sensors of the streamer 102 can be seismic sensors, such as hydrophones and/or geophones. Also, the sensors can include inclinometers. The signal sources 105 can be seismic sources, such as air guns or explosives. In some implementations, the geophones can be provided on the streamer 102 along with hydrophones, which are instruments for measuring sound received through water. Providing geophones in addition to hydrophones can be beneficial for various reasons, such as for deghosting and other applications. Note, however, in other implementations, hydrophones do not have to be employed on the streamer 102.
The sensors 202 and 204 can be geophones, and the sensor 206 can be an inclinometer. In some embodiments, each of the geophones 202 and 204 and inclinometer 206 can be implemented with accelerometers, where an accelerometer is an instrument for measuring acceleration such that vibrations can be detected and measured (e.g., as vibrations due to reflected signals from the subterranean structure 114 due to seismic source signals produced by the seismic sources 105). In one specific example, an accelerometer can include a microelectromechanical systems (MEMS) sensor, which is a sensor made using microelectronics in combination with micromachining technology. An MEMS sensor has a generally flat amplitude response versus frequency, down to DC. In other implementations, other types of accelerometers can be used.
In some embodiments, each accelerometer can record a component of a gravity field (expressed in terms of a g unit) along an axis of the accelerometer. Thus, in the example of
In some implementations, the z sensor and y sensor are used to detect seismic signals reflected from the subterranean structure, whereas the inclinometer implemented with the x sensor is used to detect an inclination of the streamer 102 with respect to the x axis.
In certain applications, it is desirable to measure the vertical component (along the z axis) of acceleration (or velocity). To be able to determine such vertical component, the orientations of the accelerometers on the streamer 102 are determined. To do so, two accelerometers (202, 204) can be mounted orthogonally in a plane perpendicular to the streamer axis (x). This plane is defined by the z and y axes. Based on the signals from the orthogonally-mounted accelerometers, the orientations of the accelerometers can be determined such that measurement data from these accelerometers can be rotated according to the orientations of the accelerometers. In another example implementation, instead of using two orthogonally mounted accelerometers, a gimbal mechanism can be used to orient an accelerometer on the streamer 102.
The streamer 102 is rotatable (with the help of the steering devices 208) about the longitudinal axis (x) in a rotational direction indicated by r in
Note also that active steering (such as with the steering devices 208) does not have to be provided to perform the calibration technique according to some embodiments. As a streamer is towed in a body of water, it will rotate over time, such as a result of varying tension in the streamer. The rotation of the streamer depends on the speed of the sea vessel and sea currents. The streamer angle will also change from deployment to deployment. By monitoring the accelerometer values over time, it is possible to find the maximum and minimum amplitude values corresponding to +1 g and −1 g as discussed.
Based on the measured signals produced by an accelerometer on the streamer 102 during rotation of the streamer 102, calibration parameters associated with the accelerometer can be calculated, where the calibration parameters are used for calibrating the accelerometer. Calibrating the accelerometer using the calibration parameters refers to applying the calibration parameters to measurement data from the accelerometer to compensate for known errors of the accelerometer. As discussed in further detail below, the calibration parameters that are computed based on signals received from an accelerometer during rotation of the streamer 102 include a sensitivity of the accelerometer and an offset of the accelerometer.
Moreover, note that the streamer is stored on a reel 104 on the vessel. The reel when unwinding can be used to cause rotation of the accelerometer (as the accelerometer is rotated with the reel as the reel is unwinding). This can be used mainly to calibrate the x sensor 206, although calibration of the DC offset of the y and z sensors can be performed using this technique.
Signals recorded by an accelerometer can be classified into three categories. A first category includes the component of gravity recorded along the axis of the accelerometer, arcsin(φ)*g, φ being the angle of the accelerometer axis with respect to horizontal. The component of gravity will be +1 g or −1 g when the accelerometer is vertical (the axis of the accelerometer is along the z direction) and 0 g when the accelerometer is horizontal (the axis of the accelerometer is in the y or x direction).
A second category of signals that can be recorded by an accelerometer includes noise, which can be acoustic, water-borne noise or streamer vibration noise. The root-mean-square (RMS) noise level of an accelerometer is usually in a range from a few milli-g's (mg) to a few tens of mg in the worst case.
A third category of signals that can be recorded by an accelerometer includes a seismic signal, which refers to seismic energy resulting from a signal produced by the seismic source. The amplitude of a seismic signal can reach a few hundred mg for the first pulse arrival at a sensor located close to the seismic source, and decay to a few mg or less with time and offset from the seismic source.
A cross-sectional view of the streamer along section 3-3 in
The maximum voltage associated with the curve 302 (representing the signal of the z accelerometer) is V1, whereas the minimum voltage of the curve 302 is V2. The voltage V1 is the amplitude of the signal recorded by the z accelerometer for a +1 g measurement, whereas the voltage V2 is the voltage measured by the z accelerometer for a −1 g measurement. The curve 304 representing the signal for they accelerometer similarly has a maximum voltage value and a minimum voltage value that corresponds to the +1 g and −1 g measurements. Maximum voltage V1 and minimum voltage V2 corresponds to +1 g and −1 g assumes that the streamer 102 is horizontal (in other words, the length of the streamer is horizontal with respect to the x axis).
Calibration parameters, including the sensitivity (S0) and offset (Offset) can be calculated as follows:
S0=(V1−V2)/2, (Eq. 1)
Offset=(V1+V2)/2. (Eq. 2)
The sensitivity S0 of an accelerometer indicates the volts-per-g sensitivity of the accelerometer. The Offset parameter represents the DC offset of the accelerometer (the voltage measured by the accelerometer when the accelerometer's axis is horizontal). For an accelerometer that has a flat amplitude response versus frequency down to DC, the calibration parameters computed using the gravity field is also valid in the seismic frequency band. Also, these calibration parameters are absolute calibration parameters for each individual accelerometer, not relative calibration parameters.
The measured V1 and V2 values will be affected by the recorded noise and seismic signal (if any). The amplitude of the first arrival of seismic energy from the seismic source can be up to a few hundred mg's, and can cause a large error when measuring +1 g and −1 g. Therefore, in some embodiments, seismic sources are not activated during calibration of an accelerometer as discussed above.
Moreover, the RMS level of noise is usually a few mg. An error of up to 5 mg when measuring +1 g and −1 g (corresponding to V1 and V2, respectively) can in some cases cause an error in the sensitivity (S0) of about 0.5%, which provides relatively good calibration accuracy.
If improved calibration accuracy is desired, several options can be used. The streamer can be fully rotated several times (rotated a full 360° several times), with V1 and V2 measured in each rotation (of 360°). The various V1 and V2 values due to multiple full rotations of the streamer can be averaged to produce an average V1 and average V2 for computing the S0 and Offset values according to Eqs. 1 and 2. Averaging V1 and V2 values across multiple rotations of the streamer can reduce error in the calibration parameters.
A second option is to slowly rotate the streamer, such that most of the vibration and ambient noise can be filtered out using a high-cut filter with a corner frequency of a few hertz (Hz), which will reduce the noise level without affecting the recording of the gravity component.
If the streamer (or the portion of the streamer containing the accelerometer to be calibrated) is not perfectly horizontal along the x axis (in other words, the length of the streamer (or streamer portion) is not horizontal with respect to the x axis), the maximum gravity field seen by the z and y accelerometers will not be exactly 1 g, but can be slightly less. For example, if the streamer's slope is 3°, the accelerometer will be at 3° from vertical incidents, and a component of gravity measured (V1) is 0.9986 g, where the 0.9986 g is equal to cos(3°).
The error on the +1 g recording in the above example is thus 1−cos(3°)=1−0.9986=1.4 mg, or less than 0.15%. This is also a relatively small error (assuming 3° streamer slope).
An incline in the streamer with respect to the x axis can be corrected if an inclinometer is used, such as an inclinometer implemented with the x sensor 206 of
Moreover, the fact that y and z are perpendicular can also be used to help find the maximum and minimum values, as they will correspond to a 0 g recording on the other axis. Because y and z are perpendicular, it is possible to calculate what absolute DC acceleration one of the y and z accelerometers should record when the acceleration of the other is known. This can be used to perform quality control of the calibration of one axis with measurements from the other axis to help find the maximum and minimum values for one axis, as they will correspond to a 0 g recording on the other axis. Alternatively, both the y and z axes can be calibrated together by combining the measurements from the y and z accelerometers, since the relationship between the y and z measurements allow for recovery of calibration parameters for both axes from several measurements at different positions, even if they do not correspond exactly to the minimum and maximum values.
In another embodiment of the invention, the sensor calibration problem is formulated as a geometrical curve fitting problem. For instance, it is observed that the tip of the vector, whose coordinates are the DC values of the x, y and z components of the acquired measurement, will trace an ellipsoid in the three dimensional coordinate system. Hence, the problem of sensor calibration can be thought of as the problem of fitting an ellipsoid to the measurement.
For the purpose of notational simplicity, it is assumed that the cable is perfectly balanced along the x axis, and there is no cross-talk between x component and the other components. Hence, the DC portion of the other two components will lie on an ellipse at all times in the y-z plane. Algebraically, the DC portion of the particle motion components should satisfy the following equation for an ellipse in the absence of noise:
a22(Yi−Yb)2+2a23(Yi−Yb)(Zi−Zb)+a33(Zi−Zb)2=g2 (Eq 3)
In this equation, Yi, Zi are the DC measurement corresponding to the two components of the particle motion sensors acquired at time i; Yb, Zb are the corresponding biases in the DC measurement; a22, a33 are scalars representing the correction on the actual sensor sensitivities; a23 is a scalar representing the crosstalk between the components of the particle motion sensors; and g=9.8 m/s2. It is to be noted that, in the absence of any bias on the DC measurement, i.e. Yb=Zb=0, the center of the ellipse is the center of the y-z plane. Furthermore, in the absence of any cross talk between individual components, i.e., y and z components are perfectly orthogonal, then a23=0, and the major axes of the ellipse are aligned with the y, and z axes of the coordinate system.
Hence it is noted that, by estimating the parameters of an ellipse that best describes the DC portion of the measurements, the calibration values are obtained for the sensors. Expanding Eq. 3 gives us a slightly different parameterization for the curve that describes the DC measurement:
a22Yi2+2a23YiZi+a33Zi2+b2Yi+b3Zi+c=0 (Eq. 4)
where b2, b3 and c are functions of a22, a23, a33 and g. Since the same equation should be satisfied by all measurements, the following relationship is derived:
where N is the index of the last DC measurement. Eq 4 describes a set of linear equations for the sensitivity corrections and the DC biases. Usually the number of measurements will far exceed the number of unknowns, i.e., N>>6. Additionally, in the presence of noise, the matrix product Md will not be exactly zero. In these cases, the unknown parameters d can be solved by a constrained minimization of a norm of Md:
{tilde over (d)}=arg min∥Md∥, (Eq. 6)
where “∥·∥” is a suitable norm, e.g., the L2 norm, and the operator “argmin” stands for “argument of the minimizer”. Since the minimization problem as stated in Eq 6 has a trivial solution of d=0, the minimization problem can be solved under some constraints. As discussed in J. M. Varah, Least Squares Data Fitting with Implicit Functions, BIT, 36, pp. 842-854 (1996), these constraints may include
It is to be noted that this way of fitting an ellipse to the data is known as the algebraic method. There exist alternative ways of fitting an ellipse to the data, such as the geometric method, the total least squares method and non-linear optimization method. More information on these alternative methods can be found in the references below:
As an illustration,
The accuracy of the described method will improve especially when long records of the DC measurement are available at various streamer orientations. In this respect, the steering of the cable to either side in the cross line directions and some measurements acquired during straight tow with little steering will provide sufficient data. It is also to be noted that, this type of a calibration can be performed once before start of the acquisition and if desired, the estimated calibration values can be continuously improved by using the measurements during seismic acquisition.
In some implementations, moving coil geophones can be used that do not measure gravity field. In such applications, inclinometers may be used in addition to the moving coil geophones to measure the orientation of the moving coil geophones to enable a processing system to rotate the data of the moving coil geophones in vertical and horizontal components. The inclinometers can be accelerometers that measure gravity, but which are not sensitive enough to be used to measure seismic signals. The calibration technique discussed above can also be applied to such inclinometers.
The process of
The CPU(s) 402 is (are) connected to a storage 404 and a communications interface 405 to communicate to a remote network. The storage 404 contains measurement data 406 (which includes data from various sensors of the streamer) as well as calibration parameters 408 calculated by the calibration software 400. The calibration parameters 408 can be communicated through the communications interface 405 to a remote device, such as over a data network.
Instructions of the calibration software 400 are loaded for execution on a processor (such as the one or more CPUs 402). The processor includes microprocessors, microcontrollers, processor modules or subsystems (including one or more microprocessors or microcontrollers), or other control or computing devices. A “processor” can refer to a single component or to plural components.
Data and instructions (of the software) are stored in respective storage devices, which are implemented as one or more computer-readable or computer-usable storage media. The storage media include different forms of memory including semiconductor memory devices such as dynamic or static random access memories (DRAMs or SRAMs), erasable and programmable read-only memories (EPROMs), electrically erasable and programmable read-only memories (EEPROMs) and flash memories; magnetic disks such as fixed, floppy and removable disks; other magnetic media including tape; and optical media such as compact disks (CDs) or digital video disks (DVDs).
While the invention has been disclosed with respect to a limited number of embodiments, those skilled in the art, having the benefit of this disclosure, will appreciate numerous modifications and variations therefrom. It is intended that the appended claims cover such modifications and variations as fall within the true spirit and scope of the invention.
This claims the benefit under 35 U.S.C. §119(e) to U.S. Provisional Application Ser. No. 60/968,496, filed Aug. 28, 2007, which is hereby incorporated by reference.
Number | Name | Date | Kind |
---|---|---|---|
4300220 | Goff et al. | Nov 1981 | A |
4520467 | Berni | May 1985 | A |
5049795 | Moulds, III | Sep 1991 | A |
5479161 | Keyes et al. | Dec 1995 | A |
5621699 | Rigsby et al. | Apr 1997 | A |
5852242 | Devolk et al. | Dec 1998 | A |
6023960 | Abrams et al. | Feb 2000 | A |
6035694 | Dupuie et al. | Mar 2000 | A |
6101864 | Abrams et al. | Aug 2000 | A |
6301195 | Faber | Oct 2001 | B1 |
6353577 | Orban et al. | Mar 2002 | B1 |
6430105 | Stephen | Aug 2002 | B1 |
6512980 | Barr | Jan 2003 | B1 |
6671223 | Bittleston | Dec 2003 | B2 |
6684160 | Ozbek et al. | Jan 2004 | B1 |
6725164 | Bednar | Apr 2004 | B1 |
6758080 | Ragan et al. | Jul 2004 | B1 |
6814179 | Corrigan et al. | Nov 2004 | B2 |
6861587 | SelVakumar et al. | Mar 2005 | B1 |
6871544 | Selvakumar et al. | Mar 2005 | B1 |
6883638 | Maxwell et al. | Apr 2005 | B1 |
6945110 | Selvakumar et al. | Sep 2005 | B2 |
7012853 | Iseli et al. | Mar 2006 | B2 |
7232701 | Gogoi et al. | Jun 2007 | B2 |
7236279 | Yu et al. | Jun 2007 | B2 |
7273762 | Gogoi | Sep 2007 | B2 |
7274079 | Selvakumar et al. | Sep 2007 | B2 |
7292504 | Luc | Nov 2007 | B2 |
7379386 | Muyzert et al. | May 2008 | B2 |
7729202 | Eperjesi et al. | Jun 2010 | B2 |
7734838 | Pavel et al. | Jun 2010 | B2 |
20030088372 | Caulfield | May 2003 | A1 |
20040073373 | Wilson | Apr 2004 | A1 |
20050119831 | Ferber | Jun 2005 | A1 |
20070286023 | Bull et al. | Dec 2007 | A1 |
20080021658 | Pavel et al. | Jan 2008 | A1 |
20080062815 | Iseli | Mar 2008 | A1 |
20080275649 | Ozdemir et al. | Nov 2008 | A1 |
20080316859 | Welker et al. | Dec 2008 | A1 |
20090022009 | Ozdemir et al. | Jan 2009 | A1 |
20100002541 | Ozdemir et al. | Jan 2010 | A1 |
20100116054 | Paulson | May 2010 | A1 |
Number | Date | Country |
---|---|---|
1 192 419 | Apr 2007 | EP |
2 086 055 | May 1982 | GB |
Number | Date | Country | |
---|---|---|---|
20090056411 A1 | Mar 2009 | US |
Number | Date | Country | |
---|---|---|---|
60968496 | Aug 2007 | US |