Method of determining the orientation of a surveying instrument in a borehole

Abstract

A method of determining the orientation of a surveying instrument in a borehole from measurements of the earth's gravitational and magnetic fields wherein the azimuth angle is calculated directly without correcting an apparent azimuth. Those measurements are used to calculate the borehole axial magnetic component at each orientation of the surveying instrument in the borehole, and the borehole axial magnetic component is used, along with those measurements, to directly calculate the proper azimuth angle. In a preferred embodiment, measurements are taken at a plurality of orientations in the borehole and the borehole axial magnetic component is calculated in simultaneous equations and then compared to insure the accuracy of the azimuth computation at each orientation. The orientations are preferably non-colinear.

Description

BACKGROUND OF THE INVENTION

The present invention relates to a method of determining the orientation of a surveying instrument in a borehole. More particularly, the present invention relates to a method of computing the azimuth of a surveying instrument at a particular orientation in a borehole from the measured gravitational and magnetic field data of the earth at that orientation without measuring and correcting an apparent azimuth.

A number of instruments have been described for determining the orientation of a borehole. See M. H. Haddock, Deep Borehole Surveys and Problems, McGraw-Hill Publishing (1931). One such instrument presently in use is mounted in the drill collar and utilizes three orthagonally positioned magnetometers and three orthagonally positioned accelerometers to provide voltage outputs which are proportional to the earth's gravitational and magnetic field vector components. See, for instance, U.S. Pat. Nos. 3,791,043, 3,862,499 and 4,163,324.

It is conventional practice for the results of a borehole survey to be expressed in terms of a series of values of an azimuth angle .theta. and an inclination angle .phi. taken along the length of the borehole. Those angles .theta. and .phi. are either measured or calculated from the magnetic and gravitational field data measured by the above-described accelerometers and magnetometers at a series of locations along the borehole. However, this type survey instrument is mounted in the drill string, and the drill string has its own magnetic field. See J. Chadderdon and W. T. Kittinger, "How to get the most from nonmagnetic collars for directional drilling," Oil and Gas Journal, Apr. 16, 1962, pp. 104-109. If azimuth is measured directly, the apparent azimuth must be corrected for the magnetic field of the drill string. If the azimuth angle is calculated, the magnetic field data must be corrected before that data is used to calculate the azimuth angle.

Methods for correcting that magnetic field data or the apparent azimuth angle are disclosed in U.S. Pat. Nos. 3,791,043 and 4,163,324. In brief, those methods involve an iterative calculation which corrects the apparent azimuth or the magnetic field measurements for the magnetic interference of the drill string.

Applicant has discovered that the azimuth angle can be calculated directly from the gravitational and magnetic field data without measuring the apparent azimuth and correcting that apparent azimuth for the magnetic field of the drill string. The result is that the azimuth can be calculated with a higher degree of accuracy and reliability than previously possible. The method of the present invention also lends itself to convenient computerized calculation. It is, therefore, an object of the present invention to provide a method of determining the orientation of a surveying instrument in a borehole which does not require the correction of the magnetic field data to calculate the azimuth angle of the instrument in the borehole.

It is another object of the present invention to provide a method for determining the orientation of a surveying instrument in a borehole which does not require the measurement and correction of an apparent azimuth.

It is another object of the present invention to provide a method of determining the orientation of a surveying instrument in a borehole comprising measuring the earth's gravitational and magnetic fields at a plurality of non-colinear orientations of a surveying instrument in a borehole, calculating the borehole axial magnetic component at each orientation directly from the gravitational and magnetic field measurements and computing the azimuth angle for each orientation directly from the gravitational and magnetic field measurements and the calculated borehole axial magnetic component.

It is another object of the present invention to provide a method of determining the orientation of a surveying instrument in a borehole width increased accuracy and reliability.

It is another object of the present invention to provide a method of determining the orientation of a surveying instrument in a borehole from the earth's measured gravitational field in the x, y and z axes and from the earth's measured magnetic field in the x and y axes.

A further object of the present invention is to provide a method of calculating the azimuth angle of a surveying instrument at a selected orientation in a borehole from the earth's gravitational and magnetic field data measured by that surveying instrument.

It is another object of the present invention to provide a method of calculating the inclination angle, dip and other survey data for a selected orientation of a survey instrument in a borehole from the earth's magnetic and gravitational field data measured by that survey instrument.

These and other objects and advantages will be clear to those skilled in the art who have the benefit of this disclosure from the following detail description of a preferred embodiment of the present invention.

SUMMARY OF THE INVENTION

Those objects and advantages are accomplished through the present invention by providing a method for determining the orientation of a surveying instrument in a borehole comprising measuring the earth's gravitational and magnetic fields at a plurality of non-colinear orientations of the surveying instrument in the borehole, calculating the borehole axial magnetic component at each orientation directly from the gravitational and magnetic field measurements, and computing the azimuth angle for each orientation directly from the gravitational and magnetic field measurements and the calculated borehole axial magnetic component. The azimuth angle is calculated directly from the gravitational and magnetic field measurements without correcting the magnetic field measurements.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

The earth's gravitational and magnetic fields are measured with a surveying instrument or probe, such as that described in U.S. Pat. No. 4,163,324, which produces six outputs in the form of variable voltages proportional to the vector components of the local gravitational and magnetic fields. The voltage output of such a probe can be represented as H.sub.x, H.sub.y, and H.sub.z for the magnetic field measurements of the earth in the x, y and z axes, where the z axis lies along the longitudinal axis of the borehole and the x and y axes lie in a radial plane perpendicular to the longitudinal axis of the borehole, and as G.sub.x, G.sub.y and G.sub.z for the gravitational field measurements of the earth in that same x, y, z coordinate system. H.sub.z is the component of the total magnetic field which is suspect due to the magnetic field of the drill string in which the probe is mounted along the longitudinal borehole axis.

Because the method of the present invention calculates H.sub.z from the measured H.sub.x, H.sub.y, G.sub.x, G.sub.y and G.sub.z data from a selected orientation in the borehole, it is not necessary that the probe produce a voltage output for the earth's magnetic field along the z axis. Although the apparatus disclosed in U.S. Pat. No. 4,163,324 does produce an output corresponding to H.sub.z, that output can simply be ignored. A probe with two orthagonally positioned magnetometers for measurement of the magnetic field in just the x and y axes, and therefore, with just five outputs, is also satisfactory for use with the method of the present invention.

The calculation of the magnetic field along the longitudinal borehole axis directly from the measured variables H.sub.x, H.sub.y, G.sub.x, G.sub.y and G.sub.z is based upon the assumption that (1) the earth's gravitational field is constant and (2) the earth's magnetic field intensity is constant. Having made those assumptions, and based upon the rule that the total field intensity at a given orientation is the square root of the sum of the squares of the magnetic field intensity, the total field intensity at one orientation must equal the total field intensity at a second selected orientation, a relationship which can be expressed as:

H.sub.x.sbsb.1.sup.2 +H.sub.y.sbsb.1.sup.2 +H.sub.z.sbsb.1.sup.2 =H.sub.x.sbsb.2.sup.2 +H.sub.y.sbsb.2.sup.2 +H.sub.z.sbsb.2.sup.2 [ 1] The same rule applies for the total gravitational field, so: H.sub.x.sbsb.1 G.sub.x.sbsb.1 +H.sub.y.sbsb.1 G.sub.y.sbsb.1 +H.sub.z.sbsb.1 G.sub.z.sbsb.1 =H.sub.x.sbsb.2 G.sub.x.sbsb.2 +H.sub.y.sbsb.2 G.sub.y.sbsb.2 +H.sub.z.sbsb.2 G.sub.z.sbsb.2 [ 2]

At this point, a second assumption is made, and the basis for this assumption is the above-referenced article by Chadderdon and Kittinger which appeared in the Apr. 16, 1962 issue of The Oil and Gas Journal to the effect that the magnetic field due to the drill string (which is the variable which affects H.sub.z) is all in the direction of the longitudinal (z) axis of the borehole, i.e., that the source of the error in H.sub.z is entirely along the longitudinal axis of the borehole or drill string having the surveying instrument mounted therein. Solving for H.sub.z.sbsb.1 in equation [2], substituting into [1], and rearranging the resulting equation gives: ##EQU1##

To simplify matters, we will do the following substitutions: ##EQU2##

Substituting a, b and c into equation [5] gives: ##EQU3## Equations [6]-[8] are then substituted back into equation [10] and the longitudinal boreholes axial magnetic component H.sub.z.sbsb.1 can be solved for directly using only the earth's gravitational field data and the earth's magnetic field data in the x and y axes: ##EQU4##

Due to the use of the quadratic equation, H.sub.z.sbsb.1 can have two values at each orientation of the survey instrument in the borehole when solved by equation [11]. That ambiguity can be removed by calculating H.sub.z in simultaneous equations at a plurality of orientations and comparing each H.sub.z to determine whether each H.sub.z is within the realm of possible values. To accomplish that comparison effectively, each H.sub.z must be calculated from an orientation which is non-colinear within the H.sub.z to which it is compared. In a presently preferred embodiment, the probe is provided with sensors of about 0.05% linearity, and the method of the present invention provides an accurate calculation of H.sub.z with as little as 5.degree. of hole turn. It is, of course, possible to calculate H.sub.z directly with only two orientations in situations in which there is no ambiguity.

Having provided a method for calculation of the borehole axial magnetic component, the next step is to calculate the actual azimuth angle using an equation which is derived as follows. The starting point is an x, y, z coordinate system oriented along the earth's actual gravitational and magnetic fields and having vector components G.sub.x, G.sub.y, G.sub.z, H.sub.x, H.sub.y, and H.sub.z, respectively: ##STR1## That coordinate system is rotated from the position .theta., .phi., .alpha. (where .theta. is the azimuth angle, .phi. is the inclination angle, and .alpha. is the gravitational toolface) to a coordinate system in which the z axis is down and the y axis is pointed along the axis of the borehole. "Toolface" refers to the angle between the y axis of the survey instrument and 90.degree. from the radial, measured and looking in the direction of the positive z axis: ##STR2## To accomplish that rotation requires a negative rotation of around the z axis to achieve 0 toolface (i.e., .alpha.=0): ##EQU5## The next step is to rotate that coordinate system around the x axis (which is horizontal as a result of the previous rotation) so that z points down. This rotation requires a positive rotation of .phi.: ##EQU6## The two matrices are then combined: ##EQU7##

The resulting matrix is then used to rotate the data from the earth's magnetic field (H.sub.x, H.sub.y, and H.sub.z) into a coordinate system in which z is down and y is along the borehole axis in the radial plane. That new coordinate system is represented as H.sub.x, H.sub.y, H.sub.z : ##EQU8##

H.sub.x =H.sub.x cos .alpha.-H.sub.y sin .alpha. [12] H.sub.y =H.sub.x sin .alpha. cos .phi.+H.sub.y cos .alpha. cos .phi.+H.sub.z sin .phi. [13] H.sub.z =-H.sub.x sin .alpha. sin .phi.-H.sub.y cos .alpha. sin .phi.+H.sub.z cos .phi. [14] Sin .alpha., cos .alpha., sin .phi. and cos .phi. can be expressed in terms of G.sub.x, G.sub.y and G.sub.z as follows: ##EQU9## Substituting [15]-[18] into equations [12] and [14]: ##EQU10## The azimuth angle can also be expressed as Atn(.sup.H x/H.sub.y), therefore equations [19] and [20] are substituted into that expression as follows to enable computation of the azimuth angle: ##EQU11## That equation is simplied as follows: ##EQU12## G.sub.z is bracketed and given a negative sign because it is common practice to install the G.sub.z accelerometer backwards, partly because gravity is considered an acceleration field pointing upwardly and partly because the G.sub.z accelerometer fits into the probe housing more easily when installed backwardly.

It is also possible to calculate other survey data directly from the measured H.sub.x, H.sub.y, G.sub.x, G.sub.y, G.sub.z, and the calculated H.sub.z. For instance, dip is calculated as follows. From substitution of equations [15]-[18] into equation [14]: ##EQU13## Further, the sum of the squares of the total magnetic field intensities of coordinate system H.sub.x, H.sub.y, H.sub.z and H.sub.x, H.sub.y, H.sub.z must be the same, therefore: ##EQU14## Dip can be expressed as Atn(.sup.H z/H.sub.x.sup.2 +H.sub.y.sup.2), and equation [23] can be substituted into that expression as follows: ##EQU15## Substituting equation [22] into equation [24]: ##EQU16##

In a presently preferred embodiment, the gravitational and magnetic field data is taken as the probe is lowered on a drill string down into the borehole to be surveyed through a series of orientations of the probe in the borehole. The data is stored in memory, the probe is retrieved from the borehole, the data retrieved from the probe and the calculations are run on a separate computer. The probe is also lowered on a wireline and data transmitted to the separate computer from each orientation. These and other variations of the method of the present invention will be apparent to those skilled in the art who have the benefit of this disclosure, and it is expected that all such variations will fall within the spirit and scope of the following claims.

Claims

1. A method of determining the orientation of a surveying instrument in a borehole comprising:

measuring the radial components of the earth's magnetic field at a plurality of orientations of a surveying instrument in a borehole;

measuring the earth's gravitational field at a plurality of orientations of a surveying instrument in a borehole;

calculating the borehole axial magnetic component at each said orientation from said gravitational and the radial magnetic field measurements; and

computing the azimuth angle of each said orientation directly from said gravitational and radial magnetic field measurements and the calculated borehole axial magnetic component.

2. The method of claim 1 wherein the orientations of the surveying instrument in the borehole are non-colinear.

3. The method of claim 1 additionally comprising calculating the dip angle directly from the gravitational and magnetic field measurements and the calculated borehole axial magnetic component.

4. The method of claim 1 additionally comprising comparing the calculated borehole axial magnetic component from a first orientation of the surveying instrument to the calculated borehole axial magnetic component from a second orientation of the surveying instrument.

5. The method of claim 4 wherein the calculated borehole axial magnetic measurement from the first orientation is compared to the calculated borehole axial magnetic component from the second orientation before computing the azimuth angle of the first orientation.

6. The method of claim 1 wherein the measurements of the earth's gravitational field are taken at the same orientation as the measurements of the earth's magnetic field at each of the plurality of orientations of the surveying instrument in the borehole.

7. A method of determining the orientation of a surveying instrument in a borehole comprising:

measuring the earth's gravitational field at a first orientation of a surveying instrument in a borehole;

measuring the radial components of the earth's magnetic field at the first orientation of the surveying instrument;

calculating the magnetic field along the borehole longitudinal axis at the first orientation of the surveying instrument from the measurements of the earth's gravitational field and the measurements of the radial components of the earth's magnetic field at the first orientation of the surveying instrument;

measuring the earth's gravitational field at a second orientation of the surveying instrument in the borehole;

measuring the radial components of the earth's magnetic field at the second orientation of the surveying instrument;

calculating the magnetic field along the borehole longitudinal axis at the second orientation of the surveying instrument from the measurements of the earth's gravitational field and the measurements of the radial components of the earth's magnetic field at the second orientation of the surveying instrument;

comparing the calculated magnetic field along the borehole longitudinal axis at the first orientation of the surveying instrument with the calculated magnetic field along the borehole longitudinal axis at second orientation of the surveying instrument; and

computing the azimuth angle at the first and second orientations of the surveying instrument from the measurements of the earth's gravitational field, the measurements of the radial components of the earth's magnetic field, and the calculated magnetic field along the borehole longitudinal axis at each of the respective first and second orientations of the surveying instrument.

8. The method of claim 7 additionally comprising calculating the dip angle directly from the measurements of the earth's gravitational field, the measurements of the radial components of the earth's magnetic field and the calculated magnetic field along the borehole longitudinal axis at each of the respective first and second orientations of the surveying instrument.

9. The method of claim 7 wherein the first and second orientations of the surveying instrument are non-colinear.

10. A method of determining the orientation of a surveying instrument in a borehole comprising:

measuring the earth's gravitational and magnetic fields at a plurality of orientations of a surveying instrument in a borehole;

calculating the borehole axial magnetic component at each said orientation from the gravitational and magnetic field measurements without measuring the magnetic field along the longitudinal borehole axis; and

computing the azimuth angle of each said orientation directly from said gravitational and magnetic field measurements and the calculated borehole axial magnetic component.

11. The method of claim 10 wherein the orientations of the surveying instrument in the borehole are non-colinear.

12. The method of claim 10 additionally comprising comparing the calculated borehole axial magnetic compound from a first orientation of the surveying instrument to the calculated borehole axial magnetic component from a second orientation of the surveying instrument.

13. The method of claim 12 wherein the calculated borehole axial magnetic component from the first orientation is compared to the calculated borehole axial magnetic component from the second orientation by computing the azimuth angle of the first orientation.

14. The method of claim 10 additionally comprising calculating the dip angle from said gravitational and magnetic field measurements and the calculated borehole axial magnetic component.

15. The method of claim 10 wherein the measurements of the earth's gravitational and magnetic fields are both taken at each of said plurality of orientations in the borehole.