CN113311495B - 一种间接测量重力的方法及系统 - Google Patents
一种间接测量重力的方法及系统 Download PDFInfo
- Publication number
- CN113311495B CN113311495B CN202110599628.8A CN202110599628A CN113311495B CN 113311495 B CN113311495 B CN 113311495B CN 202110599628 A CN202110599628 A CN 202110599628A CN 113311495 B CN113311495 B CN 113311495B
- Authority
- CN
- China
- Prior art keywords
- gravity
- acceleration
- change
- earth
- indirectly measuring
- 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
Links
- 230000005484 gravity Effects 0.000 title claims abstract description 111
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000001133 acceleration Effects 0.000 claims abstract description 86
- 230000008859 change Effects 0.000 claims abstract description 44
- 239000013598 vector Substances 0.000 claims abstract description 31
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 8
- 238000012544 monitoring process Methods 0.000 claims abstract description 8
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 8
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 20
- 230000006870 function Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 230000005358 geomagnetic field Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 6
- 238000013213 extrapolation Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 238000000691 measurement method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 2
- 229910052753 mercury Inorganic materials 0.000 description 2
- 101100460704 Aspergillus sp. (strain MF297-2) notI gene Proteins 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 229910052792 caesium Inorganic materials 0.000 description 1
- TVFDJXOCXUVLDH-UHFFFAOYSA-N caesium atom Chemical compound [Cs] TVFDJXOCXUVLDH-UHFFFAOYSA-N 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
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
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明提供一种间接测量重力的方法及系统，基于天体产生的万有引力、地球产生的万有引力以及其他惯性力影响的合成，导致待测位置的重力加速度发生变化，通过定期监测待测位置的重力加速度方向变化，反推得到该待测位置的重力测量结果；监测重力加速度的方向变化时，测量各个时刻重力的方向，得到重力加速度的方向变化观测数据；根据自身坐标求得得当前位置的地球引起的加速度矢量的概略值作为初始解，计算重力加速度方向变化估测数据，结合重力加速度方向变化观测数据，进行迭代线性化求解得到重力测量结果。本发明方案实施简单方便，实用性强，解决了相关技术存在的实用性低及实际应用不便的问题，实施时只需通过倾斜仪即可测量出重力方向的变化。
Description
技术领域
本发明涉及测量技术领域，具体涉及一种间接测量重力的技术方案。
背景技术
重力测量在大地测量、地球物理、资源勘查、海洋研究以及军事方面有着重要地位。目前的重力测量方法主要有激光干涉、原子干涉、弹簧、超导等技术，其精度不断在提高，技术也由最早的自由落体法向着量子、超导等方向拓展。但其结构复杂，使用条件苛刻。
发明内容
为了解决上述现有技术中存在的不足，本发明提出了一种间接测量重力的技术方案。
为了实现上述目的，本发明提出一种间接测量重力的方法，基于天体产生的万有引力、地球产生的万有引力以及其他惯性力影响的合成，导致待测位置的重力加速度发生变化，通过定期监测待测位置的重力加速度方向变化，反推得到该待测位置的重力测量结果；监测重力加速度的方向变化时，测量各个时刻t1,t2,…,tk+1重力的方向θ1,θ2,…,θk+1，得到重力加速度的方向变化观测数据Δθi＝θi+1–θi，i＝1…k；设根据坐标
时刻、相关天体位置，能够求出在
处各个时刻t1,t2,…,tk+1各个力的合成形成的加速度矢量为V1,V2,…,Vk+1，则加速度的方向变化估测数据ΔVi＝Vi+1-Vi，i＝1…k，k为不小于3的整数，设当前的地球引起的万有引力造成的加速度为三维矢量
建立如下的重力加速度的方向变化相应方程组，
其中，f()是方向的差值关于
的函数；
反推得到待测位置的重力测量结果实现方式为，根据自身坐标
求得得当前位置的地球引起的加速度矢量的概略值
作为初始解，计算重力加速度方向变化估测数据，结合重力加速度方向变化观测数据，通过重力加速度的方向变化相应方程组，进行迭代线性化求解直至收敛到误差小于误差限，得到重力测量结果。
而且，监测重力加速度的方向变化采用倾斜仪测量实现。
而且，定期监测待测位置的重力加速度变化时，采用原子钟获取当前时刻。
而且，所述自身坐标采用根据惯性定位设备或者地磁场定位设备得到的当前的概略位置。
而且，所述自身坐标未知时，通过观测7次以上的数据，将自身位置和当前重力矢量一共6个未知数求解。
而且，所述天体包括太阳和月亮。
而且，计算重力加速度方向变化估测数据，实现方式如下，
根据自身位置、地球自转速度计算地球自转造成的加速度矢量
根据时刻、月球质量、月球坐标、万有引力常数、地心坐标、自身位置计算在地月系统中的加速度aem及其方向vem；
根据时刻、太阳质量、太阳坐标、万有引力常数、地心坐标、自身位置计算在日地系统中的加速度ase及其方向vse；
设将以上得到的加速度矢量与
进行矢量合成，形成四个或以上时刻的加速度矢量
k为不小于3的整数；
根据
计算加速度方向的夹角ΔVi,ΔVi＝Vi+1-Vi，得到加速度的方向变化估测数据。
另一方面，本发明还提供一种间接测量重力的系统，用于实现如上所述的一种间接测量重力的方法。
而且，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种间接测量重力的方法。
或者，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如上所述的一种间接测量重力的方法。
本发明利用了天体产生的万有引力与地球产生的万有引力以及其他惯性力的关系不断在变化，导致该处的重力方向发生极其微弱的变化，力的方向就是加速度的方向，加速度的大小和力的大小成正比。因此测量加速度与测量力是等价的。而每一个位置上随着时间的变化都是不同的，通过定期监测该位置的重力加速度方向变化，即可反推出重力。实施时只需通过倾斜仪即可测量出重力方向的变化。
本发明方案实施简单方便，实用性强，解决了相关技术存在的实用性低及实际应用不便的问题，能够提高用户体验，具有重要的市场价值。
具体实施方式
以下结合实施例具体说明本发明的技术方案。
地球重力其实是万有引力和惯性力的合力。本发明即是利用了天体产生的万有引力与地球产生的万有引力以及其他惯性力的关系不断在变化，导致该处的重力方向发生极其微弱的变化，力的方向就是加速度的方向，加速度的大小和力的大小成正比。因此测量加速度与测量力是等价的。而每一个位置上随着时间的变化都是不同的，通过定期监测该位置的重力加速度方向变化，即可反推出重力。通过倾斜仪即可测量出重力方向的变化。
设在地球某处，其在地心地固坐标系下的坐标为
那么在任意时刻，其所受到的重力加速度是由这么几种力的矢量合成而形成的：地球、月球、太阳以及其他星球对其的万有引力，地球自转形成的离心力，地月系统中的离心力，日地系统中的离心力，以及其他不能忽略的物体(例如附近的重物)对其的万有引力。按照短期恒定不变与始终在变化对这些力分为两大类：
第一类，短期恒定不变的：地球形成的万有引力，地球自转形成的离心力，通常来说，附近的重物短期内也不会有很大的变化；
第二类：始终在随着时间变化的：月球、太阳以及其他星球形成的万有引力；地月系统以及日地系统中的离心力。
对于第二类力，虽然随着时间在变化，但是其运行规律早已经被人们熟知，只需要知道时间，就可以准确计算其大小和方向。
理论上来说，只需要精确知道时间，又知道自身位置，就可以精确计算重力的方向和大小，然而由于地球的引力实际上并不能直接由公式准确算得，因为地球是一个不均匀的物体，内部结构十分复杂，而这些因素对地球引力的影响是巨大的，不能使用万有引力公式直接计算其引力。但是，通过重力方向的变化来反推重力却是可行的。这样既可以不被水平方向所干扰，又能通过差的方法减弱甚至消除掉大多数其他因素的误差，例如重力异常等。
从目前掌握的知识，地球周围有很多天体，但是因为这些天体的质量以及和地球之间的距离不同，他们对重力的影响是不同的。影响最大的是月球，其次是太阳，再其次是金星、火星、木星、水星和土星。由于目前的仪器水平尚不能测量出金星等星体对地球的万有引力的影响，因此在此暂不予考虑。待以后仪器的精度水平提高后，再进行考虑，从而进一步提高观测精度。
通过计算，有如下结果：
每间隔1秒钟，月球对地面物体造成的重力方向的影响变化量为10-6角秒量级，而如果要准确提取其变化量，需要的倾斜仪的精度至少需要达到10-10角秒。
每间隔1分钟，月球对地面物体造成的重力方向的影响变化量为10-5角秒量级，而如果要准确提取其变化量，需要的倾斜仪的精度至少需要达到10-7角秒。
每间隔10分钟，月球对地面物体造成的重力方向的影响变化量为10-4角秒量级，而如果要准确提取其变化量，需要的倾斜仪的精度至少需要达到10-5角秒。
每间隔60分钟，月球对地面物体造成的重力方向的影响变化量为10-3角秒量级，而如果要准确提取其变化量，需要的倾斜仪的精度至少需要达到10-4角秒。
太阳造成的重力方向变化大约为月球的1/3。以上考虑了日地、地月系统为非惯性系。其他星球造成的影响更小，其中金星、火星、木星、水星和土星等行星的影响如果考虑则精度更高。
目前，现有技术已经能够已经有能力对这些天体的运行轨迹进行精确计算，只需要知道准确的时刻即可。而目前精密时钟技术也非常先进，现在的铯钟，精度可以轻易达到10-14/5Day的水平，可以保证1年的误差不超过100ps。对于倾斜的观测可以达到10-5角秒的精度。通过GNSS等手段，也可以进行精确定位。基于以上设备，已经具备了根据万有引力方向变化计算重力的条件。显然，倾斜仪的精度越高，则间隔观测时间可以越短，计算的周期也就越短，结果的精度也越高。从目前的技术来说，倾斜仪精度的提高要远比提高重力仪精度容易，且成本也低得多。例如电容式倾斜仪，只需要极板面积100平方毫米左右，摆长50毫米，间距0.25毫米，其大小如同保温杯，即可达到10-5角秒的精度。光学法倾斜仪的也在相似水平。只是在此以前，没有更高精度的倾斜仪这种需求。具体实施时，可以采用相应的定制倾斜仪实现本发明技术方案。因此可以预见，本发明所述的定位方法可能会促进倾斜仪市场产品的精度进一步提高。
本发明提出的测量方法为：
1)设定当前的地球引起的万有引力造成的加速度为三维矢量
短时间内，地球引起的万有引力的方向和大小变化很小，可以认为是恒定的。在一定的时间间隔，在至少4个时刻t1,t2,t3,t4观测倾斜仪的角度θ1,θ2,θ3,θ4…，然后将相邻的两个观测值作差，可以得到三个以上观测值差Δθ1,Δθ2,Δθ3…Δθ1＝θ2-θ1,
Δθ2＝θ3-θ2,
Δθ3＝θ4-θ3
…
2)再根据
时刻、星体位置，可以准确求出各个时刻前面提到的除了地球引力造成的加速度之外的各个加速度的三维矢量
(i＝1…k,k+1，k为不小于3的整数)，将其与而在这些时刻地球引起的万有引力造成的加速度做矢量合成，即得到该处的重力加速度矢量
(i＝1…k,k+1，k为不小于3的整数)。计算相邻两个重力加速度的夹角，即可以求得3个以上计算方向差ΔV1,ΔV2,ΔV3…
ΔV1＝V2-V1,
ΔV2＝V3-V2,
ΔV3＝V4-V3,
…
3)如果
就是倾斜仪所在位置的准确的加速度，那么Δθ1，Δθ2，Δθ3…应该分别等于ΔV1，ΔV2，ΔV3…而事实上，
是需要求的未知数。由于ΔV1，ΔV2，ΔV3…是
以及时间的函数，因此以上的计算角度差的算式可以表达为：
其中，t1，t2，t3，t4…均为已知的精确观测值，而ΔV1，ΔV2，ΔV3…是已知的观测值Δθ1，Δθ2，Δθ3…，f()是方向的差值关于
的函数。因此上式成为了一个有三个未知数的方程组，包含了至少三个公式，因此不出意外的话，其有唯一解。该方程是非线性方程，可以通过泰勒展开、迭代求解。得到地球引力形成的加速度
后，将其与其他的力形成的加速度进行矢量合成，即为该处的重力加速度。
为了能够缩短收敛时间，并提高定位的成功率，可以用
计算该处地球万有引力形成的加速度作为初始解。
为了达到更精确的解，可以增加观测次数，以得到超定方程组，然后使用最小二乘等方法求解。
即使不知道自身位置，也可以通过增加观测次数到至少7次，从而得到6个以上的方程，将自身位置也当成未知数一并计算得到。但这样无疑增加了观测时间以及计算量。
以上涉及的计算方法均为成熟的方法，具体过程不再赘述。
为便于实施参考起见，提供实施例提出的一种万有引力定位方法实现具体如下：
(为简化表述，该实施例仅考虑月球、太阳造成的引力影响，并假设可以得到当前的精确位置。如果要考虑其它因素，只需要仿照该例加入相应的条件即可。如果无法得到当前位置，则需要至少重复2，3，4步骤7次(即观测7次以上的数据)，并且同时将自身位置和当前重力矢量一共6个未知数求解)
1，所需的设备：高精密时钟(一般用原子钟)，计算设备(可采用计算机或者其他带处理器和存储器的设备)，精密倾斜仪，定位设备(如GNSS或者地磁场定位设备)；
2，从高精密时钟获取当前时刻(年月日时分秒)，根据当前时刻计算月球、太阳等天体的位置；由于最终要计算自身在地球上的位置，因此一般使用地心地固坐标系；该计算方法是既有方法，不再详述；
3，测量当前重力的方向(倾斜仪的倾斜角度)；
4，间隔一段时间；
5，从定位设备获取当前位置
6，重复步骤2，3，4一共执行至少4个循环，从而得到至少4组数据：
观测时刻：t1，t2，t3，t4…
月球坐标：
太阳坐标：
倾斜角度：θ1，θ2，θ3，θ4…以及当前位置
(如果有定位设备如GNSS，就可以观测得到位置)
7，由于经过了航行、晃动等过程，不能保证倾斜仪观测到的角度是绝对倾斜角度，因此需要根据倾斜角度θ1，θ2，θ3，θ4…计算倾斜角度的变化量Δθi＝θi+1–θi,(i＝1…k，k为不小于3的整数)，从而得到至少3个倾斜角度变化量Δθ1，Δθ2，Δθ3…，即重力加速度的方向变化观测数据；
8，由
以及观测时刻ti，月球坐标
太阳坐标
(i＝1…k，k为不小于3的整数)，可以求得在
处，每个观测时刻的地球引力大小Ae和方向Ve，月球引力大小Am和方向Vm，太阳引力大小As和方向Vs，以及地球自转、地月系统中、日地系统中的离心力的大小和方向。由这些可以进一步计算合力的大小和方向。这个合力的大小和方向应该与观测值一致；
9，根据
求得当前位置的地球引起的加速度矢量的概略值
此即为初始解；
10，设所处位置由地球引力形成的加速度矢量为
根据初始解
以及步骤6中得到的天体坐标数据(月球坐标和太阳坐标)，计算4个或以上时刻的万有引力(同时也是加速度)的大小和方向分别求差值(即重力加速度的方向变化估测数据)，然后基于步骤7得到的重力加速度的方向变化观测数据，根据设定的误差限，进行迭代，直至收敛到误差小于误差限，即可求得
从而可以进一步求得该处的重力加速度。具体实现过程如下：
a)根据自身位置、地球自转速度计算地球自转造成的加速度矢量
b)根据时刻、月球质量、月球坐标、万有引力常数、地心坐标、自身位置计算在地月系统中的加速度aem及其方向vem；
c)根据时刻、太阳质量、太阳坐标、万有引力常数、地心坐标、自身位置计算在日地系统中的加速度ase及其方向vse；
d)设将以上a,b,c三个步骤得到的三个加速度矢量与
进行矢量合成，形成四个或以上时刻的加速度矢量
(i＝1…k,k+1，k为不小于3的整数)；
e)根据
计算加速度方向的夹角ΔVi,ΔVi＝Vi+1-Vi，(i＝1…k，k为不小于3的整数)，即加速度的方向变化估测数据：
ΔV1＝V2-V1,
ΔV2＝V3-V2,
ΔV3＝V4-V3…
f)其中b,c两个步骤均是计算以时间为自变量的函数，而d步骤是计算以
和时间为自变量的函数，其中
是未知量，时间是已知量。而步骤e的结果是基于b,c,d四个步骤的，因此ΔVi实为
的函数，因此有如下的方程组：
其中，f()是方向的差值关于
的函数。
残差ωi＝ΔVi-Δθi(i＝1…k，k为不小于3的整数)，只有
是观测点处地球引力的准确值的时候，残差才为0，Δθi均为已知的观测量。因此f()也是残差ωi关于
的函数。
在
处，将方程1进行泰勒展开实施线性化，可以得到：
其中，Gv是对应于f()的雅各比矩阵，
是当前解与上一次解(第一次迭代则为初始解)的变化量，Δ(ΔV)为当前解与上一次的解(第一次迭代则为初始解)加速度方向值的变化量，将ωi代入方程2，即可得到
而当前解
此时已经将地球引力值得解进行了更新，
比
更接近观测点实际的地球引力。如此迭代，直到
小于预定限值则结束迭代，将当前解作为最终解，完成测量。
通常，k值会大于3，从而得到超定方程组，提高结果的可靠性，进一步以
作为初始解，通过最小二乘等方法求解，从而得到当前位置地球引力形成的加速度矢量
进而与步骤a,b,c中的其他力形成的加速度合成，即可得到该处的重力加速度。
具体实施时，本发明技术方案提出的方法可由本领域技术人员采用计算机软件技术实现自动运行流程，实现方法的系统装置例如存储本发明技术方案相应计算机程序的计算机可读存储介质以及包括运行相应计算机程序的计算机设备，也应当在本发明的保护范围内。
在一些可能的实施例中，提供一种万有引力定位系统，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种万有引力定位方法。
在一些可能的实施例中，提供一种万有引力定位系统，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如上所述的一种万有引力定位方法。
在一些可能的实施例中,可以使用多台倾斜仪做差分以消除某些误差，得到更精确和可靠的结果。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (10)
1.一种间接测量重力的方法，其特征在于：基于天体产生的万有引力、地球产生的万有引力以及其他惯性力影响的合成，导致待测位置的重力加速度发生变化，通过定期监测待测位置的重力加速度方向变化，反推得到该待测位置的重力测量结果；
监测重力加速度的方向变化时，测量各个时刻t1,t2,…,tk+1重力的方向θ1,θ2,…,θk+1，得到重力加速度的方向变化观测数据Δθi＝θi+1–θi，i＝1…k；设根据坐标
时刻、相关天体位置，能够求出在
处各个时刻t1,t2,…,tk+1各个力的合成形成的加速度矢量为V1,V2,…,Vk+1，则加速度的方向变化估测数据ΔVi＝Vi+1-Vi，i＝1…k，k为不小于3的整数，设当前的地球引起的万有引力造成的加速度为三维矢量
建立如下的重力加速度的方向变化相应方程组，
其中，f()是方向的差值关于
的函数；
反推得到待测位置的重力测量结果实现方式为，根据自身坐标
求得当前位置的地球引起的加速度矢量的概略值
作为初始解，计算重力加速度方向变化估测数据，结合重力加速度方向变化观测数据，通过重力加速度的方向变化相应方程组，进行迭代线性化求解直至收敛到误差小于误差限，得到重力测量结果。
2.根据权利要求1所述间接测量重力的方法，其特征在于：监测重力加速度的方向变化时，采用倾斜仪实现测量重力的方向。
3.根据权利要求1所述间接测量重力的方法，其特征在于：定期监测待测位置的重力加速度变化时，采用原子钟获取当前时刻。
4.根据权利要求1所述间接测量重力的方法，其特征在于：所述自身坐标采用根据惯性定位设备或者地磁场定位设备得到的当前的概略位置。
5.根据权利要求1所述间接测量重力的方法，其特征在于：所述自身坐标未知时，通过观测7次以上的数据，将自身位置和当前重力矢量一共6个未知数求解。
6.根据权利要求1或2或3或4或5所述间接测量重力的方法，其特征在于：所述天体包括太阳和月亮。
7.根据权利要求6所述间接测量重力的方法，其特征在于：计算重力加速度方向变化估测数据，实现方式如下，
根据自身位置、地球自转速度计算地球自转造成的加速度矢量
根据时刻、月球质量、月球坐标、万有引力常数、地心坐标、自身位置计算在地月系统中的加速度aem及其方向vem；
根据时刻、太阳质量、太阳坐标、万有引力常数、地心坐标、自身位置计算在日地系统中的加速度ase及其方向vse；
设将以上得到的加速度矢量与
进行矢量合成，形成四个或以上时刻的加速度矢量
k为不小于3的整数；
根据
计算加速度方向的夹角ΔVi,ΔVi＝Vi+1-Vi，得到加速度的方向变化估测数据。
8.一种间接测量重力的系统，其特征在于：用于实现如权利要求1-7任一项所述的一种间接测量重力的方法。
9.根据权利要求8所述间接测量重力的系统，其特征在于：包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如权利要求1-7任一项所述的一种间接测量重力的方法。
10.根据权利要求8所述间接测量重力的系统，其特征在于：包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如权利要求1-7任一项所述的一种间接测量重力的方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110599628.8A CN113311495B (zh) | 2021-05-31 | 2021-05-31 | 一种间接测量重力的方法及系统 |
| PCT/CN2022/096070 WO2022253198A1 (zh) | 2021-05-31 | 2022-05-30 | 一种间接测量重力的方法及系统 |
| US18/523,846 US20240142660A1 (en) | 2021-05-31 | 2023-11-29 | Method and system for indirect measurement of gravity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110599628.8A CN113311495B (zh) | 2021-05-31 | 2021-05-31 | 一种间接测量重力的方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113311495A CN113311495A (zh) | 2021-08-27 |
| CN113311495B true CN113311495B (zh) | 2022-07-19 |
Family
ID=77376365
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110599628.8A Active CN113311495B (zh) | 2021-05-31 | 2021-05-31 | 一种间接测量重力的方法及系统 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20240142660A1 (zh) |
| CN (1) | CN113311495B (zh) |
| WO (1) | WO2022253198A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113311495B (zh) * | 2021-05-31 | 2022-07-19 | 武汉大学 | 一种间接测量重力的方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU94040423A (ru) * | 1994-10-28 | 1996-12-27 | С.Г. Комаров | Способ измерения скорости распространения в пространстве гравитационных сил |
| CN101639541A (zh) * | 2009-09-07 | 2010-02-03 | 北京航天控制仪器研究所 | 一种加速度计式相对重力测量仪 |
| CN102323624A (zh) * | 2011-08-05 | 2012-01-18 | 清华大学 | 绝对重力测量系统、测量方法及自由落体下落方法 |
| CN208672814U (zh) * | 2018-09-19 | 2019-03-29 | 南风（上海）精密物理仪器有限公司 | 一种相对重力加速度勘测仪 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2600493C2 (ru) * | 2011-09-07 | 2016-10-20 | Бейкер Хьюз Инкорпорейтед | Устройство и способ определения геологических границ |
| US9939551B2 (en) * | 2012-09-24 | 2018-04-10 | Schlumberger Technology Corporation | Systems, devices and methods for borehole gravimetry |
| CN108267792B (zh) * | 2018-04-13 | 2019-07-12 | 武汉大学 | 全球重力场模型反演方法 |
| CN111366984B (zh) * | 2020-03-23 | 2022-10-14 | 东华理工大学 | 一种基于重力卫星星间激光测距系统确定引力场模型的方法 |
| CN113311495B (zh) * | 2021-05-31 | 2022-07-19 | 武汉大学 | 一种间接测量重力的方法及系统 |
-
2021
- 2021-05-31 CN CN202110599628.8A patent/CN113311495B/zh active Active
-
2022
- 2022-05-30 WO PCT/CN2022/096070 patent/WO2022253198A1/zh active Application Filing
-
2023
- 2023-11-29 US US18/523,846 patent/US20240142660A1/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU94040423A (ru) * | 1994-10-28 | 1996-12-27 | С.Г. Комаров | Способ измерения скорости распространения в пространстве гравитационных сил |
| CN101639541A (zh) * | 2009-09-07 | 2010-02-03 | 北京航天控制仪器研究所 | 一种加速度计式相对重力测量仪 |
| CN102323624A (zh) * | 2011-08-05 | 2012-01-18 | 清华大学 | 绝对重力测量系统、测量方法及自由落体下落方法 |
| CN208672814U (zh) * | 2018-09-19 | 2019-03-29 | 南风（上海）精密物理仪器有限公司 | 一种相对重力加速度勘测仪 |
Non-Patent Citations (3)
| Title |
|---|
| GNSS信号录制与回放设备的研制及应用;张涛;《全球定位系统》;20180630;第43卷(第3期);第45-50页 * |
| 基于非线性最小二乘估计的eLoran/INS/磁传感器组合导航方法;闫泓衫 等;《中国惯性技术学报》;20210228;第29卷(第1期);第84-89页 * |
| 谈重力加速度g的几种测量方法;赵歆玥;《数理化学习》;20161231(第3期);第43-44页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022253198A1 (zh) | 2022-12-08 |
| CN113311495A (zh) | 2021-08-27 |
| US20240142660A1 (en) | 2024-05-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108051866B (zh) | 基于捷联惯性/gps组合辅助水平角运动隔离的重力测量方法 | |
| Blewitt et al. | Geodetic observations and global reference frame contributions to understanding sea-level rise and variability | |
| CN105737828B (zh) | 一种基于强跟踪的相关熵扩展卡尔曼滤波的组合导航方法 | |
| Richeson | Gravity gradiometer aided inertial navigation within non-GNSS environments | |
| CN110221332A (zh) | 一种车载gnss/ins组合导航的动态杆臂误差估计和补偿方法 | |
| CN107270893A (zh) | 面向不动产测量的杆臂、时间不同步误差估计与补偿方法 | |
| CN112327379B (zh) | 一种全张量重力梯度动态测量系统及方法 | |
| CN109470241B (zh) | 一种具备重力扰动自主补偿功能的惯性导航系统及方法 | |
| CN102736118B (zh) | 一种用于全球重力场测量的综合型卫星系统 | |
| CN110017812B (zh) | 一种垂线偏差、重力加速度和重力梯度的测量装置及方法 | |
| CN111947653A (zh) | 一种月面巡视探测器双模式惯性/视觉/天文导航方法 | |
| CN104833375B (zh) | 一种借助星敏感器的imu两位置对准方法 | |
| CN112325886B (zh) | 一种基于重力梯度仪和陀螺仪组合的航天器自主定姿系统 | |
| CN113311495B (zh) | 一种间接测量重力的方法及系统 | |
| Kwon et al. | Gravity requirements for compensation of ultra-precise inertial navigation | |
| CN109084755B (zh) | 一种基于重力视速度与参数辨识的加速度计零偏估计方法 | |
| CN109029499B (zh) | 一种基于重力视运动模型的加速度计零偏迭代寻优估计方法 | |
| CN103256932A (zh) | 一种替换结合外推的着陆导航方法 | |
| Peshekhonov | Problem of the vertical deflection in high-precision inertial navigation | |
| CN110514201A (zh) | 一种惯性导航系统及适用于高转速旋转体的导航方法 | |
| CN113310486B (zh) | 一种万有引力定位方法及系统 | |
| CN112882118B (zh) | 地固坐标系下动基座重力矢量估计方法、系统及存储介质 | |
| CN114705215A (zh) | 一种捷联惯导系统的纬度估计方法 | |
| CN113551669A (zh) | 基于短基线的组合导航定位方法及装置 | |
| Colombatti et al. | Huygens probe entry dynamic model and accelerometer data analysis |
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 |