Gravity gradient measurement method and apparatus

Abstract

Provided is a gravity gradient measurement apparatus and measuring method, wherein a turntable rotates horizontally around an earth-vertical axis, a vacuum layer is arranged on the turntable defining a first chamber, a first three-axis accelerometer and a second three-axis accelerometer are located in the first chamber, the first three-axis accelerometer and the second three-axis accelerometer are arranged symmetrically on an x axis with respect to an origin of coordinates. Both the first three-axis accelerometer and the second three-axis accelerometer have a distance of R from the origin of coordinates. The first three-axis accelerometer and the second three-axis accelerometer are arranged symmetrically on an z axis with respect to the origin of coordinates, and the first three-axis accelerometer and the second three-axis accelerometer are spaced at a distance of h on the z axis. The measurement module uses measurements of the accelerometers to determine gravity gradients on the coordinate axes.

Claims

1. A gravity gradient measurement apparatus, comprising a turntable (3), rotating horizontally around an earth-vertical axis, a first three-axis accelerometer (1), a second three-axis accelerometer (2), a vacuum layer (7), arranged on the turntable and defining a first chamber with the turntable, and a measurement module; wherein, the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are located in the first chamber; the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are arranged symmetrically on an x axis with respect to an origin of coordinates; both the first three-axis accelerometer (1) and the second three-axis accelerometer (2) have a distance of R from the origin of coordinates; the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are arranged symmetrically on a z axis with respect to the origin of coordinates; and the first three-axis accelerometer and the second three-axis accelerometer are spaced at a distance of h on the z axis; and the measurement module determines a gravity gradient on the coordinate axes according to acceleration values measured by the first three-axis accelerometer (1) and the second three-axis accelerometer (2); wherein a first magnetic shield (4) is arranged in the first chamber for shielding the interference of the external geomagnetic field, the first magnetic shield (4) and the turntable (3) define a second chamber, in which the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are located; wherein a second magnetic shield (5) is arranged in the second chamber for shielding the interference of an external alternating current magnetic field, the second magnetic shield (5) and the turntable (3) define a third chamber, in which the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are located.

2. The gravity gradient measurement apparatus of claim 1, wherein, the first magnetic shield (4) is a metal shield layer.

3. The gravity gradient measurement apparatus of claim 2, wherein, the material of the first magnetic shield (4) is permalloy.

4. The gravity gradient measurement apparatus of claim 1, further comprising: a liquid nitrogen layer (6), arranged between the first magnetic shield (4) and the second magnetic shield (5).

5. The gravity gradient measurement apparatus of claim 4, wherein, the second magnetic shield (5) is a superconducting shield layer.

6. The gravity gradient measurement apparatus of claim 5, wherein, the material of the second magnetic shield (5) is a YBCO metallic compound.

7. The gravity gradient measurement apparatus of claim 1, wherein, the first three-axis accelerometer (1) and the second three-axis accelerometer (2) respectively comprises a suspended mass (8), an x-axis acceleration detection electrode (9), a y-axis acceleration detection electrode (10) and a z-axis acceleration detection electrode (11); the x-axis acceleration detection electrode (9), the y-axis acceleration detection electrode (10) and the z-axis acceleration detection electrode (11) define a cavity; and the suspended mass (8) is suspended at a central position of the cavity; the x-axis acceleration detection electrode (9), the y-axis acceleration detection electrode (10) and the z-axis acceleration detection electrode (11) are configured to measure a displacement of the suspended mass (8) in the directions of the x axis, the y axis and the z axis, respectively, so as to obtain the acceleration in the directions of the x axis, the y axis and the z axis according to the measured displacement.

8. The gravity gradient measurement apparatus of claim 7, wherein, the suspended mass (8) is a hollow metal mass with an enclosed outer surface.

9. The gravity gradient measurement apparatus in accordance with claim 1, wherein, the gravity gradient Γ.sup.g is: ${\Gamma }^g=\left[\begin{array}{ccc}{\Gamma }_{\mathrm{XX}}& {\Gamma }_{\mathrm{XY}}& {\Gamma }_{\mathrm{XZ}}\\ {\Gamma }_{\mathrm{YX}}& {\Gamma }_{\mathrm{YY}}& {\Gamma }_{\mathrm{YZ}}\\ {\Gamma }_{\mathrm{ZX}}& {\Gamma }_{\mathrm{ZY}}& {\Gamma }_{\mathrm{ZZ}}\end{array}\right]$ wherein, Γ.sub.XX, Γ.sub.YY and Γ.sub.ZZ are respectively the gravity gradient on the x axis, the y axis and the z axis, Γ.sub.XY and Γ.sub.YX are the gravity gradient on the x axis and the y axis, Γ.sub.XZ and Γ.sub.ZX the gravity gradient on the x axis and the z axis, and Γ.sub.YZ and Γ.sub.ZY are the gravity gradient on they axis and the z axis, satisfying Γ.sub.XX+Γ.sub.YY+Γ.sub.ZZ=0, Γ.sub.XY=Γ.sub.YX, Γ.sub.xz=Γ.sub.zx, Γ.sub.YZ=Γ.sub.ZY.

10. The gravity gradient measurement apparatus of claim 9, wherein, specifically, the measurement module adopts the following formula
a.sub.1X−a.sub.2X=h sin ωtΓ.sub.YZ−h cos ωtΓ.sub.XZ+2R sin 2ωtΓ.sub.XY−R cos 2ωt(Γ.sub.XX−Γ.sub.YY)−R(Γ.sub.XX+Γ.sub.YY)−2Rω.sub.z.sup.2
a.sub.1Y−a.sub.2Y=h cos ωtΓ.sub.YZ−h sin ωtΓ.sub.XZ+2R cos 2ωtΓ.sub.XY−R sin 2ωt(Γ.sub.XX−Γ.sub.YY)+2R{dot over (ω)}.sub.z
a.sub.1Z−a.sub.2Z=2R sin ωtΓ.sub.YZ−2R cos ωtΓ.sub.XZ−hΓ.sub.ZZ to determine components Γ.sub.YZ, Γ.sub.XZ, Γ.sub.XY and Γ.sub.XX−Γ.sub.YY of the gravity gradient; and determines a component Γ.sub.ZZ of the gravity gradient using the formula Γ.sub.ZZ=(a.sub.2Z−a.sub.1Z)/h; wherein, a.sub.1X, a.sub.1Y, and a.sub.1Z respectively represent measurement outputs of the first three-axis accelerometer a.sub.1 in the directions of the x axis, the y axis and the z axis; a.sub.2X, a.sub.2Y, and a.sub.2Z respectively represent measurement outputs of the second three-axis accelerometer a.sub.2 in the directions of the x axis, the y axis and the z axis; and ω is an angular velocity vector, and {dot over (ω)}.sub.z is angular acceleration around the z axis.

11. A gravity gradient measurement method, comprising: measuring acceleration using a first three-axis accelerometer (1) and a second three-axis accelerometer (2); and determining a gravity gradient on the coordinate axes according to the acceleration; wherein, the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are arranged in a first chamber defined by a turntable (3) and a vacuum layer (7); the turntable (3) rotates horizontally around an earth-vertical axis; the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are arranged symmetrically on an x axis with respect to an origin of coordinates; both the first three-axis accelerometer (1) and the second three-axis accelerometer (2) have a distance of R from the origin of coordinates; the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are arranged symmetrically on a z axis with respect to an origin of coordinates; and the first three-axis accelerometer and the second three-axis accelerometer are spaced at a distance of h on the z axis; a first magnetic shield (4) is arranged in the first chamber for shielding the interference of the external geomagnetic field, the first magnetic shield (4) and the turntable (3) define a second chamber, in which the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are located; a second magnetic shield (5) is arranged in the second chamber for shielding the interference of an external alternating current magnetic field, the second magnetic shield (5) and the turntable (3) define a third chamber, in which the first three-axis accelerometer (1) and the second three-axis accelerometer (2) are located.

12. The method of claim 11, wherein, the gravity gradient Γ.sup.g is: ${\Gamma }^g=\left[\begin{array}{ccc}{\Gamma }_{\mathrm{XX}}& {\Gamma }_{\mathrm{XY}}& {\Gamma }_{\mathrm{XZ}}\\ {\Gamma }_{\mathrm{YX}}& {\Gamma }_{\mathrm{YY}}& {\Gamma }_{\mathrm{YZ}}\\ {\Gamma }_{\mathrm{ZX}}& {\Gamma }_{\mathrm{ZY}}& {\Gamma }_{\mathrm{ZZ}}\end{array}\right]$ wherein, Γ.sub.XX, Γ.sub.YY and Γ.sub.ZZ are respectively the gravity gradient on the x axis, the y axis and the z axis, Γ.sub.XY and Γ.sub.YX are the gravity gradient on the x axis and the y axis, Γ.sub.XZ and Γ.sub.ZX the gravity gradient on the x axis and the z axis, and Γ.sub.YZ and Γ.sub.ZY are the gravity gradient on the y axis and the z axis, satisfying Γ.sub.XX+Γ.sub.YY+Γ.sub.ZZ=0, Γ.sub.XY=Γ.sub.YX, Γ.sub.xz=Γ.sub.zx, Γ.sub.YZ=Γ.sub.ZY.

13. The method of claim 12, wherein, determining the gravity gradient on the coordinate axes according to acceleration values measured by the first three-axis accelerometer (1) and the second three-axis accelerometer (2) comprises: using the formula
a.sub.1X−a.sub.2X=h sin ωtΓ.sub.YZ−h cos ωtΓ.sub.XZ+2R sin 2ωtΓ.sub.XY−R cos 2ωt(Γ.sub.XX−Γ.sub.YY)−R(Γ.sub.XX+Γ.sub.YY)−2Rω.sub.z.sup.2
a.sub.1Y−a.sub.2Y=h cos ωtΓ.sub.YZ−h sin ωtΓ.sub.XZ+2R cos 2ωtΓ.sub.XY−R sin 2ωt(Γ.sub.XX−Γ.sub.YY)+2R{dot over (ω)}.sub.z
a.sub.1Z−a.sub.2Z=2R sin ωtΓ.sub.YZ−2R cos ωtΓ.sub.XZ−hΓ.sub.ZZ to determine components Γ.sub.YZ, Γ.sub.XZ, Γ.sub.XY and Γ.sub.XX−Γ.sub.YY of the gravity gradient; and determines a component Γ.sub.ZZ of the gravity gradient using the formula Γ.sub.ZZ=(a.sub.2Z−a.sub.1Z)/h; wherein, a.sub.1X, a.sub.1Y, and a.sub.1Z respectively represent measurement outputs of the first three-axis accelerometer a.sub.1 in the directions of the x axis, the y axis and the z axis; a.sub.2X, a.sub.2Y, and a.sub.2Z respectively represent measurement outputs of the second three-axis accelerometer a.sub.2 in the directions of the x axis, the y axis and the z axis; and ω is an angular velocity vector, and {dot over (ω)}.sub.z is angular acceleration around the z axis.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The drawings illustrated below are intended to facilitate better understanding of the present invention, and constitute a part of the present application. The exemplary embodiments as well as corresponding descriptions in the present invention are used to interpret the present invention, and shall not constitute undue limitation on the present invention. In the drawings:

(2) FIG. 1 is a view of a gravity gradient measurement apparatus of one embodiment of the present invention.

(3) FIG. 2 is a view of a gravity gradient measurement apparatus of another embodiment of the present invention.

(4) FIG. 3 is a view of a gravity gradient measurement apparatus of a further embodiment of the present invention.

(5) FIG. 4 is a structural view of a three-axis accelerometer electrode of the present invention.

(6) FIG. 5 is a diagram of a gravity gradient measuring method of one embodiment of the present invention.

(7) FIGS. 6a-6c are views of position arrangement of a three-axis accelerometer of the present invention.

(8) FIG. 7 is a diagram of coordinate rotation of a three-axis accelerometer of the present invention.

DETAILED DESCRIPTION

(9) Further detailed description of the present invention is made in conjunction with the drawings and embodiments.

(10) FIG. 1 is a view of a gravity gradient measurement apparatus of one embodiment of the present invention. As shown in FIG. 1, the gravity gradient measurement apparatus includes a turntable 3, a first three-axis accelerometer 1, a second three-axis accelerometer 2, a vacuum layer 7 and a measurement module, and the turntable 3 rotates horizontally around an earth-vertical axis, wherein:

(11) a vacuum layer 7 is arranged on the turntable 3 and defines a first chamber with the turntable 3, the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are located in the first chamber; the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are arranged symmetrically on an x axis with respect to an origin of coordinates; both the first three-axis accelerometer 1 and the second three-axis accelerometer 2 have a distance of R from the origin of coordinates; the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are arranged symmetrically on an z axis with respect to the origin of coordinates, and the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are spaced at a distance of h on the z axis.

(12) The measurement module (not shown in the drawings) determines a gravity gradient on the coordinate axes according to acceleration values measured by the first three-axis accelerometer 1 and the second three-axis accelerometer 2.

(13) FIG. 2 is a view of a gravity gradient measurement apparatus of another embodiment of the present invention. Compared with the embodiment illustrated in FIG. 1, in the embodiment illustrated in FIG. 2, the gravity gradient measurement apparatus further includes an first magnetic shield 4 arranged in the first chamber for shielding the interference of the external geomagnetic field. The first magnetic shield 4 and the turntable 3 define a second chamber, in which the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are located.

(14) Optionally, the first magnetic shield 4 is a metal shield layer. For example, the material of the first magnetic shield 4 is permalloy.

(15) Optionally, the gravity gradient measurement apparatus may also include a second magnetic shield 5 arranged in the second chamber for shielding the interference of an external alternating current magnetic field. The second magnetic shield 5 and the turntable 3 define a third chamber, in which the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are located.

(16) Optionally, the second magnetic shield 5 is a superconducting shield layer. For example, the material of the second magnetic shield 5 is a YBCO (yttrium barium copper oxide) metallic compound.

(17) FIG. 3 is a view of a gravity gradient measurement apparatus of a further embodiment of the present invention. Compared with the embodiment illustrated in FIG. 2, in the embodiment illustrated in FIG. 3, the gravity gradient measurement apparatus further includes a liquid nitrogen layer 6 arranged between the first magnetic shield 4 and the second magnetic shield 5.

(18) The liquid nitrogen layer 6 is a cylindrical double-layer metal thin-wall layer, in which liquid nitrogen is filled. The liquid nitrogen layer 6 achieves an internal low temperature environment of about 77K, enabling the YBCO metallic compound of the second magnetic shield 5 to reach a superconducting state, reducing the thermal noise during Brownian motion of the first three-axis accelerometer 1 and the second three-axis accelerometer 2, and improving the measurement accuracy of the first three-axis accelerometer 1 and the second three-axis accelerometer 2.

(19) As shown in FIG. 4, the above first three-axis accelerometer 1 and the second three-axis accelerometer 2 have the same structure and functions, and respectively include a suspended mass 8, an x-axis acceleration detection electrode 9, a y-axis acceleration detection electrode 10 and a z-axis acceleration detection electrode 11. The x-axis acceleration detection electrode 9, the y-axis acceleration detection electrode 10 and the z-axis acceleration detection electrode 11 define a cavity, the suspended mass 8 is suspended at a central position of the cavity. The suspended mass 8 is a hollow metal mass with an enclosed outer surface. The x-axis acceleration detection electrode 9, the y-axis acceleration detection electrode 10 and the z-axis acceleration detection electrode 11 realize measurement of the acceleration on the x axis, y axis and z axis through high frequency capacitance bridge circuit, i.e., measure the displacement of the suspended mass 8 in the directions of the x axis, the y axis and the z axis, and obtain the acceleration in the directions of the x axis, the y axis and the z axis through second-order derivation of the displacement against time. The x-axis acceleration detection electrode 9, the y-axis acceleration detection electrode 10 and the z-axis acceleration detection electrode 11 realize the functions of electrostatic force support and feedback control in the directions of orthogonal x axis, y axis and z axis via low frequency voltage loaded on the electrode.

(20) To be noted, the gravity gradient Γ.sup.g is:

(21) ${\Gamma }^g=\left[\begin{array}{ccc}{\Gamma }_{\mathrm{XX}}& {\Gamma }_{\mathrm{XY}}& {\Gamma }_{\mathrm{XZ}}\\ {\Gamma }_{\mathrm{YX}}& {\Gamma }_{\mathrm{YY}}& {\Gamma }_{\mathrm{YZ}}\\ {\Gamma }_{\mathrm{ZX}}& {\Gamma }_{\mathrm{ZY}}& {\Gamma }_{\mathrm{ZZ}}\end{array}\right]$

(22) wherein, Γ.sub.XX, Γ.sub.YY and Γ.sub.ZZ are respectively the gravity gradient on the x axis, the y axis and the z axis, Γ.sub.XY and Γ.sub.YX are the gravity gradient on the x axis and the y axis, Γ.sub.XZ and Γ.sub.ZX the gravity gradient on the x axis and the z axis, and Γ.sub.YZ and Γ.sub.ZY, are the gravity gradient on the y axis and the z axis, satisfying Γ.sub.XX+Γ.sub.YY+Γ.sub.ZZ=0, Γ.sub.XY=Γ.sub.YX, Γ.sub.xz=Γ.sub.zx, Γ.sub.YZ=Γ.sub.ZY.

(23) Optionally, the measurement module may adopt the following formula
a.sub.1X−a.sub.2X=h sin ωtΓ.sub.YZ−h cos ωtΓ.sub.XZ+2R sin 2ωtΓ.sub.XY−R cos 2ωt(Γ.sub.XX−Γ.sub.YY)−R(Γ.sub.XX+Γ.sub.YY)−2Rω.sub.z.sup.2
a.sub.1Y−a.sub.2Y=h cos ωtΓ.sub.YZ−h sin ωtΓ.sub.XZ+2R cos 2ωtΓ.sub.XY−R sin 2ωt(Γ.sub.XX−Γ.sub.YY)+2R{dot over (ω)}.sub.z
a.sub.1Z−a.sub.2Z=2R sin ωtΓ.sub.YZ−2R cos ωtΓ.sub.XZ−hΓ.sub.ZZ

(24) to determine components Γ.sub.YZ, Γ.sub.XZ, Γ.sub.XY and Γ.sub.XX−Γ.sub.YY of the gravity gradient;

(25) and determines a component Γ.sub.ZZ of the gravity gradient using the formula Γ.sub.ZZ=(a.sub.2Z−a.sub.1Z)/h;

(26) wherein, a.sub.1X, a.sub.1Y, and a.sub.1Z respectively represent measurement outputs of the first three-axis accelerometer a.sub.1 in the directions of the x axis, the y axis and the z axis;

(27) a.sub.2X, a.sub.2Y, and a.sub.2Z respectively represent measurement outputs of the second three-axis accelerometer a.sub.2 in the directions of the x axis, the y axis and the z axis; and

(28) co is an angular velocity vector, and {dot over (ω)}.sub.z is angular acceleration around the z axis.

(29) FIG. 5 is a diagram for an embodiment of a gravity gradient measuring method of the present invention. The gravity gradient measuring method can be performed using the gravity gradient measurement apparatus involved in any embodiment in FIGS. 1-4.

(30) Step 501, measuring acceleration using a first three-axis accelerometer and a second three-axis accelerometer.

(31) Step 502, determining a gravity gradient on the coordinate axes according to the acceleration.

(32) The measuring method adopted by the present invention will be described specifically as below.

(33) As shown in FIGS. 6a-6c, the position coordinates of the first three-axis accelerometer 1 and the second three-axis accelerometer 2 on the x axis, y axis and z axis are (R, 0, h/2) and (−R, 0, −h/2) respectively. Determining for the accelerometers a coordinate system X.sub.a-O-Y.sub.a, which is referred to as a coordinate system a and is located on the same plane as the coordinate system X-O-Y as shown in FIG. 6a-6c. However, the coordinate system a is fixedly connected with the three-axis accelerometers, which is to say, the accelerometer coordinate system X.sub.a-O-Y.sub.a rotates along with the three-axis accelerometers, and the X.sub.g-O-Y.sub.g is an inertial coordinate system.

(34) FIG. 6a is a stereogram, FIG. 6b is a front view, and FIG. 6c is a top view.

(35) Firstly, the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are respectively used to measure the difference among linear accelerations in the directions of the x axis, y axis and z axis during rotation, and the linear accelerations in the directions of the x axis, y axis and z axis measured by the two three-axis accelerometers undergo subtraction, so as to obtain differential acceleration, and the relation between the differential acceleration and a gravity gradient parameter is shown with the following formula:

(36) $\begin{array}{cc}{f}_1^a-f_2^a=\left[{\Gamma }^a-\left(\stackrel{.}{\omega }\times \right)-\left({\omega }_{\mathrm{ia}}^a\times \right)\left({\omega }_{\mathrm{ia}}^a\times \right)\right]{\rho }^a=\left[{\Gamma }^a-{\stackrel{.}{\Omega }}_{\mathrm{ia}}^a-{\Omega }_{\mathrm{ia}}^a{\Omega }_{\mathrm{ia}}^a\right]{\rho }^a=L^a{\rho }^a& \left(1\right)\end{array}$

(37) wherein, f.sub.1.sup.a,f.sub.2.sup.a are respectively a measurement value matrix of the first three-axis accelerometer 1 and the second three-axis accelerometer 2 in the coordinate system a; Γ.sup.a is a gravity gradient matrix of the first three-axis accelerometer 1 and the second three-axis accelerometer 2 in the coordinate system a, the gravity gradient matrix in the inertial system may be obtained through coordinate conversion; ρ.sup.a is a difference matrix of the position vector of the first three-axis accelerometer 1 and the second three-axis accelerometer 2. The value is constant when the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are located at relatively fixed positions with respect to each other; L.sup.a is a gradient measurement matrix obtained after direct differentiation of the displacement by the measurement of the first three-axis accelerometer 1 and the second three-axis accelerometer 2. The invention further includes an angular velocity Ω.sub.ia.sup.aΩ.sub.ia.sup.a component and an angular acceleration component {dot over (Ω)}.sub.ia.sup.a; {dot over (ω)} is a rotation acceleration of the inertial coordinate system, and ω.sub.ia.sup.a is an angular velocity of the coordinate system a.

(38) wherein, the angular velocity matrix is:

(39) $\begin{array}{cc}{\Omega }_{\mathrm{ia}}^a{\Omega }_{\mathrm{ia}}^a=\left[\begin{array}{ccc}-{\omega }_y^2-{\omega }_z^2& {\omega }_x{\omega }_y& {\omega }_x{\omega }_z\\ {\omega }_x{\omega }_y& -{\omega }_x^2-{\omega }_z^2& {\omega }_y{\omega }_z\\ {\omega }_x{\omega }_z& {\omega }_y{\omega }_z& -{\omega }_x^2-{\omega }_y^2\end{array}\right]& \left(2\right)\end{array}$

(40) The angular acceleration matrix is:

(41) $\begin{array}{cc}{\stackrel{.}{\Omega }}_{\mathrm{ia}}^a=\left[\begin{array}{ccc}0& -{\stackrel{.}{\omega }}_z& {\stackrel{.}{\omega }}_y\\ {\stackrel{.}{\omega }}_z& 0& -{\stackrel{.}{\omega }}_x\\ -{\stackrel{.}{\omega }}_y& {\stackrel{.}{\omega }}_x& 0\end{array}\right]& \left(3\right)\end{array}$

(42) The measurement matrix of the first three-axis accelerometer 1 and the second three-axis accelerometer 2 in the coordinate system a may be described as:
f.sub.1.sup.a=(a.sub.1X,a.sub.1Y,a.sub.1Z).sub.T,f.sub.2.sup.a=(a.sub.2X,a.sub.2Y,a.sub.2Z).sup.T  (4)

(43) The displacement differentiation matrix is:

(44) $\begin{array}{cc}{\rho }^a=r_2^a-r_1^a=\left[\begin{array}{c}-R\\ 0\\ -h/2\end{array}\right]-\left[\begin{array}{c}R\\ 0\\ h/2\end{array}\right]=\left[\begin{array}{c}-2R\\ 0\\ -h\end{array}\right]& \left(5\right)\end{array}$

(45) (4) and (5) are substituted in formula (1) to obtain:

(46) $\begin{array}{cc}\left[\begin{array}{c}{a}_{1X}-a_{2X}\\ a_{1Y}-a_{2Y}\\ a_{1Z}-a_{2Z}\end{array}\right]=L^a\left[\begin{array}{c}-2R\\ 0\\ -h\end{array}\right]& \left(6\right)\end{array}$

(47) L.sup.a is a gradient measurement matrix, which also includes components of the angular velocity and the angular acceleration, substituting the angular velocity matrix (2) and the angular acceleration matrix (3) into formula (6) to obtain:
a.sub.1X−a.sub.2X=−2RL.sub.11.sup.a−hL.sub.13.sup.a=−2R(Γ.sub.XX.sup.a+ω.sub.y.sup.2+ω.sub.z.sup.2)−h(Γ.sub.XZ.sup.a−ω.sub.xω.sub.z−{dot over (ω)}.sub.y)
a.sub.1Y−a.sub.2Y=−2RL.sub.21.sup.a−hL.sub.23.sup.a=−2R(Γ.sub.YX.sup.a−ω.sub.xω.sub.y−{dot over (ω)}.sub.z)−h(Γ.sub.YZ.sup.a−ω.sub.yω.sub.z+{dot over (ω)}.sub.x)
a.sub.1Z−a.sub.2Z=−2RL.sub.31.sup.a−hL.sub.33.sup.a=−2R(Γ.sub.ZX.sup.a−ω.sub.xω.sub.z+{dot over (ω)}.sub.y)−h(Γ.sub.ZZ.sup.a+ω.sub.x.sup.2+ω.sub.y.sup.2)  (7)

(48) Hereby direct relation between the measurement of the first three-axis accelerometer 1 and the second three-axis accelerometer 2 and the gravity gradient in the coordinate system a is obtained, and the coordinate conversion from the gravity gradient in the coordinate system a into that in the inertial system is conducted as bellow.

(49) The coordinate conversion from the gravity gradient in the coordinate system a into the gravity gradient in the inertial coordinate system is expressed in the following formula:
Γ.sup.a=C.sub.g.sup.aΓ.sup.gC.sub.a.sup.g  (8)

(50) As shown in FIG. 7, when the coordinate system a rotates around the inertial system in an angular velocity of ω, the transposed matrix of the conversion formula is:

(51) 0$\begin{array}{cc}{C}_g^a=\left[\begin{array}{ccc}\cos \omega t& -\sin \omega t& 0\\ \sin \omega t& \cos \omega t& 0\\ 0& 0& 1\end{array}\right],C_a^g=\left[\begin{array}{ccc}\cos \omega t& \sin \omega t& 0\\ -\sin \omega t& \cos \omega t& 0\\ 0& 0& 1\end{array}\right]& \left(9\right)\end{array}$

(52) Substituting (9) in (8):

(53) $\begin{array}{cc}\begin{array}{c}{\Gamma }^a=\left[\begin{array}{ccc}{\Gamma }_{\mathrm{XX}}^a& {\Gamma }_{\mathrm{XY}}^a& {\Gamma }_{\mathrm{XZ}}^a\\ {\Gamma }_{\mathrm{YX}}^a& {\Gamma }_{\mathrm{YY}}^a& {\Gamma }_{\mathrm{YZ}}^a\\ {\Gamma }_{\mathrm{ZX}}^a& {\Gamma }_{\mathrm{ZY}}^a& {\Gamma }_{\mathrm{ZZ}}^a\end{array}\right]\\ =\left[\begin{array}{ccc}\cos \omega t& -\sin \omega t& 0\\ \sin \omega t& \cos \omega t& 0\\ 0& 0& 1\end{array}\right].Math.\left[\begin{array}{ccc}{\Gamma }_{\mathrm{XX}}& {\Gamma }_{\mathrm{XY}}& {\Gamma }_{\mathrm{XZ}}\\ {\Gamma }_{\mathrm{YX}}& {\Gamma }_{\mathrm{YY}}& {\Gamma }_{\mathrm{YZ}}\\ {\Gamma }_{\mathrm{ZX}}& {\Gamma }_{\mathrm{ZY}}& {\Gamma }_{\mathrm{ZZ}}\end{array}\right].Math.\\ \left[\begin{array}{ccc}\cos \omega t& \sin \omega t& 0\\ -\sin \omega t& \cos \omega t& 0\\ 0& 0& 1\end{array}\right]\end{array}\mathrm{Thus}\text{:}{\Gamma }_{\mathrm{XX}}^a=\left(1/\right)\left({\Gamma }_{\mathrm{XX}}-{\Gamma }_{\mathrm{YY}}\right)\cos 2\omega t-{\Gamma }_{\mathrm{XY}}\sin 2\omega t+\left(1/2\right)\left({\Gamma }_{\mathrm{XX}}+{\Gamma }_{\mathrm{YY}}\right){\Gamma }_{\mathrm{XY}}^a=\left(1/2\right)\left({\Gamma }_{\mathrm{XX}}-{\Gamma }_{\mathrm{YY}}\right)\sin 2\omega t-{\Gamma }_{\mathrm{XY}}\cos 2\omega t{\Gamma }_{\mathrm{XZ}}^a=\cos \omega t{\Gamma }_{\mathrm{XZ}}-\sin \omega t{\Gamma }_{\mathrm{YZ}}{\Gamma }_{\mathrm{YZ}}^a=\sin \omega t{\Gamma }_{\mathrm{XZ}}+\cos \omega t{\Gamma }_{\mathrm{YZ}}{\Gamma }_{\mathrm{ZZ}}^a={\Gamma }_{\mathrm{ZZ}}& \left(10\right)\end{array}$

(54) Because the first three-axis accelerometer 1 and the second three-axis accelerometer 2 only rotate around the z axis, ω.sub.x=ω.sub.y=0, ω.sub.z is constantly equal to ω, (10) is substituted in (7), therefore:
a.sub.1X−a.sub.2X=h sin ωtΓ.sub.YZ−h cos ωtΓ.sub.XZ+2R sin 2ωtΓ.sub.XY−R cos 2ωt(Γ.sub.XX−Γ.sub.YY)−R(Γ.sub.XX+Γ.sub.YY)−2Rω.sub.z.sup.2
a.sub.1Y−a.sub.2Y=h cos ωtΓ.sub.YZ−h sin ωtΓ.sub.XZ+2R cos 2ωtΓ.sub.XY−R sin 2ωt(Γ.sub.XX−Γ.sub.YY)+2R{dot over (ω)}.sub.z
a.sub.1Z−a.sub.2Z=2R sin ωtΓ.sub.YZ−2R cos ωtΓ.sub.XZ−hΓ.sub.ZZ  (11)

(55) After analysing the results of Formula (11), it can be seen that four gravity gradient components Γ.sub.YZ, Γ.sub.XZ, Γ.sub.XY, (Γ.sub.XX−Γ.sub.YY) are respectively added to 1 frequency and double-frequency orthogonal signals of rotational angular velocity ω of combined outputs of the first three-axis accelerometer 1 and the second three-axis accelerometer 2. First, the four gravity gradient components Γ.sub.YZ, Γ.sub.XZ, Γ.sub.XY, (Γ.sub.XX−Γ.sub.YY) undergo frequency division to respectively obtain the 1 frequency and double-frequency orthogonal signals which are respectively amplified and demodulated, following by separation via two detection signals with a phase difference of 90°, thus obtaining four 4 gravity gradient components. The Γ.sub.ZZ component is a constant of the combined acceleration output a.sub.1Z-a.sub.2Z in the z-axis direction, i.e., Γ.sub.ZZ=(a.sub.2Z−a.sub.1Z)/h, the difference between the outputs of a.sub.1Z and a.sub.2Z is directly filtered to obtain the component. Hereby all five independent components of the gravity gradient are obtained, thus obtaining a full-tensor gravity gradient.

(56) The symbols in the above formula and representation have the following meanings:

(57) X,Y,Z: three axes in the space;

(58) O: the origin of the axes in the space;

(59) a.sub.1, a.sub.2: respectively represent the first three-axis accelerometer 1 and the second three-axis accelerometer 2;

(60) a.sub.1X, a.sub.1Y, a.sub.1Z: respectively represent measurement outputs of the first three-axis accelerometer a.sub.1 in the directions of the X,Y,Z axis;

(61) a.sub.2X, a.sub.2Y, a.sub.2Z: respectively represent measurement outputs of the second three-axis accelerometer a.sub.2 in the directions of the X, Y, Z axis;

(62) R: the distance of the first three-axis accelerometer a.sub.1 and the second three-axis accelerometer a.sub.2 on the z axis with respect to the origin O;

(63) h: the distance of the first three-axis accelerometer a.sub.1 and the second three-axis accelerometer a.sub.2 on the z axis;

(64) ω: the angular velocity vector, generally has three components ω.sub.x, ω.sub.y, ω.sub.z, and is a 3×1 vector, with a scalar value equal to ω.sub.z, because the apparatus of the present invention only rotates around the Z axis;

(65) {dot over (ω)}: the angular acceleration vector, is a first-order derivative of the angular velocity ω, generally has three components {dot over (ω)}.sub.x, {dot over (ω)}.sub.y, {dot over (ω)}.sub.z, and is a 3×1 vector, wherein {dot over (ω)}.sub.z=0 when ideal uniform angular velocity rotation is conducted;

(66) ω.sub.x, ω.sub.y, ω.sub.z: are respectively the angular velocity around the x axis, y axis, and z axis. The direction of the angular velocity vector can be determined by the right-hand rule, i.e., the axis around which the rotation is conducted is the direction of the angular velocity of the axis;

(67) {dot over (ω)}.sub.x, {dot over (ω)}.sub.y, {dot over (ω)}.sub.z: are respectively the angular acceleration around the x axis, y axis and z axis;

(68) X.sub.a-O-Y.sub.a: is the coordinate system a on the X-O-Y plane fixedly connected with the three-axis accelerometers, and rotating along with the three-axis accelerometers, i.e., the three-axis accelerometers have a position vector remaining constant in the X.sub.a-O-Y.sub.a coordinate system;

(69) X.sub.g-O-Y.sub.g: is the coordinate system g on the X-O-Y plane, and is a static inertial system not rotating along with the three-axis accelerometers;

(70) f.sub.1.sup.a,f.sub.2.sup.a: are respectively an measurement matrix of the two three-axis accelerometer in the coordinate system a;

(71) Γ.sup.a: is the gravity gradient matrix of the three-axis accelerometers in the coordinate system a;

(72) Γ.sup.g: is the coordinate system g, i.e., the gravity gradient matrix in the inertial coordinate system;

(73) r.sub.1.sup.a: the position vector of the first three-axis accelerometer a.sub.1 in the coordinate system a, and is (R, 0, ½h);

(74) r.sub.2.sup.a: the position vector of the second three-axis accelerometer a.sub.2 in the coordinate system a, and is (−R, 0, −½h);

(75) ρ.sup.a: is the difference matrix of the position vector of the first three-axis accelerometer a.sub.1 and the second three-axis accelerometer a.sub.2 in the coordinate system a;

(76) Ω: the angular velocity matrix;

(77) {dot over (Ω)}: the angular acceleration matrix;

(78) L.sup.a: the gradient measurement matrix, which is obtained through direct division of combinations of the outputs from the accelerometer by the differential displacement, including the angular velocity and the angular acceleration and other components, and is endowed with a superscript a on the top right, as it is in the a coordinate system;

(79) The gradient measurement matrix L.sup.a includes:

(80) $\hspace{1em}\left[\begin{array}{ccc}{L}_{11}^a& L_{12}^a& L_{13}^a\\ L_{21}^a& L_{22}^a& L_{23}^a\\ L_{31}^a& L_{32}^a& L_{33}^a\end{array}\right]$

(81) The gradient matrix Γ.sup.a in the coordinate system a includes:

(82) $\hspace{1em}\left[\begin{array}{ccc}{\Gamma }_{\mathrm{XX}}^a& {\Gamma }_{\mathrm{XY}}^a& {\Gamma }_{\mathrm{XZ}}^a\\ {\Gamma }_{\mathrm{YX}}^a& {\Gamma }_{\mathrm{YY}}^a& {\Gamma }_{\mathrm{YZ}}^a\\ {\Gamma }_{\mathrm{ZX}}^a& {\Gamma }_{\mathrm{ZY}}^a& {\Gamma }_{\mathrm{ZZ}}^a\end{array}\right]$

(83) The gradient matrix Γ.sup.g in the inertial coordinate system g includes:

(84) $\left[\begin{array}{ccc}{\Gamma }_{\mathrm{XX}}& {\Gamma }_{\mathrm{XY}}& {\Gamma }_{\mathrm{XZ}}\\ {\Gamma }_{\mathrm{YX}}& {\Gamma }_{\mathrm{YY}}& {\Gamma }_{\mathrm{YZ}}\\ {\Gamma }_{\mathrm{ZX}}& {\Gamma }_{\mathrm{ZY}}& {\Gamma }_{\mathrm{ZZ}}\end{array}\right],$
wherein Γ.sup.g is also a target gradient matrix in the measuring method of the present invention.

(85) The present invention enables convenient measurement of a full-tensor gravity gradient through two three-axis accelerometers.

(86) The description of the present invention is made in an attempt to facilitate exemplary description thereof, and is not exhaustive or intended to confine the present invention to the disclosed implementations. Multiple modifications and changes are obvious to those skilled in the art. And selecting and describing the embodiments are aimed at better illustrate the principle and practical applications of the present invention, and enabling one with ordinary skill in the art to understand the present invention and to devise various embodiments with various modifications suitable for specific uses.