CN112213796B - 一种绝对重力测量方法及绝对重力仪 - Google Patents

一种绝对重力测量方法及绝对重力仪 Download PDF

Info

Publication number
CN112213796B
CN112213796B CN202010981101.7A CN202010981101A CN112213796B CN 112213796 B CN112213796 B CN 112213796B CN 202010981101 A CN202010981101 A CN 202010981101A CN 112213796 B CN112213796 B CN 112213796B
Authority
CN
China
Prior art keywords
sensor
measuring element
translational motion
acceleration
motion
Prior art date
Legal status (The legal status 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 status listed.)
Active
Application number
CN202010981101.7A
Other languages
English (en)
Other versions
CN112213796A (zh
Inventor
陈彦钧
曾卫益
阳春霞
李正斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking University
Original Assignee
Peking University
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 Peking University filed Critical Peking University
Priority to CN202010981101.7A priority Critical patent/CN112213796B/zh
Publication of CN112213796A publication Critical patent/CN112213796A/zh
Application granted granted Critical
Publication of CN112213796B publication Critical patent/CN112213796B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种绝对重力测量方法及绝对重力仪。本发明的测量方法步骤包括:1)在待测位置设置一平移运动测量元件和一旋转运动测量元件;2)利用该平移运动测量元件测量得到在该待测位置因平移运动产生的相对位移d,利用旋转运动测量元件测量得到平移运动测量元件敏感轴因旋转运动产生的旋转角度θ;3)根据公式
Figure DDA0002687542040000011
确定该待测位置的绝对重力加速度g;其中f为观测频率。本方法可直接从天然地震波中观测绝对重力,因而无需主动或被动隔振装置,且可以实现地震、重力的同时观测,便于地球物理联合勘探的实施;在理论上,随着观测频率的降低,绝对重力的观测精度将依照平方关系增加,有利于提升绝对重力观测的精度。

Description

一种绝对重力测量方法及绝对重力仪
技术领域
本发明涉及一种绝对重力测量方法及绝对重力仪,可用于绝对重力加速度的观测范畴,如地球重力场监测、大地测量、地球动力学研究、密度体及异常监测监控、建筑抗震实验、矿产资源开发、地球物理勘探、地震及地质灾害预警等工程及科研领域。
背景技术
地球上的物体均会影响到重力加速度值(以下所称的“重力”均等价于“重力加速度”),据研究,地表的山川及大洋底部的海沟对重力的影响为小数点后第3位,而大型建筑物对重力的影响为第7位。因此,提高重力测量的精度,可极大提升重力对于周边物体的敏感性。
在重力测量方面,目前有两大类测量方式:绝对重力测量和相对重力测量,前者指直接测定重力加速度的真值,而后者则是测定某点相对于参考点的重力加速度的差值。在绝对重力测量方面,典型的重力仪为基于自由落体原理的美国劳雷公司FG5型、意大利计量院IMGC型和基于原子干涉的法国巴黎天文台可移动式重力仪,上述仪器最高精度约为小数点后第8位。这类仪器体积大、结构复杂、不易使用。在相对重力测量方面,典型的重力仪为基于弹簧结构的LCR重力仪,其原理简单,但精度低;采用超导是目前的主流研究方向,其噪声低、稳定性好,精度高,可到小数点后11位,但需要体积庞大的低温冷却设备。从激光干涉到冷原子干涉仪,绝对重力测量手段不断进步,但测量原理始终是自由落体原理,不利于、甚至限制了绝对重力仪的发展。
发明内容
针对目前绝对重力测量原理较为单一的问题,本发明的目的在于提供一种绝对重力测量方法,采用的是一种不同于自由落体原理的绝对重力测量原理。该原理在理论上无需目前高精度重力仪所需的低温冷却设备,降低了设备的体积;由于其可直接从天然地震波中观测绝对重力,因而无需主动或被动隔振装置,且可以实现地震、重力的同时观测,便于地球物理联合勘探的实施;在理论上,随着观测频率的降低,绝对重力的观测精度将依照平方关系增加,有利于提升绝对重力观测的精度。
本发明的技术方案为:
一种绝对重力测量方法,其步骤包括:
1)在待测位置设置一平移运动测量元件和一旋转运动测量元件;
2)利用该平移运动测量元件测量得到在该待测位置因平移运动产生的相对位移d,利用旋转运动测量元件测量得到平移运动测量元件敏感轴因旋转运动产生的旋转角度θ;
3)根据公式
Figure BDA0002687542020000021
确定该待测位置的绝对重力加速度g;其中f为观测频率。
进一步的,步骤1)所述平移运动测量元件为位移传感器或速度传感器或加速度传感器或加加速度传感器,步骤1)所述旋转运动测量元件为角度传感器或角速度传感器或角加速度传感器;步骤2)中采用位移传感器或速度传感器或加速度传感器或加加速度传感器测量平移运动时,分别采集在该待测位置因平移运动产生的相对位移d或相对运动速度v或加速度a或加加速度b;步骤2)中采用角度传感器或角速度传感器或角加速度传感器测量旋转运动时,分别采集平移运动测量元件敏感轴因旋转运动产生的旋转角度θ或角速度ω或角加速度β;步骤3)中根据公式
Figure BDA0002687542020000022
v=j2πfd、a=(j2πf)2d、b=(j2πf)3d、ω=j2πfθ、β=(j2πf)2θ确定该待测位置的绝对重力加速度g;其中f为观测频率。
进一步的,所述平移运动测量元件适用的理论模型为质量-弹簧-活塞模型。
一种绝对重力仪,其特征在于,包括一平移运动测量元件、一旋转运动测量元件和数据处理单元;其中
该平移运动测量元件,用于测量得到其在待测位置因平移运动产生的相对位移d;
旋转运动测量元件,用于测量得到平移运动测量元件敏感轴因旋转运动产生的旋转角度θ;
数据处理单元,用于根据公式
Figure BDA0002687542020000024
确定该待测位置的绝对重力加速度g;其中f为观测频率。
进一步的,所述平移运动测量元件为位移传感器或速度传感器或加速度传感器或加加速度传感器,旋转运动测量元件为角度传感器或角速度传感器或角加速度传感器;当采用位移传感器或速度传感器或加速度传感器或加加速度传感器测量平移运动时,分别采集在该待测位置因平移运动产生的相对位移d或相对运动速度v或加速度a或加加速度b;当采用角度传感器或角速度传感器或角加速度传感器测量旋转运动时,分别采集平移运动测量元件敏感轴因旋转运动产生的旋转角度θ或角速度ω或角加速度β;然后根据公式
Figure BDA0002687542020000023
v=j2πfd、a=(j2πf)2d、b=(j2πf)3d、ω=j2πfθ、β=(j2πf)2θ确定该待测位置的绝对重力加速度g;其中f为观测频率。
附图2展示了该绝对重力测量方法的基本原理,其中,有一固定在待测位置(如地面)处的平移运动传感器(类型包括位移传感器、速度传感器、加速度传感器、加加速度传感器等等),平移运动传感器示意图如附图1所示。由于实际平移运动传感器系统均具有弹性阻尼和粘滞阻尼,为不失一般性,在理论上采用质量-弹簧-活塞模型描述具有以上阻尼的平移运动传感器:
Figure BDA0002687542020000031
其中:H(s)表示平移运动传感器在拉普拉斯域中的单位加速度的冲激位移响应;s为拉普拉斯复频率;k为弹簧系数,由附图1所示的弹簧所代表的弹性阻尼力产生;c为阻尼系数,由附图1所示的活塞所代表的粘滞阻尼力产生;m为附图1中质量块的质量。
在待测位置(如地面),同时存在附图2所示方向(即平移运动传感器底面与xoz面相交方向)的平移运动d和法向为y方向的旋转运动θ。其中,平移运动d可由天然存在的平移地震波、人工或各类设备激发的平移地震波等等原因产生,旋转运动θ可由天然存在的旋转地震波、人工或各类设备激发的旋转地震波等等导致。不妨不失一般性地假设平移运动的单位为位移(国际单位制下的单位为米)、旋转运动的单位为角度(国际单位制下的单位为弧度)。在该平移运动d的激励下,平移运动传感器内部质量块产生相对运动,记为δ,其方向与平移运动d的方向一致,如附图2所示。则在拉普拉斯域中,仅由于平移运动d产生的相对位移响应ψ为:
ψ(s)=s2d·H(s) (2)
由于平移运动传感器固定于地面,平移运动传感器受到该旋转运动影响,也将产生旋转运动,使得平移运动传感器与水平面(xoy面)的夹角为θ。则在拉普拉斯域中,平移运动传感器由于平移运动d和旋转运动θ共同产生的响应δ为:
δ(s)=[s2d+gθ]·H(s) (3)
其中,g为重力加速度。
由式(3)可得:
Figure BDA0002687542020000032
其中用到了s=j2πf,f为观测频率。
式(4)说明:本发明给出的这一绝对重力测量原理是在频率域中观测重力,而目前的重力测量原理均为时域测量。本原理可从平移及旋转运动中观测绝对重力,即绝对重力加速度g可通过测量平移运动d和旋转运动θ获得,因而无需主动或被动隔振装置,具备地震、重力同时观测的理论基础,为地球物理联合勘探的实施提供了一种高效的理论。另外,当旋转运动传感器与平移运动传感器的精度一定时,降低观测频率f可提升绝对重力的探测精度,且是依照平方关系增长的,这有利于提升绝对重力观测的精度。
需要指出的是,这里为描述的方便,以位移为单位描述平移运动,以角度为单位描述旋转运动。但事实上,无论平移运动单位为速度、加速度还是加加速度,旋转运动单位为角速度或角加速度,都可以通过求导数获得本质与式(4)相同的测量式,仍属于本发明所提出的这一绝对重力测量原理的派生。
具体地,由于在时间域求一次导数相当于在拉普拉斯域乘s=j2πf(其中j为单位虚数根),故当平移运动的单位变为速度v时,在拉普拉斯域中与前述位移量纲的平移运动d的关系为:
v=j2πfd (5)
用式(5)替换式(4)中的平移运动d项即可得到平移运动单位为速度、旋转运动单位为角度时的重力观测方程。
依此类推,记平移运动单位为加速度、加加速度时,拉普拉斯域的表示分别为a、b,则:
a=(j2πf)2d (6)
b=(j2πf)3d (7)
记旋转运动单位为角速度、角加速度时,拉普拉斯域的表示分别为ω、β,则:
ω=j2πfθ (8)
β=(j2πf)2θ (9)
根据式(6)至式(9),结合式(4)即可得到不同单位表示的平移运动、旋转运动测定绝对重力的观测方程。
本发明的实施方式如附图3的框图所示。如前述,本发明提供的绝对重力测量原理基于平移运动和旋转运动的测量,而后再经由式(4)或其派生式(结合式(5)至式(9))获得绝对重力。目前平移运动测量元件包括位移传感器(即敏感量为以位移为单位的平移运动)、速度传感器(即敏感量为以速度为单位的平移运动)、加速度传感器(即敏感量为以加速度为单位的平移运动)、加加速度传感器(即敏感量为以加加速度为单位的平移运动),旋转运动测量元件包括角度传感器(即敏感量为以角度为单位的旋转运动)、角速度(即敏感量为以角速度为单位的旋转运动)、角加速度(即敏感量为以角加速度为单位的旋转运动)。
与现有技术相比,本发明的积极效果为:
本发明提供了一种不同于自由落体原理的绝对重力测量原理,进行重力测量时无需目前高精度重力仪所需的低温冷却设备,降低了设备的体积;由于其可直接从天然地震波中观测绝对重力,因而无需主动或被动隔振装置,且可以实现地震、重力的同时观测,便于地球物理联合勘探的实施;而且随着观测频率的降低,绝对重力的观测精度将依照平方关系增加,有利于提升绝对重力观测的精度。
附图说明
图1为平移运动传感器示意图,其中质量块、弹簧、活塞及其外部黑色方块所示的框架整体组成了一个平移运动传感器(类型包括位移传感器、速度传感器、加速度传感器、加加速度传感器等等),用质量-弹簧-活塞模型表示。该传感器固定于地面。
图2为绝对重力观测原理图,该平移运动传感器的固定位置同时存在旋转运动θ及平移运动d,δ表示平移运动传感器内部质量块的相对位移响应。通过测量旋转运动θ及平移运动d即可通过式(4)得到重力加速度。
图3为绝对重力观测原理实施框图,首先通过角度传感器、角速度传感器或者角加速度传感器获得相应旋转运动量,通过位移传感器、速度传感器、加速度传感器或者加加速度传感器获得相应平移运动量,最后按式(4)-式(9)获得绝对重力加速度值。
图4为绝对重力观测原理实施示例,将旋转运动传感器固定于框架侧面,平移运动传感器固定于底面,整体框架固定于地面上。
具体实施方式
下面结合附图对本发明进行进一步详细描述。
基于本绝对重力测量原理的传感器代表示例如附图4所示。其中,平移运动传感器(类型包括位移传感器、速度传感器、加速度传感器、加加速度传感器等等)、旋转运动传感器(类别包括角度传感器、角速度传感器、角加速度传感器等等)被固定于固定框架上,固定框架被固定在地面。在地面存在天然存在的平移地震波、人工或各类设备激发的平移地震波等等即框架受到的平移运动来源,天然存在的旋转地震波、人工或各类设备激发的旋转地震波等等即框架受到的旋转运动来源。如附图3所示,通过角度传感器、角速度传感器或者角加速度传感器获得相应旋转运动量,通过位移传感器、速度传感器、加速度传感器或者加加速度传感器获得相应平移运动量,最后按式(4)-式(9)获得绝对重力加速度值。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种绝对重力测量方法,其步骤包括:
1)在待测位置设置一平移运动测量元件和一旋转运动测量元件;
2)利用该平移运动测量元件测量得到在该待测位置因平移运动产生的相对位移d,利用旋转运动测量元件测量得到平移运动测量元件敏感轴因旋转运动产生的旋转角度θ;
3)根据公式
Figure FDA0002687542010000011
确定该待测位置的绝对重力加速度g;其中f为观测频率。
2.如权利要求1所述的方法,其特征在于,步骤1)所述平移运动测量元件为位移传感器或速度传感器或加速度传感器或加加速度传感器,步骤1)所述旋转运动测量元件为角度传感器或角速度传感器或角加速度传感器;步骤2)中采用位移传感器或速度传感器或加速度传感器或加加速度传感器测量平移运动时,分别采集在该待测位置因平移运动产生的相对位移d或相对运动速度v或加速度a或加加速度b;步骤2)中采用角度传感器或角速度传感器或角加速度传感器测量旋转运动时,分别采集平移运动测量元件敏感轴因旋转运动产生的旋转角度θ或角速度ω或角加速度β;步骤3)中根据公式
Figure FDA0002687542010000012
v=j2πfd、a=(j2πf)2d、b=(j2πf)3d、ω=j2πfθ、β=(j2πf)2θ确定该待测位置的绝对重力加速度g;其中f为观测频率。
3.如权利要求1所述的方法,其特征在于,所述平移运动测量元件适用的理论模型为质量-弹簧-活塞模型。
4.一种绝对重力仪,其特征在于,包括一平移运动测量元件、一旋转运动测量元件和数据处理单元;其中
该平移运动测量元件,用于测量得到其在待测位置因平移运动产生的相对位移d;
旋转运动测量元件,用于测量得到平移运动测量元件敏感轴因旋转运动产生的旋转角度θ;
数据处理单元,用于根据公式
Figure FDA0002687542010000013
确定该待测位置的绝对重力加速度g;其中f为观测频率。
5.如权利要求4所述的绝对重力仪,其特征在于,所述平移运动测量元件为位移传感器或速度传感器或加速度传感器或加加速度传感器,旋转运动测量元件为角度传感器或角速度传感器或角加速度传感器;当采用位移传感器或速度传感器或加速度传感器或加加速度传感器测量平移运动时,分别采集在该待测位置因平移运动产生的相对位移d或相对运动速度v或加速度a或加加速度b;当采用角度传感器或角速度传感器或角加速度传感器测量旋转运动时,分别采集平移运动测量元件敏感轴因旋转运动产生的旋转角度θ或角速度ω或角加速度β;然后根据公式
Figure FDA0002687542010000021
v=j2πfd、a=(j2πf)2d、b=(j2πf)3d、ω=j2πfθ、β=(j2πf)2θ确定该待测位置的绝对重力加速度g;其中f为观测频率。
CN202010981101.7A 2020-09-17 2020-09-17 一种绝对重力测量方法及绝对重力仪 Active CN112213796B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010981101.7A CN112213796B (zh) 2020-09-17 2020-09-17 一种绝对重力测量方法及绝对重力仪

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010981101.7A CN112213796B (zh) 2020-09-17 2020-09-17 一种绝对重力测量方法及绝对重力仪

Publications (2)

Publication Number Publication Date
CN112213796A CN112213796A (zh) 2021-01-12
CN112213796B true CN112213796B (zh) 2021-09-28

Family

ID=74049966

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010981101.7A Active CN112213796B (zh) 2020-09-17 2020-09-17 一种绝对重力测量方法及绝对重力仪

Country Status (1)

Country Link
CN (1) CN112213796B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113672002B (zh) * 2021-08-23 2022-03-22 九江学院 一种基于名义模型的冷原子重力仪主动隔振控制方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101639541A (zh) * 2009-09-07 2010-02-03 北京航天控制仪器研究所 一种加速度计式相对重力测量仪
CN102301257A (zh) * 2008-11-28 2011-12-28 格雷维泰克仪器有限公司 重力梯度仪
CN102608668A (zh) * 2011-12-19 2012-07-25 华中科技大学 重力梯度测量系统及其测量方法
WO2014115848A1 (ja) * 2013-01-25 2014-07-31 マグネデザイン株式会社 回転情報演算方法、回転情報演算プログラム、磁気型ジャイロスコープおよび移動体
CN107831551A (zh) * 2017-11-28 2018-03-23 济南大学 一种测量重力加速度的装置和方法
CN110850497A (zh) * 2019-11-21 2020-02-28 中国船舶重工集团公司第七0七研究所 基于原子干涉效应的绝对重力仪、陀螺仪敏感器及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107015287B (zh) * 2017-03-09 2018-12-11 中国科学院电工研究所 一种重力梯度测量装置及测量方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102301257A (zh) * 2008-11-28 2011-12-28 格雷维泰克仪器有限公司 重力梯度仪
CN101639541A (zh) * 2009-09-07 2010-02-03 北京航天控制仪器研究所 一种加速度计式相对重力测量仪
CN102608668A (zh) * 2011-12-19 2012-07-25 华中科技大学 重力梯度测量系统及其测量方法
WO2014115848A1 (ja) * 2013-01-25 2014-07-31 マグネデザイン株式会社 回転情報演算方法、回転情報演算プログラム、磁気型ジャイロスコープおよび移動体
CN107831551A (zh) * 2017-11-28 2018-03-23 济南大学 一种测量重力加速度的装置和方法
CN110850497A (zh) * 2019-11-21 2020-02-28 中国船舶重工集团公司第七0七研究所 基于原子干涉效应的绝对重力仪、陀螺仪敏感器及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Vertical deformation and absolute gravity;Ming Fang 等;《Geophys》;20011231;第146卷;第539-548页 *
世界重力梯度仪的研究现状;吴琼 等;《物探与化探》;20131031;第37卷(第5期);第761-768页 *

Also Published As

Publication number Publication date
CN112213796A (zh) 2021-01-12

Similar Documents

Publication Publication Date Title
Agnew Strainmeters and tiltmeters
US9645267B2 (en) Triaxial accelerometer assembly and in-situ calibration method for improved geodetic and seismic measurements
CN107015287B (zh) 一种重力梯度测量装置及测量方法
Jekeli A review of gravity gradiometer survey system data analyses
US3926054A (en) Gravity gradiometer
Evstifeev The state of the art in the development of onboard gravity gradiometers
CN109001829B (zh) 一种捷联式水下动态重力测量仪
Chapin Gravity instruments: Past, present, future
Jekeli Airborne gradiometry error analysis
CN113885098A (zh) 一种重力敏感器低频频率响应误差在线建模及补偿方法
CN112213796B (zh) 一种绝对重力测量方法及绝对重力仪
CN110068876A (zh) 基于载体自振动航空重力梯度仪运动误差补偿方法
Santos et al. Understanding the fundamentals of earthquake signal sensing networks
Wyatt et al. A comparison of tiltmeters at tidal frequencies
CN112649889A (zh) 一种六分量地震数据及绝对重力测量仪、测量方法
Wooldridge Review of modern airborne gravity focusing on results from GT-1A surveys
Pang et al. Low-cost IMU error intercorrection method for verticality measurement
Korobiichuk et al. Design of the novel double-ring dynamical gravimeter
CN104111063B (zh) 一种基于磁场的无线三维倾角传感器及其检测方法
CN106917621B (zh) 小孔径单陀螺水平井旋转定向测斜装置及方法
Tkachuk et al. Single-rotor integrating gyroscopic gravimeter
CN210488007U (zh) 三分量重力梯度场、三分量磁力梯度场获取装置
CN211741611U (zh) 基于减震超材料的井中三分量重力磁力复合仪器系统
Harrison Tilt observations in the Poorman mine near Boulder, Colorado
Chen et al. Surveying method of measurement while drilling based on the inertial sensor

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant