WO2018161474A1 - 重力梯度测量方法及装置 - Google Patents

重力梯度测量方法及装置 Download PDF

Info

Publication number
WO2018161474A1
WO2018161474A1 PCT/CN2017/091107 CN2017091107W WO2018161474A1 WO 2018161474 A1 WO2018161474 A1 WO 2018161474A1 CN 2017091107 W CN2017091107 W CN 2017091107W WO 2018161474 A1 WO2018161474 A1 WO 2018161474A1
Authority
WO
WIPO (PCT)
Prior art keywords
axis
gravity gradient
accelerometer
axis accelerometer
acceleration
Prior art date
Application number
PCT/CN2017/091107
Other languages
English (en)
French (fr)
Inventor
胡新宁
王秋良
陆锦焱
王晖
崔春艳
Original Assignee
中国科学院电工研究所
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 中国科学院电工研究所 filed Critical 中国科学院电工研究所
Priority to US16/491,936 priority Critical patent/US11402538B2/en
Publication of WO2018161474A1 publication Critical patent/WO2018161474A1/zh

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
    • G01V7/02—Details
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
    • H—ELECTRICITY
    • H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05K—PRINTED CIRCUITS; CASINGS OR CONSTRUCTIONAL DETAILS OF ELECTRIC APPARATUS; MANUFACTURE OF ASSEMBLAGES OF ELECTRICAL COMPONENTS
    • H05K9/00—Screening of apparatus or components against electric or magnetic fields
    • H05K9/0073—Shielding materials
    • H05K9/0075—Magnetic shielding materials
    • H05K9/0077—Magnetic shielding materials comprising superconductors

Definitions

  • the invention relates to the field of gravity measurement, and in particular to a gravity gradient measuring method and device.
  • Gravity gradient measurement is a measure of the change in Earth's gravitational acceleration with space. Since the gravity gradient is the spatial differentiation of the Earth's gravitational field, reflecting the rate of change of gravity along different directions of the space, the gravity gradient measurement can reflect the details of the field source.
  • the main advantage of gravity gradient measurement is that gravity gradient values or gravity higher order derivatives have higher resolution than gravity itself. Therefore, by using the gravity gradiometer measuring the second derivative of the gravitational potential to measure the gravity gradient tensor component in real time, more accurate gravity values and vertical deviations can be obtained, which plays an important role in the development of science and technology such as space science, earth science and geological science. .
  • the gravity gradient is the first-order spatial derivative of the gravitational acceleration vector, that is, the second-order spatial derivative of the gravitational potential, which can be expressed by the following tensor:
  • the gravity gradient measurement method is mainly based on the principle of differential acceleration measurement, and the respective components of the gravitational gradient tensor can be obtained by different combinations of the observations between the accelerometers.
  • Full tensor gravity gradient measurement requires three orthogonal rotation axes to separate the gradient signal, especially the non-diagonal component of the gravity gradient; or an additional angular accelerometer to measure the non-diagonal component and eliminate the dynamic error, generally requires 6
  • the combined structure of the above accelerometers greatly increases the complexity and cost of the instrumentation system.
  • a gravity gradient measuring apparatus comprising a rotating table, a first three-axis accelerometer, a second three-axis accelerometer, a vacuum layer and a measuring module, wherein the rotating table is horizontally rotated about a vertical axis, wherein:
  • the vacuum layer is disposed on the rotating table, the vacuum layer and the rotating table form a first chamber, and the first three-axis accelerometer and the second three-axis accelerometer are located inside the first chamber, the first three-axis accelerometer and the second three-axis accelerometer
  • the accelerometer is symmetrically arranged on the x-axis with respect to the origin of the coordinate axis, and the distance between the first three-axis accelerometer and the second three-axis accelerometer and the origin of the coordinate axis is R; the first three-axis accelerometer and the second three-axis accelerometer Arranging symmetrically with respect to the origin of the coordinate axis on the z-axis, the distance between the first three-axis accelerometer and the second three-axis accelerometer on the z-axis is h;
  • the measurement module determines the gravity gradient value under the coordinate axis using the acceleration values measured by the first three-axis accelerometer and the second three-axis accelerometer.
  • the gravity gradient measuring device further includes:
  • first magnetic shielding layer disposed inside the first chamber to shield external magnetic field interference, the first magnetic shielding layer and the rotating table forming a second chamber, the first three-axis accelerometer and the second three-axis accelerometer Located inside the second chamber.
  • the first magnetic shielding layer is a metal shielding layer.
  • the material of the first magnetic shielding layer is permalloy.
  • the gravity gradient measuring device further includes:
  • a second magnetic shielding layer disposed inside the second chamber for shielding external AC magnetic field interference, the second magnetic shielding layer and the rotating table forming a third chamber, the first three-axis accelerometer and the second three-axis accelerometer are located Inside the third chamber.
  • a liquid nitrogen layer disposed between the first magnetic shield layer and the second magnetic shield layer.
  • the second magnetic shielding layer is a superconducting shielding layer.
  • the material of the second magnetic shielding layer is a YBCO metal compound.
  • the first three-axis accelerometer and the second three-axis accelerometer both include a floating mass body, an x-axis acceleration detecting electrode, a y-axis acceleration detecting electrode, and a z-axis acceleration detecting electrode, an x-axis acceleration detecting electrode, and a y-axis acceleration.
  • the detecting electrode and the z-axis acceleration detecting electrode form a cavity, and the suspended mass body is suspended at a center position of the cavity;
  • the x-axis acceleration detecting electrode, the y-axis acceleration detecting electrode, and the z-axis acceleration detecting electrode measure the displacement of the suspended mass body in the x-axis, the y-axis, and the z-axis direction, and then obtain the x-axis, the y-axis, and the z-axis according to the measured displacement. Acceleration in the direction.
  • the suspended mass is a hollow outer surface closed metal mass.
  • the gravity gradient value ⁇ g is:
  • ⁇ XX , ⁇ YY and ⁇ ZZ are the gravity gradient values on the x-axis, y-axis and z-axis, respectively
  • ⁇ XY and ⁇ YX are the gravity gradient values on the x-axis and y-axis
  • ⁇ XY ⁇ YX
  • ⁇ XY ⁇ YX
  • ⁇ YZ ⁇ ZY .
  • the measurement module specifically uses a formula
  • a 1Z -a 2Z 2R sin ⁇ t ⁇ YZ -2R cos ⁇ t ⁇ XZ -h ⁇ ZZ
  • a 1X , a 1Y and a 1Z respectively represent measured outputs of the first three-axis accelerometer a 1 in the x-axis, y-axis and z-axis directions;
  • a 2X , a 2Y and a 2Z respectively represent measured outputs of the second three-axis accelerometer a 2 in the x-axis, y-axis and z-axis directions;
  • ⁇ is the angular velocity vector, Is the angular acceleration around the z-axis.
  • a gravity gradient measuring method comprising:
  • the first three-axis accelerometer and the second three-axis accelerometer are disposed in a chamber formed by the rotating table and the vacuum layer, and the rotating table rotates horizontally around the vertical axis, the first three-axis accelerometer and the second three-axis accelerometer Symmetrically set on the x-axis relative to the origin of the coordinate axis, the distance between the first three-axis accelerometer and the second three-axis accelerometer and the origin of the coordinate axis are both R; the first three-axis accelerometer and the second three-axis accelerometer are in z
  • the axis is symmetrically arranged with respect to the origin of the coordinate axis, and the distance between the first three-axis accelerometer and the second three-axis accelerometer on the z-axis is h.
  • the gravity gradient value ⁇ g is:
  • ⁇ XX , ⁇ YY and ⁇ ZZ are the gravity gradient values on the x-axis, y-axis and z-axis, respectively
  • ⁇ XY and ⁇ YX are the gravity gradient values on the x-axis and y-axis
  • ⁇ XY ⁇ YX
  • ⁇ XY ⁇ YX
  • ⁇ YZ ⁇ ZY .
  • determining the gravity gradient values under the nominal axis according to the acceleration values measured by the first three-axis accelerometer and the second three-axis accelerometer includes:
  • a 1Z -a 2Z 2R sin ⁇ t ⁇ YZ -2R cos ⁇ t ⁇ XZ -h ⁇ ZZ
  • a 1X , a 1Y and a 1Z respectively represent measured outputs of the first three-axis accelerometer a 1 in the x-axis, y-axis and z-axis directions;
  • a 2X , a 2Y and a 2Z respectively represent measured outputs of the second three-axis accelerometer a 2 in the x-axis, y-axis and z-axis directions;
  • ⁇ is the angular velocity vector, Is the angular acceleration around the z-axis.
  • FIG. 1 is a schematic view of an embodiment of a gravity gradient measuring device of the present invention.
  • FIG. 2 is a schematic view of another embodiment of the gravity gradient measuring device of the present invention.
  • FIG. 3 is a schematic view of still another embodiment of the gravity gradient measuring device of the present invention.
  • FIG. 4 is a schematic view showing the electrode structure of the three-axis accelerometer of the present invention.
  • FIG. 5 is a schematic diagram of an embodiment of a gravity gradient measuring method of the present invention.
  • 6a-6c are schematic diagrams showing the positional arrangement of a three-axis accelerometer according to the present invention.
  • Figure 7 is a schematic diagram showing the coordinate rotation of the three-axis accelerometer of the present invention.
  • the gravity gradient measuring apparatus includes a rotating table 3, a first three-axis accelerometer 1, a second three-axis accelerometer 2, a vacuum layer 7, and a measuring module, and the rotating table 3 is horizontally rotated about a vertical axis. among them:
  • the vacuum layer 7 is disposed on the rotary table 3, and the vacuum layer 7 and the rotary table 3 form a first chamber.
  • the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are located inside the first chamber, and the first three-axis acceleration
  • the meter 1 and the second three-axis accelerometer 2 are symmetrically arranged on the x-axis with respect to the origin of the coordinate axis, and the distance between the first three-axis accelerometer 1 and the second three-axis accelerometer 2 and the origin of the coordinate axis are both R;
  • the axis accelerometer 1 and the second three-axis accelerometer 2 are symmetrically arranged on the z-axis with respect to the origin of the coordinate axis, and the distance between the first three-axis accelerometer and the second three-axis accelerometer on the z-axis is h.
  • a measurement module determines the gravity gradient value under the coordinate axis using the acceleration values measured by the first three-axis accelerometer 1 and the second three-axis accelerometer 2.
  • the gravity gradient measuring device further comprises a first magnetic shielding layer 4 disposed inside the first chamber to shield the external magnetic field interference.
  • the first magnetic shield layer 4 and the rotary table 3 form a second chamber, and the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are located inside the second chamber.
  • the first magnetic shield layer 4 is a metal shield layer.
  • the material of the first magnetic shield layer 4 is permalloy.
  • the gravity gradient measuring device may further include a second magnetic shielding layer 5 disposed inside the second chamber for shielding external AC magnetic field interference.
  • the second magnetic shielding layer 5 and the rotating table 3 form a third chamber, and the first three-axis accelerometer 1 and the second three-axis accelerometer 2 are located inside the third chamber.
  • the second magnetic shield layer 5 is a superconducting shield layer.
  • the material of the second magnetic shield layer 5 is a YBCO (yttrium copper oxide) metal compound.
  • the gravity gradient measuring device further includes a liquid nitrogen layer 6 disposed between the first magnetic shield layer 4 and the second magnetic shield layer 5.
  • the liquid nitrogen layer 6 is a cylindrical double-layer metal thin wall layer filled with liquid nitrogen in the metal thin wall layer.
  • the liquid nitrogen layer 6 internally realizes a low temperature environment of about 77K, so that the YBCO metal compound of the second magnetic shield layer 5 realizes a superconducting state, and can reduce the Brownian motion heat of the first three-axis accelerometer 1 and the second three-axis accelerometer 2 Noise Sound, improve the measurement accuracy of the two first three-axis accelerometer 1 and the second three-axis accelerometer 2.
  • the first three-axis accelerometer 1 and the second three-axis accelerometer 2 have the same structure and function, and each includes a floating mass body 8, an x-axis acceleration detecting electrode 9, a y-axis acceleration detecting electrode 10, and The z-axis acceleration detecting electrode 11.
  • the x-axis acceleration detecting electrode 9, the y-axis acceleration detecting electrode 10, and the z-axis acceleration detecting electrode 11 constitute a cavity in which the suspended mass body 8 is suspended.
  • the suspended mass 8 is a hollow outer surface closed metal mass.
  • the x-axis acceleration detecting electrode 9, the y-axis acceleration detecting electrode 10, and the z-axis acceleration detecting electrode 11 measure the x-axis, y-axis, and z-axis acceleration values through a high-frequency capacitor bridge circuit, that is, by measuring the suspended mass body 8 on the x-axis.
  • the displacement in the y-axis and z-axis directions is then derived by displacement versus time second order to obtain the magnitude of the acceleration in the x-axis, y-axis, and z-axis directions.
  • the x-axis acceleration detecting electrode 9, the y-axis acceleration detecting electrode 10, and the z-axis acceleration detecting electrode 11 realize electrostatic force support and feedback control in three orthogonal directions of the x-axis, the y-axis, and the z-axis by the low-frequency voltage applied to the electrodes The function.
  • ⁇ XX , ⁇ YY and ⁇ ZZ are the gravity gradient values on the x-axis, y-axis and z-axis, respectively
  • ⁇ XY and ⁇ YX are the gravity gradient values on the x-axis and y-axis
  • ⁇ XY ⁇ YX
  • ⁇ XY ⁇ YX
  • ⁇ YZ ⁇ ZY .
  • the measurement module can utilize formulas
  • a 1Z -a 2Z 2R sin ⁇ t ⁇ YZ -2R cos ⁇ t ⁇ XZ -h ⁇ ZZ
  • a 1X , a 1Y and a 1Z respectively represent measured outputs of the first three-axis accelerometer a 1 in the x-axis, y-axis and z-axis directions;
  • a 2X , a 2Y and a 2Z respectively represent measured outputs of the second three-axis accelerometer a 2 in the x-axis, y-axis and z-axis directions;
  • ⁇ is the angular velocity vector, Is the angular acceleration around the z-axis.
  • FIG. 5 is a schematic diagram of an embodiment of a gravity gradient measuring method of the present invention.
  • the gravity gradient measuring method can be performed by using the gravity gradient measuring device according to any one of the embodiments of FIG. 1 to FIG. 4.
  • step 501 the acceleration value is measured by the first three-axis accelerometer and the second three-axis accelerometer.
  • Step 502 Determine a gravity gradient value under the coordinate axis according to the acceleration value.
  • the position coordinates of the first three-axis accelerometer 1 and the second three-axis accelerometer 2 on the x-axis, the y-axis, and the z-axis are (R, 0, h/2) and (- R, 0, -h/2).
  • An accelerometer coordinate system X a -O-Y a is determined , referred to as a coordinate system a, on the same plane as the coordinate system X-O-Y shown in Figures 6a-6c.
  • the coordinate system a and the three-axis accelerometer are fixed, that is, the accelerometer coordinate system X a -O-Y a is a coordinate system that rotates together with the three-axis accelerometer, and X g -O-Y g is the inertial coordinate. system.
  • FIG. 6a is a perspective view
  • FIG. 6b is a front view
  • FIG. 6c is a top view.
  • the difference between the linear accelerations in the x-axis, y-axis, and z-axis directions is measured by the first three-axis accelerometer 1 and the second three-axis accelerometer 2, respectively, and the x-axis measured by the two three-axis accelerometers.
  • the linear accelerations in the y-axis and z-axis directions are subtracted to obtain differential acceleration.
  • the relationship between differential acceleration and gravity gradient parameters is:
  • the gravity gradient matrix below, the gravity gradient matrix under the inertial system can be obtained by coordinate transformation; ⁇ a is the difference matrix of the first three-axis accelerometer 1 and the second three-axis accelerometer 2 position vector, when the first triaxial acceleration 1 meter and a second three-axis accelerometer when placed in position between the relatively fixed, this value is constant;
  • L a is a first three-axis accelerometer 1 second triaxial accelerometer 2 and the measured value obtained after the direct displacement difference Gradient measurement matrix, which also contains angular velocity components Angular acceleration component Rotating the acceleration for the inertial coordinate system, The angular velocity of the coordinate system a.
  • the angular velocity matrix is:
  • the angular acceleration matrix is:
  • the matrix of measured values of the first three-axis accelerometer 1 and the second three-axis accelerometer 2 in the coordinate system a can be written as:
  • the displacement difference matrix is:
  • L a is a gradient measurement matrix, which also includes angular velocity and angular acceleration components, and the angular velocity matrix (2) and the angular acceleration matrix (3) are substituted into the equation (6), and the expansion is obtained:
  • the separation signal of ° is processed to obtain four gravity gradient components.
  • X, Y, Z three coordinate axes of space
  • a 1 , a 2 respectively representing the first three-axis accelerometer 1 and the second three-axis accelerometer 2;
  • a 1X , a 1Y , a 1Z respectively represent the measured output of the first three-axis accelerometer a 1 in the X, Y, Z axis directions;
  • a 2X , a 2Y , a 2Z respectively represent the measured output of the second three-axis accelerometer a 2 in the X, Y, Z axis directions;
  • R the distance between the first three-axis accelerometer a 1 or the second three-axis accelerometer a 2 and the zero point O on the X-axis;
  • ⁇ angular velocity vector. Normally there are three components of ⁇ x , ⁇ y , ⁇ z , which are 3 ⁇ 1 vectors, but since the device of the present invention rotates only around the Z axis, its scalar value is equal to ⁇ z ;
  • the angular acceleration vector is the first derivative of the angular velocity ⁇ . Normally The three components are a 3 ⁇ 1 vector. Rotate at an ideal leveling speed,
  • ⁇ x , ⁇ y , ⁇ z angular velocities around the X axis, the Y axis, and the Z axis, respectively.
  • the direction of the angular velocity vector can be judged by the right hand rule, that is, the axis to be wound is the angular velocity direction of the axis;
  • X a -OY a is the coordinate system a on the XOY plane, is fixed together with the three-axis accelerometer, and rotates with the three-axis accelerometer, that is, the position vector of the three-axis accelerometer in the X a -OY a coordinate system constant;
  • X g -OY g is the coordinate system g on the XOY plane, which is inertial, stationary, and does not rotate with the three-axis accelerometer;
  • ⁇ a is the gravity gradient matrix under the coordinate system a of the three-axis accelerometer
  • ⁇ g is the coordinate system g system, that is, the gravity gradient matrix in the inertial coordinate system;
  • the position vector of the first three-axis accelerometer a 1 in the coordinate system a is
  • the position vector of the second three-axis accelerometer a 2 in the coordinate system a is
  • ⁇ a a difference matrix of position vectors of the first three-axis accelerometer a 1 and the second three-axis accelerometer a 2 in the coordinate system a;
  • ⁇ angular velocity matrix
  • the gradient measurement matrix is a quantity obtained by directly dividing the combination of the accelerometer output values by the displacement difference, and includes components such as angular velocity and angular acceleration. Because it is in the a coordinate system, the upper right corner is marked a;
  • the gradient measurement matrix L a contains:
  • the gradient matrix ⁇ a under the coordinate system a contains:
  • the gradient matrix ⁇ g under the inertial coordinate system g contains: Where ⁇ g is also the target gradient matrix of the measurement method of the present invention.
  • the measurement of the full tensor gravity gradient can be conveniently achieved by two three-axis accelerometers.

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Gyroscopes (AREA)

Abstract

一种重力梯度测量方法及装置，其中在重力梯度测量装置中，旋转台（3）绕地垂轴水平旋转，真空层（7）设置在旋转台（3）上，真空层（7）和旋转台（3）形成第一腔室，第一三轴加速度计（1）和第二三轴加速度计（2）位于第一腔室内部，第一三轴加速度计（1）和第二三轴加速度计（2）在x轴上相对于坐标轴原点对称设置，第一三轴加速度计（1）和第二三轴加速度计（2）与坐标轴原点的距离均为R，第一三轴加速度计（1）和第二三轴加速度计（2）在z轴上相对于坐标轴原点对称设置，第一三轴加速度计（1）和第二三轴加速度计（2）在z轴上的距离为h；测量模块利用第一三轴加速度计（1）和第二三轴加速度计（2）测量的加速度值确定在坐标轴下的重力梯度值。重力梯度测量方法及装置可通过两个三轴加速度计实现全张量重力梯度的测量。

Description

重力梯度测量方法及装置 技术领域
本发明涉及重力测量领域，特别是涉及一种重力梯度测量方法及装置。
背景技术
重力梯度测量始于油气普查勘探应用。在20年代的美国，扭秤重力梯度仪是油气勘探普查的唯一有效工具。但由于仪器笨重，效率低，梯度数据的解释方法研究又没有跟上，因此逐渐被淘汰。但是随着时代的发展，高精度重力梯度数据的重要性逐渐显现出来，无论是对于地质勘探，地球重力场的精细模拟还是高精度的惯性导航，它都不可或缺。
重力梯度测量是测量地球重力加速度随空间的变化。由于重力梯度是地球重力场的空间微分，反映重力沿空间不同方向的变化率，因此，重力梯度测量能够反映场源的细节。
重力梯度测量最主要的优点在于，重力梯度值或重力高次导数具有比重力本身更高的分辨率。因此通过利用测量重力势二阶导数的重力梯度仪实时测量重力梯度张量分量，就能获得更准确的重力值和垂直偏差，对空间科学、地球科学和地质科学等科学技术发展起着重要作用。
虽然惯性导航系统定位精度高，但缺陷在于误差随时间积累不断增加，必须定期重调。重力梯度测量是水下修正或限定无源自主惯导系统误差积累的一种重要方法，不仅可以提高惯性导航的精度，而且解决了水下导航的长期隐蔽性问题。
重力梯度是重力加速度矢量的一阶空间导数，即重力位的二阶空间导数，可由下述张量表示：
Figure PCTCN2017091107-appb-000001
其中，由于重力梯度场的对称性和无旋性，又有Γxy＝Γyx,Γxz＝Γzx,Γyz＝Γzy， Γxx+Γyy+Γzz＝0。因此，重力梯度的9个张量可以简化为5个独立张量。如果测量得到5个独立张量，即可得到所有9个张量信息，即全张量重力梯度测量。
目前重力梯度测量方法主要是基于差分加速度测量原理，通过加速度计之间观测量的不同组合能够获取引力梯度张量的各个分量。
全张量重力梯度测量需要三个正交旋转轴来分离梯度信号，特别是重力梯度的非对角分量；或者需要额外使用角加速度计来测量非对角分量以及消除动态误差，一般需要6个以上加速度计的组合结构，因此大大增加了仪器系统的复杂性和成本。
发明内容
本发明的目的是：提供一种重力梯度测量方法及装置，通过简单的方法就能够进行全张量重力梯度测量。
根据本发明的一个方面，提供一种重力梯度测量装置，包括旋转台、第一三轴加速度计、第二三轴加速度计、真空层和测量模块，旋转台绕地垂轴水平旋转，其中：
真空层设置在旋转台上，真空层和旋转台形成第一腔室，第一三轴加速度计和第二三轴加速度计位于第一腔室内部，第一三轴加速度计和第二三轴加速度计在x轴上相对于坐标轴原点对称设置，第一三轴加速度计和第二三轴加速度计与坐标轴原点的距离均为R；第一三轴加速度计和第二三轴加速度计在z轴上相对于坐标轴原点对称设置，第一三轴加速度计和第二三轴加速度计在z轴上的距离为h；
测量模块利用第一三轴加速度计和第二三轴加速度计所测量的加速度值确定在坐标轴下的重力梯度值。
可选地，重力梯度测量装置还包括：
设置在第一腔室内部内部的第一磁屏蔽层，以便于屏蔽外界地磁场干扰，第一磁屏蔽层和旋转台形成第二腔室，第一三轴加速度计和第二三轴加速度计位于在第二腔室内部。
可选地，第一磁屏蔽层为金属屏蔽层。
可选地，第一磁屏蔽层的材料为坡莫合金。
可选地，重力梯度测量装置还包括：
设置在第二腔室内部的第二磁屏蔽层，用于屏蔽外界交流磁场干扰，第二磁屏蔽层和旋转台形成第三腔室，第一三轴加速度计和第二三轴加速度计位于在第三腔室内部。
重力梯度测量装置，其特征在于，还包括：
设置在第一磁屏蔽层和第二磁屏蔽层之间的液氮层。
可选地，第二磁屏蔽层是超导屏蔽层。
可选地，第二磁屏蔽层的材料为YBCO金属化合物。
可选地，第一三轴加速度计和第二三轴加速度计均包括悬浮质量体、x轴加速度检测电极、y轴加速度检测电极和z轴加速度检测电极，x轴加速度检测电极、y轴加速度检测电极和z轴加速度检测电极组成空腔，悬浮质量体悬浮在空腔的中心位置；
x轴加速度检测电极、y轴加速度检测电极和z轴加速度检测电极测量悬浮质量体在x轴、y轴和z轴方向上的位移，进而根据所测量的位移得到x轴、y轴和z轴方向上的加速度。
可选地，悬浮质量体为一个空心的外表面封闭的金属质量体。
可选地，重力梯度值Γg为：
Figure PCTCN2017091107-appb-000002
其中，ΓXX、ΓYY和ΓZZ分别为在x轴、y轴和z轴上的重力梯度值，ΓXY和ΓYX为在x轴和y轴上的重力梯度值，ΓXZ和ΓZX为在x轴和z轴上的重力梯度值，ΓYZ和ΓZY为在y轴和z轴上的重力梯度值，且满足ΓXX+ΓYY+ΓZZ＝0，ΓXY＝ΓYX，ΓXY＝ΓYX，ΓYZ＝ΓZY。
可选地，测量模块具体利用公式
Figure PCTCN2017091107-appb-000003
Figure PCTCN2017091107-appb-000004
a1Z-a2Z＝2R sinωtΓYZ-2R cosωtΓXZ-hΓZZ
确定重力梯度的分量ΓYZ、ΓXZ、ΓXY和ΓXX-ΓYY；
利用公式ΓZZ＝(a2Z-a1Z)/h确定重力梯度的分量ΓZZ；
其中，a1X、a1Y和a1Z分别表示第一三轴加速度计a1在x轴、y轴和z轴方向上的测量输出；
a2X、a2Y和a2Z分别表示第二三轴加速度计a2在x轴、y轴和z轴方向上的测量输出；
ω为角速度矢量，
Figure PCTCN2017091107-appb-000005
为绕z轴的角加速度。
根据本发明的另一方面，提供一种重力梯度测量方法，包括：
利用第一三轴加速度计和第二三轴加速度计测量加速度值；
根据加速度值确定在坐标轴下的重力梯度值；
其中，第一三轴加速度计、第二三轴加速度计设置在由旋转台和真空层形成的腔室内，旋转台绕地垂轴水平旋转，第一三轴加速度计和第二三轴加速度计在x轴上相对于坐标轴原点对称设置，第一三轴加速度计和第二三轴加速度计与坐标轴原点的距离均为R；第一三轴加速度计和第二三轴加速度计在z轴上相对于坐标轴原点对称设置，第一三轴加速度计和第二三轴加速度计在z轴上的距离为h。
可选地，重力梯度值Γg为：
Figure PCTCN2017091107-appb-000006
其中，ΓXX、ΓYY和ΓZZ分别为在x轴、y轴和z轴上的重力梯度值，ΓXY和ΓYX为在x轴和y轴上的重力梯度值，ΓXZ和ΓZX为在x轴和z轴上的重力梯度值，ΓYZ和ΓZY为在y轴和z轴上的重力梯度值，且满足ΓXX+ΓYY+ΓZZ＝0，ΓXY＝ΓYX，ΓXY＝ΓYX，ΓYZ＝ΓZY。
可选地，根据第一三轴加速度计和第二三轴加速度计所测量的加速度值确定在标轴下的重力梯度值包括：
利用公式
Figure PCTCN2017091107-appb-000007
Figure PCTCN2017091107-appb-000008
a1Z-a2Z＝2R sinωtΓYZ-2R cosωtΓXZ-hΓZZ
确定重力梯度的分量ΓYZ、ΓXZ、ΓXY和ΓXX-ΓYY；
利用公式ΓZZ＝(a2Z-a1Z)/h确定重力梯度的分量ΓZZ；
其中，a1X、a1Y和a1Z分别表示第一三轴加速度计a1在x轴、y轴和z轴方向上的测量输出；
a2X、a2Y和a2Z分别表示第二三轴加速度计a2在x轴、y轴和z轴方向上的测量输出；
ω为角速度矢量，
Figure PCTCN2017091107-appb-000009
为绕z轴的角加速度。
通过以下参照附图对本发明的示例性实施例的详细描述，本发明的其它特征及其优点将会变得清楚。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明重力梯度测量装置一个实施例的示意图。
图2为本发明重力梯度测量装置另一实施例的示意图。
图3为本发明重力梯度测量装置又一实施例的示意图。
图4为本发明三轴加速度计电极结构示意图。
图5为本发明重力梯度测量方法一个实施例的示意图。
图6a-图6c为本发明三轴加速度计位置布置示意图。
图7为本发明三轴加速度计的坐标旋转示意图。
具体实施方式
以下结合附图和实施例对本发明做进一步的详细说明。
图1为本发明重力梯度测量装置一个实施例的示意图。如图1所示，重力梯度测量装置包括旋转台3、第一三轴加速度计1、第二三轴加速度计2、真空层7和测量模块，旋转台3绕地垂轴水平旋转。其中：
真空层7设置在旋转台3上，真空层7和旋转台3形成第一腔室，第一三轴加速度计1和第二三轴加速度计2位于第一腔室内部，第一三轴加速度计1和第二三轴加速度计2在x轴上相对于坐标轴原点对称设置，第一三轴加速度计1和第二三轴加速度计2与坐标轴原点的距离均为R；第一三轴加速度计1和第二三轴加速度计2在z轴上相对于坐标轴原点对称设置，第一三轴加速度计和第二三轴加速度计在z轴上的距离为h。
测量模块(图中未示出)利用第一三轴加速度计1和第二三轴加速度计2所测量的加速度值确定在所述坐标轴下的重力梯度值。
图2为本发明重力梯度测量装置另一实施例的示意图。与图1所示实施例相比，在图2所示实施例中，重力梯度测量装置还包括设置在第一腔室内部内部的第一磁屏蔽层4，以便于屏蔽外界地磁场干扰。其中第一磁屏蔽层4和旋转台3形成第二腔室，第一三轴加速度计1和第二三轴加速度计2位于在第二腔室内部。
可选地，第一磁屏蔽层4为金属屏蔽层。例如，第一磁屏蔽层4的材料为坡莫合金。
可选地，重力梯度测量装置还可包括设置在第二腔室内部的第二磁屏蔽层5，用于屏蔽外界交流磁场干扰。其中，第二磁屏蔽层5和旋转台3形成第三腔室，第一三轴加速度计1和第二三轴加速度计2位于在第三腔室内部。
可选地，第二磁屏蔽层5是超导屏蔽层。例如，第二磁屏蔽层5的材料为YBCO(氧化钇钡铜)金属化合物。
图3为本发明重力梯度测量装置又一实施例的示意图。与图2所示实施例相比，在图3所示实施例中，重力梯度测量装置还包括设置在第一磁屏蔽层4和第二磁屏蔽层5之间的液氮层6。
其中，液氮层6为一个圆柱形的双层金属薄壁层，在金属薄壁层内充入液体氮。液氮层6内部实现77K左右的低温环境，使得第二磁屏蔽层5的YBCO金属化合物实现超导态，并可以降低第一三轴加速度计1和第二三轴加速度计2的布朗运动热噪 声，提高两个第一三轴加速度计1和第二三轴加速度计2的测量精度。
如图4所示，上述的第一三轴加速度计1和第二三轴加速度计2的结构和功能相同，均包括悬浮质量体8、x轴加速度检测电极9、y轴加速度检测电极10和z轴加速度检测电极11。x轴加速度检测电极9、y轴加速度检测电极10和z轴加速度检测电极11组成一个空腔，悬浮质量体8悬浮在该空腔的中心位置。悬浮质量体8为一个空心的外表面封闭的金属质量体。x轴加速度检测电极9、y轴加速度检测电极10和z轴加速度检测电极11通过高频电容电桥电路实现x轴、y轴和z轴加速度值测量，即通过测量悬浮质量体8在x轴、y轴和z轴方向上的位移，然后通过位移对时间二阶求导以得到x轴、y轴和z轴方向上的加速度大小。x轴加速度检测电极9、y轴加速度检测电极10和z轴加速度检测电极11通过加载在电极上的低频电压实现x轴、y轴、z轴三个正交方向上的静电力支承和反馈控制的功能。
这里需要说明的是，重力梯度值Γg为：
Figure PCTCN2017091107-appb-000010
其中，ΓXX、ΓYY和ΓZZ分别为在x轴、y轴和z轴上的重力梯度值，ΓXY和ΓYX为在x轴和y轴上的重力梯度值，ΓXZ和ΓZX为在x轴和z轴上的重力梯度值，ΓYZ和ΓZY为在y轴和z轴上的重力梯度值，且满足ΓXX+ΓYY+ΓZZ＝0，ΓXY＝ΓYX，ΓXY＝ΓYX，ΓYZ＝ΓZY。
可选地，测量模块可利用公式
Figure PCTCN2017091107-appb-000011
Figure PCTCN2017091107-appb-000012
a1Z-a2Z＝2R sinωtΓYZ-2R cosωtΓXZ-hΓZZ
确定重力梯度的分量ΓYZ、ΓXZ、ΓXY和ΓXX-ΓYY；
利用公式ΓZZ＝(a2Z-a1Z)/h确定重力梯度的分量ΓZZ；
其中，a1X、a1Y和a1Z分别表示第一三轴加速度计a1在x轴、y轴和z轴方向上的测量输出；
a2X、a2Y和a2Z分别表示第二三轴加速度计a2在x轴、y轴和z轴方向上的测量输出；
ω为角速度矢量，
Figure PCTCN2017091107-appb-000013
为绕z轴的角加速度。
图5为本发明重力梯度测量方法一个实施例的示意图。其中，该重力梯度测量方法可采用图1-图4中任一实施例涉及的重力梯度测量装置来执行。
步骤501，利用第一三轴加速度计和第二三轴加速度计测量加速度值。
步骤502，根据加速度值确定在坐标轴下的重力梯度值。
下面对本发明所采用的测量方法进行具体说明。
如图6a-6c所示，第一三轴加速度计1和第二三轴加速度计2在x轴、y轴、z轴上的位置坐标分别为(R、0、h/2)和(-R、0、-h/2)。确定一个加速度计坐标系Xa－O－Ya，称之为坐标系a，与图6a-6c所示的坐标系X－O－Y在同一平面上。但是坐标系a与三轴加速度计是固连的，即加速度计坐标系Xa－O－Ya是随三轴加速度计一起旋转的坐标系，而Xg－O－Yg则是惯性坐标系。
其中，图6a为立体图，图6b为正视图，图6c为俯视图。
首先通过第一三轴加速度计1和第二三轴加速度计2分别在旋转过程中测量x轴、y轴、z轴方向上的直线加速度的差值，两个三轴加速度计测量的x轴、y轴、z轴方向上的直线加速度相减，得到差分加速度，差分加速度与重力梯度参数的关系公式为：
Figure PCTCN2017091107-appb-000014
其中，
Figure PCTCN2017091107-appb-000015
分别为第一三轴加速度计1和第二三轴加速度计2在坐标系a内的测量值矩阵；Γa为第一三轴加速度计1和第二三轴加速度计2所在坐标系a系下的重力梯度矩阵，惯性系下的重力梯度矩阵可以通过坐标变换得到；ρa为第一三轴加速度计1和第二三轴加速度计2位置向量的差值矩阵，当第一三轴加速度计1和第二三轴加速度计2之间的摆放位置相对固定时，此值恒定；La为第一三轴加速度计1和第二三轴加速度计2测量值对位移直接差分后得到的梯度测量矩阵，其中还包含角速度分量
Figure PCTCN2017091107-appb-000016
及角加速度分量
Figure PCTCN2017091107-appb-000017
Figure PCTCN2017091107-appb-000018
为惯性坐标系旋转加速度，
Figure PCTCN2017091107-appb-000019
为坐标系a的角速度。
其中，角速度矩阵为：
Figure PCTCN2017091107-appb-000020
角加速度矩阵为：
Figure PCTCN2017091107-appb-000021
第一三轴加速度计1和第二三轴加速度计2在坐标系a内的测量值矩阵可写为：
Figure PCTCN2017091107-appb-000022
位移差分矩阵为：
Figure PCTCN2017091107-appb-000023
将式(4)、式(5)代入公式(1)可得：
Figure PCTCN2017091107-appb-000024
La为梯度测量矩阵，其中还包含角速度及角加速度分量，将角速度矩阵(2)和角加速度矩阵(3)代入式(6)，展开可得：
Figure PCTCN2017091107-appb-000025
至此获得了第一三轴加速度计1和第二三轴加速度计2测量值与坐标系a下重力梯度的直接关系，以下进行坐标系a下的重力梯度与惯性系的坐标变换。
坐标系a下的重力梯度与惯性坐标系下的重力梯度的坐标变换公式为：
Figure PCTCN2017091107-appb-000026
如图7所示，当坐标系a绕惯性系旋转的角速度为ω时，变换公式的转置矩阵为：
Figure PCTCN2017091107-appb-000027
将式(9)代入式(8)，展开：
Figure PCTCN2017091107-appb-000028
由此得到：
Figure PCTCN2017091107-appb-000029
由于第一三轴加速度计1和第二三轴加速度计2仅绕Z轴旋转，因此ωx＝ωy＝0，ωz恒等于ω，故将式(10)代入式(7)，整理可得：
Figure PCTCN2017091107-appb-000030
ΓYZ、ΓXZ、ΓXY和ΓXX-ΓYY；
利用公式ΓZZ＝(a2Z-a1Z)/h确定重力梯度的分量ΓZZ；
对式(11)所示的所得结果进行分析，可以发现四个重力梯度分量ΓYZ、ΓXZ、ΓXY和ΓXX-ΓYY分别附加在第一三轴加速度计1和第二三轴加速度计2输出组合中的旋转角速度ω的一倍频和二倍频正交信号上。首先通过对四个重力梯度分量ΓYZ、ΓXZ、ΓXY和ΓXX-ΓYY分频，分别获取一倍频和二倍频的信号，再分别通过放大解调，并经由两个相差90°的检波信号分离处理，即可获得4个重力梯度分量。而ΓZZ分量为Z轴方向上的加速度输出组合a1Z-a2Z中的常量，即ΓZZ＝(a2z-a1z)/h，直接通过a1Z和a2Z输出值的差值进行滤波即可获得该分量。至此获得了重力梯度的全部5个独立分量，即可获得全张量重力梯度。
上述公式和表达式中的符号的含义为:
X,Y,Z：空间的三个坐标轴；
O：空间坐标轴的零点；
a1,a2：分别表示第一三轴加速度计1和第二三轴加速度计2；
a1X,a1Y,a1Z：分别表示第一三轴加速度计a1在X,Y,Z轴方向上的测量输出；
a2X,a2Y,a2Z：分别表示第二三轴加速度计a2在X,Y,Z轴方向上的测量输出；
R：第一三轴加速度计a1或第二三轴加速度计a2与零点O在X轴上的距离；
h：第一三轴加速度计a1和第二三轴加速度计a2在z轴上的距离；
ω：角速度矢量。正常有ωx,ωy,ωz三个分量，是个3×1向量，但是由于本发明所涉及的装置是只绕Z轴旋转，所以其标量值等于ωz；
Figure PCTCN2017091107-appb-000031
角加速度矢量，为角速度ω的一阶导。正常有
Figure PCTCN2017091107-appb-000032
三个分量，是个3×1向量。按理想匀角速度旋转，
Figure PCTCN2017091107-appb-000033
ωx,ωy,ωz：分别为绕X轴,Y轴,Z轴的角速度。角速度矢量的方向可通过右手定则判断，即所绕的轴为该轴的角速度方向；
Figure PCTCN2017091107-appb-000034
分别为绕X轴,Y轴,Z轴的角加速度；
Xa-O-Ya：为X-O-Y平面上的坐标系a，与三轴加速度计一起固连，并随三轴加速度计一起旋转，即三轴加速度计在Xa-O-Ya坐标系中的位置向量保持不变；
Xg-O-Yg：为X-O-Y平面上的坐标系g，为惯性系，静止，不与三轴加速度计一起旋转；
Figure PCTCN2017091107-appb-000035
分别为两个三轴加速度计在坐标系a内的测量值矩阵；
Γa：为三轴加速度计所在坐标系a系下的重力梯度矩阵；
Γg：为坐标系g系，即惯性坐标系下的重力梯度矩阵；
Figure PCTCN2017091107-appb-000036
第一三轴加速度计a1在坐标系a中的位置向量，为
Figure PCTCN2017091107-appb-000037
Figure PCTCN2017091107-appb-000038
第二三轴加速度计a2在坐标系a中的位置向量，为
Figure PCTCN2017091107-appb-000039
ρa：为第一三轴加速度计a1和第二三轴加速度计a2在坐标系a中位置向量的差值矩阵；
Ω：角速度矩阵；
Figure PCTCN2017091107-appb-000040
角加速度矩阵；
La：梯度测量矩阵，为加速度计输出值的组合直接除以位移差分得到的量，其中含有角速度及角加速度等分量。因为是在a坐标系下，所以右上角标a；
梯度测量矩阵La内包含：
Figure PCTCN2017091107-appb-000041
坐标系a下的梯度矩阵Γa包含：
Figure PCTCN2017091107-appb-000042
惯性坐标系g下的梯度矩阵Γg包含：
Figure PCTCN2017091107-appb-000043
其中Γg也是本发明测量方法的目标梯度矩阵。
通过实施本发明，可方便地通过两个三轴加速度计实现全张量重力梯度的测量。
本发明的描述是为了示例和描述起见而给出的，而并不是无遗漏的或者将本发明限于所公开的形式。很多修改和变化对于本领域的普通技术人员而言是显然的。选择和描述实施例是为了更好说明本发明的原理和实际应用，并且使本领域的普通技术人员能够理解本发明从而设计适于特定用途的带有各种修改的各种实施例。

Claims (15)

  1. 一种重力梯度测量装置，其特征在于，包括旋转台(3)、第一三轴加速度计(1)、第二三轴加速度计(2)、真空层(7)和测量模块，旋转台(3)绕地垂轴水平旋转，其中：
    真空层(7)设置在旋转台(3)上，真空层(7)和旋转台(3)形成第一腔室，第一三轴加速度计(1)和第二三轴加速度计(2)位于第一腔室内部，第一三轴加速度计(1)和第二三轴加速度计(2)在x轴上相对于坐标轴原点对称设置，第一三轴加速度计(1)和第二三轴加速度计(2)与坐标轴原点的距离均为R；第一三轴加速度计(1)和第二三轴加速度计(2)在z轴上相对于坐标轴原点对称设置，第一三轴加速度计和第二三轴加速度计在z轴上的距离为h；
    测量模块利用第一三轴加速度计(1)和第二三轴加速度计(2)所测量的加速度值确定在所述坐标轴下的重力梯度值。
  2. 根据权利要求1所述的重力梯度测量装置，其特征在于，还包括：
    设置在第一腔室内部内部的第一磁屏蔽层(4)，以便于屏蔽外界地磁场干扰，第一磁屏蔽层(4)和旋转台(3)形成第二腔室，第一三轴加速度计(1)和第二三轴加速度计(2)位于在第二腔室内部。
  3. 根据权利要求2所述的重力梯度测量装置，其特征在于，
    第一磁屏蔽层(4)为金属屏蔽层。
  4. 根据权利要求3所述的重力梯度测量装置，其特征在于，
    第一磁屏蔽层(4)的材料为坡莫合金。
  5. 根据权利要求2所述的重力梯度测量装置，其特征在于，还包括：
    设置在第二腔室内部的第二磁屏蔽层(5)，用于屏蔽外界交流磁场干扰，第二磁屏蔽层(5)和旋转台(3)形成第三腔室，第一三轴加速度计(1)和第二三轴加 速度计(2)位于在第三腔室内部。
  6. 根据权利要求5所述的重力梯度测量装置，其特征在于，还包括：
    设置在第一磁屏蔽层(4)和第二磁屏蔽层(5)之间的液氮层(6)。
  7. 根据权利要求6所述的重力梯度测量装置，其特征在于，
    第二磁屏蔽层(5)是超导屏蔽层。
  8. 根据权利要求7所述的重力梯度测量装置，其特征在于，
    第二磁屏蔽层(5)的材料为YBCO金属化合物。
  9. 根据权利要求1所述的重力梯度测量装置，其特征在于，
    第一三轴加速度计(1)和第二三轴加速度计(2)均包括悬浮质量体(8)、x轴加速度检测电极(9)、y轴加速度检测电极(10)和z轴加速度检测电极(11)，x轴加速度检测电极(9)、y轴加速度检测电极(10)和z轴加速度检测电极(11)组成空腔，悬浮质量体(8)悬浮在所述空腔的中心位置；
    x轴加速度检测电极(9)、y轴加速度检测电极(10)和z轴加速度检测电极(11)测量悬浮质量体(8)在x轴、y轴和z轴方向上的位移，进而根据所测量的位移得到x轴、y轴和z轴方向上的加速度。
  10. 根据权利要求9所述的重力梯度测量装置，其特征在于，
    悬浮质量体(8)为一个空心的外表面封闭的金属质量体。
  11. 根据权利要求1-10中任一项所述的重力梯度测量装置，其特征在于，重力梯度值Γg为：
    Figure PCTCN2017091107-appb-100001
    其中，ΓXX、ΓYY和ΓZZ分别为在x轴、y轴和z轴上的重力梯度值，ΓXY和ΓYX为在x轴和y轴上的重力梯度值，ΓXZ和ΓZX为在x轴和z轴上的重力梯度值，ΓYZ和ΓZY为在y轴和z轴上的重力梯度值，且满足ΓXX+ΓYY+ΓZZ＝0，ΓXY＝ΓYX，ΓXY＝ΓYX，ΓYZ＝ΓZY。
  12. 根据权利要求11所述的重力梯度测量装置，其特征在于，
    测量模块具体利用公式
    Figure PCTCN2017091107-appb-100002
    Figure PCTCN2017091107-appb-100003
    a1Z-a2Z＝2R sinωtΓYZ-2R cosωtΓXZ-hΓZZ
    确定重力梯度的分量ΓYZ、ΓXZ、ΓXY和ΓXX-ΓYY；
    利用公式ΓZZ＝(a2Z-a1Z)/h确定重力梯度的分量ΓZZ；
    其中，a1X、a1Y和a1Z分别表示第一三轴加速度计a1在x轴、y轴和z轴方向上的测量输出；
    a2X、a2Y和a2Z分别表示第二三轴加速度计a2在x轴、y轴和z轴方向上的测量输出；
    ω为角速度矢量，
    Figure PCTCN2017091107-appb-100004
    为绕z轴的角加速度。
  13. 一种重力梯度测量方法，其特征在于，包括：
    利用第一三轴加速度计(1)和第二三轴加速度计(2)测量加速度值；
    根据所述加速度值确定在所述坐标轴下的重力梯度值；
    其中，第一三轴加速度计(1)、第二三轴加速度计(2)设置在由旋转台(3)和真空层(7)形成的腔室内，旋转台(3)绕地垂轴水平旋转，第一三轴加速度计(1)和第二三轴加速度计(2)在x轴上相对于坐标轴原点对称设置，第一三轴加速度计(1)和第二三轴加速度计(2)与坐标轴原点的距离均为R；第一三轴加速度计(1)和第二三轴加速度计(2)在z轴上相对于坐标轴原点对称设置，第一三轴加速度计和第二三轴加速度计在z轴上的距离为h。
  14. 根据权利要求13所述的方法，其特征在于，
    重力梯度值Γg为：
    Figure PCTCN2017091107-appb-100005
    其中，ΓXX、ΓYY和ΓZZ分别为在x轴、y轴和z轴上的重力梯度值，ΓXY和ΓYX为在x轴和y轴上的重力梯度值，ΓXZ和ΓZX为在x轴和z轴上的重力梯度值，ΓYZ和ΓZY为在y轴和z轴上的重力梯度值，且满足ΓXX+ΓYY+ΓZZ＝0，ΓXY＝ΓYX，ΓXY＝ΓYX，ΓYZ＝ΓZY。
  15. 根据权利要求14所述的方法，其特征在于，
    根据第一三轴加速度计(1)和第二三轴加速度计(2)所测量的加速度值确定在所述坐标轴下的重力梯度值包括：
    利用公式
    Figure PCTCN2017091107-appb-100006
    Figure PCTCN2017091107-appb-100007
    a1Z-a2Z＝2R sinωtΓYZ-2R cosωtΓXZ-hΓZZ
    确定重力梯度的分量ΓYZ、ΓXZ、ΓXY和ΓXX-ΓYY；
    利用公式ΓZZ＝(a2Z-a1Z)/h确定重力梯度的分量ΓZZ；
    其中，a1X、a1Y和a1Z分别表示第一三轴加速度计a1在x轴、y轴和z轴方向上的测量输出；
    a2X、a2Y和a2Z分别表示第二三轴加速度计a2在x轴、y轴和z轴方向上的测量输出；
    ω为角速度矢量，
    Figure PCTCN2017091107-appb-100008
    为绕z轴的角加速度。
PCT/CN2017/091107 2017-03-09 2017-06-30 重力梯度测量方法及装置 WO2018161474A1 (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US16/491,936 US11402538B2 (en) 2017-03-09 2017-06-30 Gravity gradient measurement method and apparatus

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201710137515.X 2017-03-09
CN201710137515.XA CN107015287B (zh) 2017-03-09 2017-03-09 一种重力梯度测量装置及测量方法

Publications (1)

Publication Number Publication Date
WO2018161474A1 true WO2018161474A1 (zh) 2018-09-13

Family

ID=59439828

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2017/091107 WO2018161474A1 (zh) 2017-03-09 2017-06-30 重力梯度测量方法及装置

Country Status (3)

Country Link
US (1) US11402538B2 (zh)
CN (1) CN107015287B (zh)
WO (1) WO2018161474A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110967041A (zh) * 2019-12-18 2020-04-07 自然资源部国土卫星遥感应用中心 基于张量不变理论的卫星引力梯度数据精度的验证方法
CN114089441A (zh) * 2021-11-26 2022-02-25 华中科技大学 一种重力梯度仪测量系统数值仿真方法

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109062157B (zh) * 2018-08-09 2021-02-02 中电九天智能科技有限公司 无人电视机组装线的智能控制系统
CN109902695B (zh) * 2019-03-01 2022-12-20 辽宁工程技术大学 一种面向像对直线特征匹配的线特征矫正与提纯方法
CN110006633A (zh) * 2019-04-11 2019-07-12 南京聚科光电技术有限公司 一种降低激光相频噪声测试用光纤干涉仪热噪声的装置
CN111830589B (zh) * 2020-06-18 2022-07-26 中国船舶重工集团公司第七0七研究所 旋转加速度计重力梯度敏感器冗余测量信号离线融合方法
CN112213796B (zh) * 2020-09-17 2021-09-28 北京大学 一种绝对重力测量方法及绝对重力仪
CN112230295B (zh) * 2020-09-18 2021-11-09 北京大学 基于Sagnac效应角加速度计的重力梯度探测方法
CN112505793B (zh) * 2020-12-07 2021-11-23 北京大学 一种绝对重力测量装置
CN112698416B (zh) * 2020-12-16 2022-05-03 湖南航天机电设备与特种材料研究所 一种适用于捷联惯导系统的重力加速度测量方法及系统
CN113267821B (zh) * 2021-04-30 2022-05-03 北京大学 一种基于角运动的重力梯度测量方法及系统
CN113447994B (zh) * 2021-06-25 2022-04-12 中国船舶重工集团公司第七0七研究所 一种基于旋转加速度计原理的重力梯度全张量测量方法
CN113484921B (zh) * 2021-09-02 2021-12-24 华中光电技术研究所（中国船舶重工集团公司第七一七研究所） 一种四频双拉曼激光系统及冷原子水平重力梯度测量方法

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3926054A (en) * 1974-11-11 1975-12-16 Draper Lab Charles S Gravity gradiometer
US5357802A (en) * 1993-05-07 1994-10-25 Textron, Incorporated Rotating accelerometer gradiometer
US20060207326A1 (en) * 2005-03-18 2006-09-21 Moody Martin V Cross-component superconducting gravity gradiometer with improved linearity and sensitivity and method for gravity gradient sensing
US20090235740A1 (en) * 2008-03-18 2009-09-24 Symphony Acoustics, Inc. Gravity Gradient Sensor
CN201477216U (zh) * 2009-09-07 2010-05-19 北京航天控制仪器研究所 一种加速度计式相对重力测量仪
CN102099710A (zh) * 2008-09-25 2011-06-15 技术资源有限公司 用于检测重力梯度的检测器
CN102227651A (zh) * 2008-10-24 2011-10-26 M·V·穆迪 带有扭转挠曲枢轴的重力梯度仪
CN102608668A (zh) * 2011-12-19 2012-07-25 华中科技大学 重力梯度测量系统及其测量方法
CN104133252A (zh) * 2014-07-05 2014-11-05 中国科学院电工研究所 重力梯度测量辅助定位的超导磁悬浮定位装置及定位方法
CN104597520A (zh) * 2015-01-05 2015-05-06 中国船舶重工集团公司第七0七研究所 一种重力梯度仪重力梯度敏感器输入激励产生方法
CN105652334A (zh) * 2016-01-05 2016-06-08 华中科技大学 一种基于位移差分的mems重力梯度仪
CN105674972A (zh) * 2014-11-17 2016-06-15 中国航空工业第六八研究所 小型化组合单轴冷原子惯性传感器及其测量方法
CN105717553A (zh) * 2016-01-29 2016-06-29 东南大学 一种旋转加速度计重力梯度仪标定方法
CN105738962A (zh) * 2016-02-23 2016-07-06 武汉天琴引力科技有限公司 全张量重力梯度测量装置和测量方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3965753A (en) * 1970-06-01 1976-06-29 Browning Jr Alva Laroy Electrostatic accelerometer and/or gyroscope radioisotope field support device
EP0547742B1 (en) * 1991-12-19 1995-12-13 Motorola, Inc. Triaxial accelerometer
CA2690958C (en) * 2006-11-20 2015-06-30 Technological Resources Pty Limited A gravity gradiometer
US8347711B2 (en) * 2008-09-15 2013-01-08 Lockheed Martin Corporation Atom-interferometric, stepped gravity gradient measuring system
US8061048B2 (en) * 2008-09-29 2011-11-22 Schlumberger Technology Corporation Apparatus for azimuth measurements using gyro sensors
US8909496B2 (en) * 2011-02-17 2014-12-09 Gedex Inc. Gravity gradiometer and methods for measuring gravity gradients
GB201415087D0 (en) * 2014-08-26 2014-10-08 Univ Glasgow Measurement of acceleration
CN204613420U (zh) * 2015-05-04 2015-09-02 中国科学院测量与地球物理研究所 一种基于高压悬浮的便携式相对重力仪
CN105738963B (zh) * 2016-03-02 2018-08-03 中国科学院电工研究所 一种重力加速度测量装置
CN205594182U (zh) * 2016-03-14 2016-09-21 中国科学院测量与地球物理研究所 海洋重力仪稳定平台

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3926054A (en) * 1974-11-11 1975-12-16 Draper Lab Charles S Gravity gradiometer
US5357802A (en) * 1993-05-07 1994-10-25 Textron, Incorporated Rotating accelerometer gradiometer
US20060207326A1 (en) * 2005-03-18 2006-09-21 Moody Martin V Cross-component superconducting gravity gradiometer with improved linearity and sensitivity and method for gravity gradient sensing
US20090235740A1 (en) * 2008-03-18 2009-09-24 Symphony Acoustics, Inc. Gravity Gradient Sensor
CN102099710A (zh) * 2008-09-25 2011-06-15 技术资源有限公司 用于检测重力梯度的检测器
CN102227651A (zh) * 2008-10-24 2011-10-26 M·V·穆迪 带有扭转挠曲枢轴的重力梯度仪
CN201477216U (zh) * 2009-09-07 2010-05-19 北京航天控制仪器研究所 一种加速度计式相对重力测量仪
CN102608668A (zh) * 2011-12-19 2012-07-25 华中科技大学 重力梯度测量系统及其测量方法
CN104133252A (zh) * 2014-07-05 2014-11-05 中国科学院电工研究所 重力梯度测量辅助定位的超导磁悬浮定位装置及定位方法
CN105674972A (zh) * 2014-11-17 2016-06-15 中国航空工业第六八研究所 小型化组合单轴冷原子惯性传感器及其测量方法
CN104597520A (zh) * 2015-01-05 2015-05-06 中国船舶重工集团公司第七0七研究所 一种重力梯度仪重力梯度敏感器输入激励产生方法
CN105652334A (zh) * 2016-01-05 2016-06-08 华中科技大学 一种基于位移差分的mems重力梯度仪
CN105717553A (zh) * 2016-01-29 2016-06-29 东南大学 一种旋转加速度计重力梯度仪标定方法
CN105738962A (zh) * 2016-02-23 2016-07-06 武汉天琴引力科技有限公司 全张量重力梯度测量装置和测量方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110967041A (zh) * 2019-12-18 2020-04-07 自然资源部国土卫星遥感应用中心 基于张量不变理论的卫星引力梯度数据精度的验证方法
CN110967041B (zh) * 2019-12-18 2021-09-14 自然资源部国土卫星遥感应用中心 基于张量不变理论的卫星引力梯度数据精度的验证方法
CN114089441A (zh) * 2021-11-26 2022-02-25 华中科技大学 一种重力梯度仪测量系统数值仿真方法

Also Published As

Publication number Publication date
CN107015287A (zh) 2017-08-04
US20200073010A1 (en) 2020-03-05
US11402538B2 (en) 2022-08-02
CN107015287B (zh) 2018-12-11

Similar Documents

Publication Publication Date Title
WO2018161474A1 (zh) 重力梯度测量方法及装置
Moody et al. Three-axis superconducting gravity gradiometer for sensitive gravity experiments
US20230097313A1 (en) Electromagnetic gradiometers
US20160091620A1 (en) Triaxial accelerometer assembly and in-situ calibration method for improved geodetic and seismic measurements
CN109709628B (zh) 一种旋转加速度计重力梯度仪标定方法
Wu et al. A nano-g micromachined seismic sensor for levelling-free measurements
US8578771B2 (en) Compact gradiometer with accelerometer based rotational control
AU2010248152B2 (en) Improved gravity sensing instrument
US11175305B2 (en) Accelerometer based on diamagnetic levitation in a ring magnet
CN108445547A (zh) 一种三分量海洋重力磁力复合测量装置
CN108710001B (zh) 一种两轴一体陀螺加速度计及伺服控制方法
CN208207239U (zh) 一种三分量海洋重力磁力复合测量装置
CN108845365A (zh) 一种三分量海洋重力测量装置
Zhang et al. Cross-coupling coefficient estimation of a nano-g accelerometer by continuous rotation modulation on a tilted rate table
CN105738962B (zh) 全张量重力梯度测量装置和测量方法
Graizer The response to complex ground motions of seismometers with Galperin sensor configuration
CN112213796B (zh) 一种绝对重力测量方法及绝对重力仪
Fremd Symmetric triaxial rotational piezoelectric seismometers
Zhi et al. Digital fluxgate magnetometer for detection of microvibration
Tsuchiya et al. Dynamic sensitivity matrix measurement for single-mass SOI 3-axis accelerometer
Ali et al. Issues and Challenges of Sensor Technologies in Microelectromechanical System (MEMS) in Smartphones for Motion Tracking Applications
Paik et al. Development of a sensitive superconducting gravity gradiometer for geological and navigational applications
CN113267821B (zh) 一种基于角运动的重力梯度测量方法及系统
CN113466492B (zh) 一种测量加速度计模型参数磁矢量灵敏度的装置及方法
Kivikangas Estimating three-dimensional motion using IMU sensors

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17899622

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 17899622

Country of ref document: EP

Kind code of ref document: A1