Borehole survey system utilizing strapdown inertial navigation
First Claim
1. A signal processing method for use in borehole surveys of the type wherein a probe supplies inertial acceleration signals representative of probe inertial acceleration relative to the three axes of a first coordinate system that is referenced to said probe and angular rate signals representative of angular rotation of the probe about the three axes of the first coordinate system as the probe is moved along a borehole by a cable, said signal processing method comprising the steps of:
- (a) transforming said acceleration signals in said first coordinate system to obtain inertial signals representative of movement of said probe in a second coordinate system that is fixed relative to the earth, said inertial signals in said second coordinate system including probe velocity signals;
(b) generating a signal representative of the amount of cable being fed into the entrance opening of said borehole;
(c) processing said signal representative of the amount of cable being fed into said entrance opening of said borehole to obtain a signal that is representative of the progress of said probe along said borehole and is corrected for cable stretch and misalignment of the probe in the borehole and of the probe relative to the first coordinate system;
(d) transforming said inertial signals representative of movement of said probe in said second coordinate system into inertial signals representative of movement of said probe in said first coordinate system;
(e) combining said signal representative of said progress of said probe along said borehole with said inertial signals representative of movement of said probe in said first coordinate system to obtain error signals;
(f) transforming said error signals into said second coordinate system to obtain error correction signals;
(g) combining said error correction signals with said inertial signals representative of movement of said probe in said second coordinate system to obtain corrected probe velocity signals; and
(h) integrating said corrected probe velocity signals to obtain signals representative of the course of said borehole relative to said second coordinate system.
1 Assignment
0 Petitions
Accused Products
Abstract
Disclosed is a borehole survey system that utilizes strapdown intertial navigation techniques for mapping a borehole while the system probe (10) is continuously moved along a borehole (12) by means of a cable (14) that is wound on a cable reel (16). Signals representative of the acceleration of the probe (10) relative to the three axes of a probe body coordinate system (34)) and signals representative of angular rotation of the probe (10) about the three axes of the probe body coordinate system are processed within a signal processor (24) to obtain signals that represent probe velocity and probe position in a level coordinate system (36) that is fixed in orientation relative to the geographic location of the borehole (12). Precise and continuous surveys are accommodated by correction of the level coordinate probe velocity signals and probe position signals with error correction signals that are based on the difference between inertially derived probe body coordinate position signals representative of the distance traveled by the probe (10 ) along the borehole (12) and a cable length signal that is derived from a cable measurement apparatus (26), which indicates the amount of cable (14) fed into or retrieved from the borehole (12). Error correction also is provided to correct for Coriolis effect, centrifugal acceleration and variations in the earth'"'"'s gravitational field as a function of probe depth.
36 Citations
34 Claims
-
1. A signal processing method for use in borehole surveys of the type wherein a probe supplies inertial acceleration signals representative of probe inertial acceleration relative to the three axes of a first coordinate system that is referenced to said probe and angular rate signals representative of angular rotation of the probe about the three axes of the first coordinate system as the probe is moved along a borehole by a cable, said signal processing method comprising the steps of:
-
(a) transforming said acceleration signals in said first coordinate system to obtain inertial signals representative of movement of said probe in a second coordinate system that is fixed relative to the earth, said inertial signals in said second coordinate system including probe velocity signals; (b) generating a signal representative of the amount of cable being fed into the entrance opening of said borehole; (c) processing said signal representative of the amount of cable being fed into said entrance opening of said borehole to obtain a signal that is representative of the progress of said probe along said borehole and is corrected for cable stretch and misalignment of the probe in the borehole and of the probe relative to the first coordinate system; (d) transforming said inertial signals representative of movement of said probe in said second coordinate system into inertial signals representative of movement of said probe in said first coordinate system; (e) combining said signal representative of said progress of said probe along said borehole with said inertial signals representative of movement of said probe in said first coordinate system to obtain error signals; (f) transforming said error signals into said second coordinate system to obtain error correction signals; (g) combining said error correction signals with said inertial signals representative of movement of said probe in said second coordinate system to obtain corrected probe velocity signals; and (h) integrating said corrected probe velocity signals to obtain signals representative of the course of said borehole relative to said second coordinate system. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23)
-
8. The signal processing method of claim 7 wherein said step of updating said transformation matrix further includes the steps of:
-
(a) generating signals representative of the time rate of change in said transformation matrix; and (b) integrating with respect to time said signals representative of said time rate of change in said transformation matrix.
-
- 9. The signal processing method of claim 8 wherein said step of generating said signals representative of said time rate of change in said transformation matrix mathematically corresponds to
- space="preserve" listing-type="equation">C.sub.B.sup.L =C.sub.B.sup.L (ω
.sub.IB.sup.B)-(ω
.sub.IE.sup.L +ω
.sub.EL.sup.L)C.sub.B.sup.L
where CBL represents said signals representative of said time rate of change in said transformation matrix; CBL represents said transformation matrix; parenthenses enclosing a vector represent a skew symmetric matrix formed from the enclosed vector; and
##EQU21## with ω
Bx, ω
By, and ω
Bz respectively indicating said angular rate signals representative of angular rotation of said probe about the three axes of said first coordinate system. - space="preserve" listing-type="equation">C.sub.B.sup.L =C.sub.B.sup.L (ω
-
-
10. The signal processing method of claim 9 further comprising the step of correcting said transformed acceleration signals representative of probe acceleration in said second coordinate system for Coriolis effect, the effects of centrifugal acceleration and the effect of variation in gravitational field as a function of probe depth, prior to the step of integrating said transformed acceleration signals.
-
11. The signal processing method of claim 10 wherein of a time rate of change in said probe velocity signals that represent said velocity of said probe in said second coordinate system corresponds to the mathematical expression:
-
space="preserve" listing-type="equation">V.sup.L =TC.sub.B.sup.L A.sup.B -(2ω
.sub.IE.sup.L +ω
.sub.EL.sup.L)xV.sup.L -G.sup.Lwhere VL represents said time rate of change in said probe velocity signals; VL represents said probe velocity signals; AB represents said acceleration signals representative of said probe acceleration relative to said three axes of said first coordinate system; and
##EQU22## with gzL denoting acceleration due to gravity for the current depth position of said probe.
-
-
12. The signal processing method of claim 7 wherein said step of generating said earth rate signals further comprises the steps of:
-
(a) scaling said probe position error signal by a third set of predetermined scaling factors; and (b) generating a signal representative of the integral with respect to time of the scaled probe position error signals; and
,(c) combining said signal representative of said integral with said earth rate signal to update said earth rate signal.
-
-
13. The signal processing method of claim 12 wherein said step of updating said transformation matrix further includes the steps of:
-
(a) generating signals representative of the time rate of change in said transformation matrix; and (b) integrating with respect to time said signals representative of said time rate of change in said transformation matrix.
-
- 14. The signal processing method of claim 13 wherein said step of generating said signals representative of said time rate of change in said transformation matrix mathematically corresponds to
- space="preserve" listing-type="equation">C.sub.B.sup.L =C.sub.B.sup.L (ω
.sub.IB.sup.B)-(ω
.sub.IE.sup.L +ω
.sub.EL.sup.L)C.sub.B.sup.L
where CBL represents said signal representative of said time rate of change in said transformation matrix; CBL represents said transformation matrix; parenthenses enclosing a vector represent a skew symmetric matrix formed from the enclosed vector; and
##EQU23## with ω
Bx, ω
By, and ω
Bz respectively indicating said angular rate signals representative of angular rotation of said probe about the three axes of said first coordinate system. - space="preserve" listing-type="equation">C.sub.B.sup.L =C.sub.B.sup.L (ω
-
said step of processing said signal representative of said progress of said probe along said borehole comprising the steps of; (a) detecting whether the magnitude of the difference between said probe velocity signal and said cable feed rate signal exceeds a predetermined value; (b) selecting said cable feed rate signal when said magnitude of said difference between said probe velocity signal and said cable feed rate signal does not exceed said predetermined value; (c) selecting said probe velocity signal when said magnitude of said difference between said probe velocity signal and said cable feed rate signal exceeds said predetermined value; and (d) integrating with respect to time the selected one of said probe velocity signals and cable feed rate signals.
-
(a) generating earth rate signals; (b) generating transport rate signals based on said probe velocity signal representative of the velocity of said probe in said second coordinate system; (c) scaling said transformed probe position error signal by a third set of predetermined scaling factors to obtain tilt error correction signals; and
,(d) combining said earth rate signals, said transport rate signals, said tilt error correction signals and said angular rate signals representative of angular rotation of said probe relative to said first coordinate system.
-
IEE
where ##EQU24## with l representing the latitude of the location of said borehole and α
representing the wander angle associated with said location;and where ##EQU25## with ω
E representing the rotation of the earth.
- space="preserve" listing-type="equation">ω
.sub.EL.sup.L =(U.sub.R.sup.L xv.sup.L)/R
,
-
(a) generating signals representative of the time rate of change in said transformation matrix; and (b) integrating with respect to time said signals representative of said time rate of change in said transformation matrix.
- space="preserve" listing-type="equation">C.sub.B.sup.L =C.sub.B.sup.L (ω
.sub.IB.sup.B)-(ω
.sub.IE.sup.L +ω
.sub.EL.sup.L)C.sub.B.sup.L
##EQU26## with ω
Bx, ω
By, and ω
Bz respectively indicating said angular rate signals representative of angular rotation of said probe about the three axes of said first coordinate system.
-
space="preserve" listing-type="equation">V.sup.L =TC.sub.B.sup.L A.sup.B -(2ω
.sub.IE.sup.L +ω
.sub.EL.sup.L)xV.sup.L -G.sup.Lwhere VL represents said time rate of change in said probe velocity signals; VL represents said probe velocity signals; AB represents said acceleration signals representative of said probe acceleration relative to said three axes of said first coordinate system; and
##EQU27## with gzL denoting acceleration due to gravity for the current depth position of said probe.
-
24. A borehole survey system comprising:
-
a probe configured and arranged for passage along said borehole, said probe including acceleration sensing means for supplying inertial acceleration signals representative of the acceleration of said probe relative to three axes of a first coordinate system of fixed orientation relative to said probe, said probe further including angular rate sensing means for supplying angular rate signals representative of angular rotation of said probe about said axes of said first coordinate system; a cable affixed to said probe for raising and lowering said probe through said borehole; cable control means for paying out and retrieving said cable to lower said probe into and retrieve said probe from said borehole; cable measurement means for supplying a cable measurement signal representative of the amount of cable being paid out and retrieved by said cable control means; and
,signal processing means connected for receiving said acceleration signals and said angular rate signals from said probe and connected for receiving said cable measurement signals from said cable measurement means, said signal processing means providing; (a) means for transforming said acceleration signals supplied by said probe into level coordinate system acceleration signals representing acceleration of said probe with respect to a second coordinate system that is fixed relative to the earth; (b) means for converting said level coordinate system acceleration signals into level coordinate system signals representing movement of said probe in said second coordinate system, said signals representing said movement of said probe including signals representing the velocity of said probe relative to said second coordinate system; (c) means for converting said level coordinate system signals representing movement of said probe into inertially derived signals representative of the movement of said probe in said first coordinate system; (d) means responsive to said cable measurement signals for supplying signals representative of the progress of said probe along said borehole and correcting those signals for cable stretch and misalignment of the probe in the borehole and of the probe relative to the first coordinate system; (e) means for combining said signal representative of progress of said probe along said borehole and said signals representative of movement of said probe in said first coordinate system to provide probe movement error signals; (f) means for transforming said probe movement error signals from said first coordinate system to said second coordinate system to obtain probe movement correction signals; (g) means for combining probe movement correction signals with said signals representing movement of said probe in said second coordinate system to provide corrected probe velocity signals; and (h) means for integrating said corrected probe velocity signals to provide signals representative of the coordinates of said borehole relative to said second coordinate system. - View Dependent Claims (25, 26, 27, 28, 29, 30, 31, 32, 33, 34)
-
33. The borehole survey system of claim 32 further comprising means for correcting said probe acceleration signals representative of probe acceleration in said second coordinate system for Coriolis effect, the effects of centrifugal acceleration and the effect of variation in gravitational field as a function of probe depth.
-
34. The borehole survey system of claim 33 wherein a time rate of change in said level coordinate system velocity signals is characterized by the mathematical expression:
-
space="preserve" listing-type="equation">V.sup.L =TC.sub.B.sup.L A.sup.B -(2ω
.sub.IE.sup.L +ω
.sub.EL.sup.L)xV.sup.L -G.sup.Lwhere VL represents said time rate of change in said velocity signals; VL represents said velocity signals; AB represents said acceleration signals representative of said probe acceleration relative to said three axes of said first coordinate system; and
##EQU31## with gzL denoting acceleration due to gravity for the current depth position of said probe.
-
-
Specification