JP4316777B2

JP4316777B2 - Gravity measuring device and method

Info

Publication number: JP4316777B2
Application number: JP2000225568A
Authority: JP
Inventors: 爾朗瀬川; 貞臣佐久間; 悟史吉田
Original assignee: Tokyo Keiki Inc
Current assignee: Tokyo Keiki Inc
Priority date: 2000-07-26
Filing date: 2000-07-26
Publication date: 2009-08-19
Anticipated expiration: 2020-07-26
Also published as: JP2002040155A

Description

【０００１】
【発明の属する技術分野】
本発明は地球表面における重力を検出するための重力測定装置及び方法に関し、特に、航行体に搭載される重力測定装置及びそれを使用する重力測定方法に関する。
【０００２】
【従来の技術】
地球表面における重力の値は、地下の岩石の密度構造、地心距離の変化等の理由により、測定地によって僅かに変動する。重力の測定値と標準重力の差は重力異常と呼ばれる。重力の測定は、この重力異常の測定に外ならない。重力を測定する場合、平地では携帯型の重力測定装置が用いられるが、極地方、砂漠地帯、山岳地帯、海洋等では車両、航空機等の航行体又は移動体に搭載した重力測定装置が使用される。
【０００３】
航行体又は移動体に搭載した重力測定装置によって重力を測定する場合、測定点の位置を精密に測定するために、ＤＧＰＳ（differential global positioning system）が使用される。
【０００４】
ＤＧＰＳはＧＰＳの測定精度を改善するために開発されたものである。ＧＰＳでは移動局（利用者）のみによって測地されるのに対して、ＤＧＰＳでは移動局（利用者）と基準局（固定局）によって測地される。ＤＧＰＳでは、予め正確な位置が測定されている基準局において、ＧＰＳ衛星からの擬似距離、時刻情報及び軌道データより補正値を演算し、それを移動局に送信する。移動局は、基準局からの補正値を用いてＧＰＳ衛星からの信号の補正を行う。
【０００５】
ＤＧＰＳにはコード位相測位法と搬送波位相測位法の２つの方法がある。コード位相測位法は測定値としてコード位相を用い、リアルタイム測位のために使用される。精度は１〜１０ｍ程度である。搬送波位相測位法は測定値として搬送波位相を用い、高精度の位置測定に使用する。精度はｍｍ〜ｃｍ程度である。
【０００６】
図４を参照して従来の重力測定装置及び方法の例を説明する。重力測定装置２０はヘリコプタ等の航行体１０に搭載される。ＤＧＰＳは、少なくとも４つのＧＰＳ衛星（図４では２つのＧＰＳ衛星１２、１３のみ図示）と航行体１０に搭載された２つの移動局（利用者）１４、１５と地上の２つの基準局（固定局）１６、１７を使用する。
【０００７】
第１の移動局（利用者）１４及び基準局（固定局）１６は搬送波位相測位用の受信装置であり、第２の移動局（利用者）１５及び基準局（固定局）１７はコード位相測位用の受信装置である。
【０００８】
図５を参照して従来の重力測定装置の構成及び動作を説明する。重力測定装置２０は航行体１０の基台１０Ａに装着されたジャイロコンパス２１と重力計２２とを有する。ジャイロコンパス２１は鉛直ジャイロであってよい。重力計２２は基台１０Ａに装着された水平安定台２３上に装着されている。水平安定台２３は重力計２２の入力軸を鉛直方向に保持するように構成され、その構造は後に説明する。重力測定装置２０は更に移動局１４、１５からの信号を入力する記録演算装置３１、３２を有する。
【０００９】
図６を参照して従来の水平安定台２３の構造を説明する。水平安定台は、プラットフォーム１０１とプラットフォーム１０１を水平に保持するためのジンバル装置を有する。プラットフォーム１０１上には重力計２２が配置されている。ジンバル装置は、プラットフォーム１０１の両端に装着されたロール軸１０３Ａ、１０３Ｂとロール軸の外端に設けられたロール軸トルカ１０５Ａ、１０５Ｂとロール軸１０３Ａ、１０３Ｂ及びロール軸トルカ１０５Ａ、１０５Ｂを支持するジンバル環１０７とジンバル環１０７の両端に装着されたピッチ軸１０９Ａ、１０９Ｂとピッチ軸の外端に設けられたピッチ軸トルカ１１１Ａ、１１１Ｂとピッチ軸１０９Ａ、１０９Ｂ及びピッチ軸トルカ１１１Ａ、１１１Ｂを支持する支持部材１１３Ａ、１１３Ｂとを有する。
【００１０】
ロール軸１０３Ａ、１０３Ｂとピッチ軸１０９Ａ、１０９Ｂは互いに直交するように装着される。ロール軸１０３Ａ、１０３Ｂとピッチ軸１０９Ａ、１０９Ｂは図示しない軸受けによって支持されている。支持部材１１３Ａ、１１３Ｂは航行体１０の基台１０Ａに装着される。水平安定台２３は、ロール軸１０３Ａ、１０３Ｂが航行体１０の首尾線に沿って配置されるように、航行体１０の基台１０Ａに装着される。
【００１１】
再び図５を参照して説明する。コード位相測地用の移動局１５からの速度信号はジャイロコンパス２１に供給される。ジャイロコンパス２１はこの速度信号を速度ログとして姿勢角（ロール角及びピッチ角）を正確に演算する。
【００１２】
ジャイロコンパス２１からの姿勢角（ロール角及びピッチ角）はロール軸トルカ１０５Ａ、１０５Ｂ及びピッチ軸トルカ１１１Ａ、１１１Ｂに供給される。それによってプラットフォーム１０１は水平に保持され、重力計２２の入力軸は鉛直方向に保持される。
【００１３】
搬送波位相測地用の移動局１４からは高精度の位置信号及び時刻信号が出力され、第１の記録演算装置３１に記録される。コード位相測地用の移動局１５からは時刻信号が出力され、第２の記録演算装置３２に記録される。重力計２２からは重力信号Ｇｔが出力され、リアルタイムにて第２の記録演算装置３２に記録される。重力計２２の出力は次の式によって表される。
【００１４】
【数１】
Ｇｔ＝ΔＧ＋Ａ＋Ｃ＋Ｚ＋Ｅ
【００１５】
Ｇｔは重力計２２の出力、ΔＧは重力異常、Ａは基準楕円体重力、Ｃはフリーエア補正、Ｚは鉛直加速度、Ｅはエトベス効果である。重力異常ΔＧは測定地の重力値と標準重力の偏差であり、１〜２００ｍＧａｌ程度である。重力測定装置２０はこの重力異常ΔＧを測定する。
【００１６】
基準楕円体重力Ａは測定点の緯度によって決まり、例えば緯度３５度では９７９．７Ｇａｌである。フリーエア補正Ｃは高度によって変化し、−０．３ｍＧａｌ／ｍである。鉛直加速度Ｚは航行体の上下運動による鉛直方向加速度である。これは航行体の位置を精密に測定し、その２階微分より求められる。その値は１〜１０Ｇａｌ程度である。エトベス効果Ｅは航行体の速度及び高度、基準楕円体の形状の影響によるコリオリ力、遠心力等によって変化し、０．１〜０．５Ｇａｌ程度である。
【００１７】
記録演算装置３１、３２は、先ず、入力信号に対して、航行体１０の振動等の雑音を除去するためのフィルタ処理をする。次に、数１の式の右辺の第２〜５項を演算し、それを重力計の出力より減算し、重力異常ΔＧを求める。記録演算装置３１、３２における演算はオフラインにてなされる。
【００１８】
【発明が解決しようとする課題】
従来の重力測定では、数１の式に示されているように、重力計に作用する鉛直方向の加速度Ｚのみが考慮されていた。これは重力計を支持するプラットフォームが水平であると仮定している。しかしながら、実際には、プラットフォームは必ずしも正確に水平に保持されているとは限らない。重力異常ΔＧの精度が１ｍＧａｌ以下の場合、プラットフォームの許容傾斜角Δθは次のようになる。
【００１９】
【数２】
Δθ＝ｃｏｓ^−１（１−１／９８００００）＝０．０８２［度］＝４．９［分］
【００２０】
プラットフォームが傾斜していると、重力計の出力には重力加速度の鉛直方向成分ばかりでなく水平加速度の鉛直方向成分が含まれる。この鉛直方向成分は誤差である。従って重力異常ΔＧを正確に測定することができない。
【００２１】
従って本発明は、航行体に搭載した重力測定装置によって重力を測定する場合、重力計を支持するプラットフォームが傾斜していても、正確に重力を測定することができることを目的とする。
【００２２】
本発明は、航行体に搭載した重力測定装置によって重力を測定する場合、航行体の加速度が大きい場合でも、正確に重力を測定することができることを目的とする。
【００２３】
【課題を解決するための手段】
本発明の重力測定装置は、重力計と、該重力計を支持するプラットフォームと、該プラットフォームを水平に支持するジンバル機構を有する水平安定台と、該水平安定台と同一面内に装着され上記水平安定台に姿勢角信号を供給するジャイロコンパスと、を有し、航行体に搭載される重力測定装置において、
上記プラットフォームには水平方向の加速度を検出するための加速度計が装着され、ＤＧＰＳによる搬送波位相測位用の移動局と該移動局からの位置情報及び時刻情報を記録するための第１の記録演算装置が設けられ、ＤＧＰＳによるコード位相測位用の移動局と該移動局からの位置情報及び時刻情報を記録するための第２の記録演算装置が設けられ、上記第２の記録演算装置は上記重力計の出力信号と上記加速度計の出力信号と上記ジャイロコンパスの出力信号を記録し、上記搬送波位相測位用の移動局からの位置信号によって東西方向加速度及び南北方向加速度を演算し、該東西方向加速度及び南北方向加速度と上記ジャイロコンパスからの方位角信号によってロール方向加速度及びピッチ方向加速度を演算し、上記加速度計からの水平加速度信号によって上記プラットフォームのロール角誤差及びピッチ角誤差をフィルタ演算し、それによって上記重力計の出力を補正するように構成されていることを特徴とする。
【００２４】
従って、航行体に搭載した重力測定装置によって重力を測定する場合、重力測定装置を支持するプラットフォームが傾斜しても、正確に重力異常を測定することができる。
【００２６】
従って、航行体に搭載した重力測定装置によって重力（重力異常）を測定する場合、オフライン処理によって重力異常を演算することができる。
【００２８】
本発明の重力測定方法は、重力測定装置を用いて、上記航行体の上記プラットフォームに装着された重力計により重力を測定する重力測定方法であって、ＤＧＰＳによる上記搬送波位相測位用の移動局からの位置情報によって東西方向加速度及び南北方向加速度を演算することと、該東西方向加速度及び南北方向加速度と上記ジャイロコンパスからの方位角信号によってロール方向加速度及びピッチ方向加速度を演算することと、上記加速度計からの水平加速度信号によって上記プラットフォームのロール角誤差及びピッチ角誤差を演算することと、上記ロール角誤差及びピッチ角誤差によって上記重力計の出力を補正することと、を含む。
【００２９】
【発明の実施の形態】
図１及び図２を参照して本発明の重力測定装置の例を説明する。本例の重力測定装置は航行体１０の基台１０Ａに装着されたジャイロコンパス２１と重力計２２と加速度計２４、２５とを有する。
【００３０】
図２に示すように、重力計２２及び加速度計２４、２５は基台１０Ａに装着された水平安定台２３のプラットフォーム１０１上に装着されている。重力測定装置２０は更に移動局１４、１５からの信号を入力する記録演算装置３１、３２を有する。記録演算装置３１、３２は図示のように別個に設けてもよいが、１個の記録演算装置であってよい。
【００３１】
本例の重力測定装置は図５に示した従来の重力測定装置と比較して、加速度計２４、２５が付加的に装着されている点が異なり、それ以外の構成は同様であってよい。
【００３２】
第１の加速度計２４は、ロール方向の加速度を検出し、第２の加速度計２５は、ピッチ方向の加速度を検出する。
【００３３】
図３を参照して本例の重力測定装置及び方法における信号の流れ及び処理を説明する。搬送波位相測地用の移動局１４によって航行体１０の東西加速度αＥ及び南北加速度αＮが検出され、第１の記録演算装置３１に記録される。ジャイロコンパス２１によって航行体１０の方位角φが検出され、第２の記録演算装置３２に記録される。これらの信号より、次の式によって、航行体のロール方向加速度α１及びピッチ方向加速度α２が求められる。
【００３４】
【数３】
α１＝αＥｃｏｓφ−αＮｓｉｎφ
α２＝αＥｓｉｎφ＋αＮｃｏｓφ
【００３５】
尚、図３の右下の絵に示すように、ロール方向加速度は航行体のピッチ軸線方向に沿った加速度であり、ピッチ方向加速度は航行体のロール軸線方向に沿った加速度である。
【００３６】
重力計２２の出力Ｇｔは第２の記録演算装置３２に記録される。ここで、重力計２２の出力は、数１の式の代わりに次の式によって表される。
【００３７】
【数４】

【００３８】
ここに、Ｈはプラットフォームの傾斜に起因した誤差である。それ以外の項ΔＧ、Ａ、Ｃ、Ｚ、Ｅは数１の式の項と同様である。Ｇ₀をここでは重力方向加速度と称する。
【００３９】
加速度計２４、２５の出力β１、β２は第２の記録演算装置３２に記録される。加速度計２４、２５の出力β１、β２は、次の式によって表される。
【００４０】
【数５】
β１＝α１’＋Ｇｔ・ｓｉｎθ１
β２＝α２’＋Ｇｔ・ｓｉｎθ２
【００４１】
数５の式の右辺の第１項α１’、α２’は航行体の運動に起因したプラットフォームのロール方向加速度及びピッチ方向加速度である。その大きさは５０Ｇａｌ程度であり、周期は２〜３分程度である。数５の式の右辺の第２項はプラットフォームの傾斜に起因した誤差である。その大きさは１Ｇａｌ程度であり、かなり長い周期を有する。
【００４２】
ロール角誤差及びピッチ誤差θ１、θ２は、それぞれプラットフォームのロール軸線周りの傾斜角及びピッチ軸線周りの傾斜角である。従って、プラットフォームが水平であるときには、この第２項は、ゼロとなる。以下に、加速度計２４、２５の出力β１、β２よりロール角誤差及びピッチ誤差θ１、θ２を求める手順を説明する。
【００４３】
数５の式の右辺の第１項と第２項では、その周期及びスペクトラムが異なるから適当なフィルタ処理によって両者を容易に分離することができる。従って、数５の式の右辺の第２項のみを取り出すことができる。
【００４４】
一方、数４の式の右辺において、プラットフォーム傾斜誤差Ｈは重力方向加速度Ｇ₀に比べて十分小さいとして無視する。更に、重力方向加速度Ｇ₀において、重力異常ΔＧ、フリーエア補正Ｃ、鉛直加速度Ｚ及びエトベス効果Ｅは、基準楕円体重力Ａに比べて十分小さいため、ここではそれらを無視する。従って次の式が成り立つ。
【００４５】
【数６】
Ｇｔ≒Ｇ₀≒Ａ
【００４６】
数５の式の右辺の第２項を、それぞれγ１、γ２とすると、次のように表される。
【００４７】
【数７】
γ１＝Ｇｔ・ｓｉｎθ１≒Ｇ₀・ｓｉｎθ１≒Ａ・ｓｉｎθ１
γ２＝Ｇｔ・ｓｉｎθ２≒Ｇ₀・ｓｉｎθ２≒Ａ・ｓｉｎθ２
【００４８】
基準楕円体重力Ａは予め求められているから、ロール角誤差及びピッチ誤差θ１、θ２が求まる。ロール角誤差及びピッチ誤差θ１、θ２が求められると、プラットフォーム傾斜誤差Ｈは次の式によって演算される。
【００４９】
【数８】
Ｈ＝α１ｓｉｎθ１＋α２ｓｉｎθ２
【００５０】
基準楕円体重力Ａ、フリーエア補正Ｃ、鉛直加速度Ｚ、エトベス効果Ｅは予め求められている。従って、プラットフォーム傾斜誤差Ｈが求められると、数４の式より、重力異常ΔＧが求められる。尚、数４の式の右辺の鉛直加速度Ｚは、搬送波位相測地用の移動局１４によって検出される航行体の位置情報及び時刻情報より、積分演算によって求められる。
【００５１】
図３に示すように、搬送波位相測地用の移動局１４からの東西加速度αＥ及び南北加速度αＮ、ジャイロコンパス２１からの方位角φ、加速度計２４、２５の出力β１、β２、重力計２２からの出力Ｇｔは、航行体１０上にて、測定点の位置情報及び時刻情報と共にリアルタイムにて、得られる。
【００５２】
航行体のロール方向加速度α１及びピッチ方向加速度α２、ロール角誤差θ１及びピッチ角誤差θ２、プラットフォーム傾斜誤差Ｈ及び重力異常ΔＧはオフラインにて地上にて求められる。重力計２２の出力Ｇｔに含まれる各種の変数及び定数Ａ、Ｃ、Ｚ、Ｅ等はオフラインにて地上にて求められる。
【００５３】
以上本発明の例について説明したが、本発明は上述の例に限定されることなく特許請求の範囲に記載された本発明の範囲にて様々な変更が可能であることは当業者にとって容易に理解されよう。
【００５４】
【発明の効果】
本発明によると、航行体に搭載した重力測定装置によって重力を測定する場合、重力測定装置を支持するプラットフォームが傾斜しても、正確に重力異常を測定することができる利点を有する。
【００５５】
本発明によると、航行体に搭載した重力測定装置によって重力を測定する場合、オフライン処理によって重力異常を演算することができる利点を有する。
【図面の簡単な説明】
【図１】本発明による重力測定装置の主要部を示す図である。
【図２】本発明による重力測定装置の水平安定台の構造を示す図である。
【図３】本発明による重力測定装置における演算の手順を示す図である。
【図４】従来の重力測定装置及び方法の概略を示す図である。
【図５】従来の重力測定装置の主要部を示す図である。
【図６】従来の重力測定装置の水平安定台の構造を示す図である。
【符号の説明】
１０・・・航行体、１０Ａ・・・基台、１２，１３・・・ＧＰＳ衛星、１４，１５・・・移動局、１６，１７・・・基準局、２０・・・重力測定装置、２１・・・ジャイロコンパス、２２・・・重力計、２３・・・水平安定台、２４，２５・・・加速度計、３１，３２・・・記録演算装置[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a gravity measuring apparatus and method for detecting gravity on the earth surface, and more particularly to a gravity measuring apparatus mounted on a navigation body and a gravity measuring method using the same.
[0002]
[Prior art]
The value of gravity on the earth's surface varies slightly depending on the measurement site, for reasons such as the density structure of underground rocks and changes in the geocentric distance. The difference between measured gravity and standard gravity is called gravity anomaly. The measurement of gravity is no different from the measurement of gravity anomaly. When measuring gravity, a portable gravity measurement device is used on flat ground, but in the polar regions, desert regions, mountainous regions, oceans, etc., gravity measurement devices mounted on navigational or moving objects such as vehicles and aircraft are used. The
[0003]
When the gravity is measured by the gravity measuring device mounted on the navigation body or the moving body, a DGPS (differential global positioning system) is used to accurately measure the position of the measurement point.
[0004]
DGPS was developed to improve the measurement accuracy of GPS. In GPS, the location is determined only by the mobile station (user), whereas in DGPS, the location is determined by the mobile station (user) and the reference station (fixed station). In DGPS, a correction value is calculated from a pseudo distance from a GPS satellite, time information, and orbit data at a reference station whose accurate position is measured in advance, and is transmitted to a mobile station. The mobile station corrects the signal from the GPS satellite using the correction value from the reference station.
[0005]
There are two methods for DGPS: code phase positioning and carrier phase positioning. The code phase positioning method uses the code phase as a measurement value and is used for real-time positioning. The accuracy is about 1 to 10 m. The carrier phase positioning method uses the carrier phase as a measurement value and is used for highly accurate position measurement. The accuracy is about mm to cm.
[0006]
An example of a conventional gravity measuring apparatus and method will be described with reference to FIG. The gravity measuring device 20 is mounted on a navigation body 10 such as a helicopter. The DGPS includes at least four GPS satellites (only two GPS satellites 12 and 13 are shown in FIG. 4), two mobile stations (users) 14 and 15 mounted on the navigation body 10, and two reference stations (fixed stations) on the ground. ) Uses 16 and 17.
[0007]
The first mobile station (user) 14 and the reference station (fixed station) 16 are receiving devices for carrier phase positioning, and the second mobile station (user) 15 and the reference station (fixed station) 17 are for code phase positioning. The receiving device.
[0008]
The configuration and operation of a conventional gravity measuring device will be described with reference to FIG. The gravity measuring device 20 includes a gyro compass 21 and a gravimeter 22 mounted on the base 10 </ b> A of the navigation body 10. The gyro compass 21 may be a vertical gyro. The gravimeter 22 is mounted on a horizontal stabilizer 23 mounted on the base 10A. The horizontal stabilizer 23 is configured to hold the input shaft of the gravimeter 22 in the vertical direction, and the structure will be described later. The gravity measuring device 20 further includes

recording operation devices

31 and 32 for inputting signals from the

mobile stations

14 and 15.
[0009]
The structure of the conventional horizontal stabilizer 23 will be described with reference to FIG. The horizontal stabilizer has a platform 101 and a gimbal device for holding the platform 101 horizontally. A gravity meter 22 is disposed on the platform 101. The gimbal device includes roll shafts 103A and 103B mounted on both ends of the platform 101,

roll shaft torquers

105A and 105B provided on the outer ends of the roll shafts, roll shafts 103A and 103B, and

roll shaft torquers

105A and 105B. Support for supporting pitch shafts 109A and 109B mounted on both ends of the ring 107 and the gimbal ring 107,

pitch shaft torquers

111A and 111B provided on the outer ends of the pitch shaft, pitch shafts 109A and 109B, and

pitch shaft torquers

111A and 111B.

Members

113A and 113B.
[0010]
The roll shafts 103A and 103B and the pitch shafts 109A and 109B are mounted so as to be orthogonal to each other. The roll shafts 103A and 103B and the pitch shafts 109A and 109B are supported by bearings (not shown). The

support members

113A and 113B are mounted on the base 10A of the navigation body 10. The horizontal stabilizer 23 is mounted on the base 10A of the navigation body 10 such that the roll shafts 103A and 103B are arranged along the tail line of the navigation body 10.
[0011]
A description will be given with reference to FIG. 5 again. The velocity signal from the code phase geodetic mobile station 15 is supplied to the gyrocompass 21. The gyrocompass 21 accurately calculates the posture angle (roll angle and pitch angle) using this speed signal as a speed log.
[0012]
The attitude angles (roll angle and pitch angle) from the gyrocompass 21 are supplied to the

roll axis torquers

105A and 105B and the

pitch axis torquers

111A and 111B. Thereby, the platform 101 is held horizontally, and the input shaft of the gravimeter 22 is held in the vertical direction.
[0013]
A highly accurate position signal and time signal are output from the mobile station 14 for carrier wave phase geodetic measurement and recorded in the first recording arithmetic unit 31. A time signal is output from the code phase geodetic mobile station 15 and recorded in the second recording arithmetic unit 32. A gravity signal Gt is output from the gravimeter 22 and recorded in the second recording arithmetic device 32 in real time. The output of the gravimeter 22 is expressed by the following equation.
[0014]
[Expression 1]
Gt = ΔG + A + C + Z + E
[0015]
Gt is the output of the gravimeter 22, ΔG is gravity anomaly, A is reference ellipsoid gravity, C is free air correction, Z is vertical acceleration, and E is the Etobes effect. Gravity abnormality (DELTA) G is a deviation of the gravity value of a measurement place, and standard gravity, and is about 1-200 mGal. The gravity measuring device 20 measures the gravity abnormality ΔG.
[0016]
The reference ellipsoid gravity A is determined by the latitude of the measurement point, and is 979.7 Gal at a latitude of 35 degrees, for example. Free air correction C varies with altitude and is -0.3 mGal / m. The vertical acceleration Z is the vertical acceleration due to the vertical movement of the navigation body. This is obtained by measuring the position of the navigation object precisely and calculating its second derivative. The value is about 1 to 10 Gal. The Etobes effect E varies depending on the speed and altitude of the navigation body, the Coriolis force due to the influence of the shape of the reference ellipsoid, the centrifugal force, etc., and is about 0.1 to 0.5 Gal.
[0017]
First, the

recording arithmetic devices

31 and 32 perform filter processing for removing noise such as vibration of the navigation body 10 on the input signal. Next, the second to fifth terms on the right side of Equation 1 are calculated and subtracted from the output of the gravimeter to obtain a gravity anomaly ΔG. The calculations in the

recording calculation devices

31 and 32 are performed off-line.
[0018]
[Problems to be solved by the invention]
In the conventional gravity measurement, as shown in the equation (1), only the vertical acceleration Z acting on the gravimeter is considered. This assumes that the platform supporting the gravimeter is horizontal. In practice, however, the platform is not necessarily held exactly level. When the accuracy of the gravity anomaly ΔG is 1 mGal or less, the allowable tilt angle Δθ of the platform is as follows.
[0019]
[Expression 2]
Δθ = cos ⁻¹ (1-1 / 9800000) = 0.082 [degrees] = 4.9 [minutes]
[0020]
When the platform is tilted, the output of the gravimeter includes not only the vertical component of gravitational acceleration but also the vertical component of horizontal acceleration. This vertical component is an error. Therefore, the gravity anomaly ΔG cannot be measured accurately.
[0021]
Accordingly, an object of the present invention is to accurately measure the gravity even when the platform supporting the gravimeter is inclined when the gravity is measured by the gravity measuring device mounted on the navigation body.
[0022]
An object of the present invention is to measure gravity accurately even when acceleration of a navigation body is large when measuring gravity with a gravity measurement device mounted on the navigation body.
[0023]
[Means for Solving the Problems]
The gravity measuring device of the present invention includes a gravimeter, a platform for supporting the gravimeter, a horizontal stabilizer having a gimbal mechanism for horizontally supporting the platform, and the horizontal stabilizer mounted on the same plane as the horizontal stabilizer. A gyrocompass for supplying a posture angle signal to a stable platform, and a gravity measuring device mounted on a navigation body,
The platform is equipped with an accelerometer for detecting acceleration in the horizontal direction, and a mobile station for carrier phase positioning by DGPS and a first recording operation for recording position information and time information from the mobile station A mobile station for code phase positioning by DGPS, and a second recording operation device for recording position information and time information from the mobile station, and the second recording operation device The output signal of the meter, the output signal of the accelerometer, and the output signal of the gyrocompass are recorded, and the east-west acceleration and the north-south acceleration are calculated from the position signal from the mobile station for carrier phase positioning, The roll direction acceleration and pitch direction acceleration are calculated from the north-south direction acceleration and the azimuth angle signal from the gyrocompass, and the horizontal direction from the accelerometer is calculated. The roll angle error and pitch angle errors of the platform filter operation by the speed signal, thereby characterized in that it is configured to correct the output of the gravimeter.
[0024]
Therefore, when the gravity is measured by the gravity measuring device mounted on the navigation body, the gravity abnormality can be accurately measured even when the platform supporting the gravity measuring device is inclined.
[0026]
Therefore, when measuring gravity (gravity abnormality) with the gravity measurement device mounted on the navigation body, the gravity abnormality can be calculated by offline processing.
[0028]
Gravity measurement method of the present invention, by using the gravimetric device, a gravity measurement method for measuring gravity by gravity meter mounted on the platform of the navigation body, from the mobile station for the carrier phase positioning by DGPS and computing the east-west direction acceleration and the north-south direction acceleration by the location information of the computing the roll direction acceleration and the pitch direction acceleration azimuthally signal from該東west direction acceleration and the north-south direction acceleration and the gyrocompass, the acceleration Calculating a roll angle error and a pitch angle error of the platform based on a horizontal acceleration signal from the meter, and correcting an output of the gravimeter based on the roll angle error and the pitch angle error.
[0029]
DETAILED DESCRIPTION OF THE INVENTION
An example of the gravity measuring device of the present invention will be described with reference to FIGS. The gravity measuring device of this example includes a gyrocompass 21, a gravimeter 22, and

accelerometers

24 and 25 that are mounted on a base 10 </ b> A of the navigation body 10.
[0030]
As shown in FIG. 2, the gravimeter 22 and the

accelerometers

24 and 25 are mounted on the platform 101 of the horizontal stabilization table 23 mounted on the base 10A. The gravity measuring device 20 further includes

recording operation devices

31 and 32 for inputting signals from the

mobile stations

14 and 15. The

recording operation devices

31 and 32 may be provided separately as shown in the figure, but may be a single recording operation device.
[0031]
The gravity measuring device of this example is different from the conventional gravity measuring device shown in FIG. 5 in that

accelerometers

24 and 25 are additionally attached, and the other configurations may be the same.
[0032]
The first accelerometer 24 detects the acceleration in the roll direction, and the second accelerometer 25 detects the acceleration in the pitch direction.
[0033]
With reference to FIG. 3, the flow and processing of signals in the gravity measuring apparatus and method of this example will be described. The mobile station 14 for detecting the carrier phase detects the east-west acceleration αE and the north-south acceleration αN of the navigation body 10 and records them in the first recording arithmetic unit 31. The azimuth angle φ of the navigation body 10 is detected by the gyrocompass 21 and recorded in the second recording arithmetic device 32. From these signals, the roll direction acceleration α1 and the pitch direction acceleration α2 of the navigation body are obtained by the following equations.
[0034]
[Equation 3]
α1 = αEcosφ-αNsinφ
α2 = αEsinφ + αNcosφ
[0035]
3, the roll direction acceleration is an acceleration along the pitch axis direction of the navigation body, and the pitch direction acceleration is an acceleration along the roll axis direction of the navigation body.
[0036]
The output Gt of the gravimeter 22 is recorded in the second recording calculation device 32. Here, the output of the gravimeter 22 is expressed by the following equation instead of the equation (1).
[0037]
[Expression 4]

[0038]
Here, H is an error due to the tilt of the platform. The other terms ΔG, A, C, Z, and E are the same as the terms in the equation (1). Here, G ₀ is referred to as gravity direction acceleration.
[0039]
The outputs β1 and β2 of the

accelerometers

24 and 25 are recorded in the second recording arithmetic device 32. The outputs β1 and β2 of the

accelerometers

24 and 25 are expressed by the following equations.
[0040]
[Equation 5]
β1 = α1 ′ + Gt · sin θ1
β2 = α2 ′ + Gt · sin θ2
[0041]
The first terms α1 ′ and α2 ′ on the right side of Equation 5 are the roll direction acceleration and the pitch direction acceleration of the platform due to the movement of the navigation body. Its size is about 50 Gal, and the period is about 2 to 3 minutes. The second term on the right side of Equation 5 is an error due to the tilt of the platform. Its size is about 1 Gal and has a fairly long period.
[0042]
The roll angle error and the pitch errors θ1 and θ2 are an inclination angle around the roll axis of the platform and an inclination angle around the pitch axis, respectively. Thus, this second term is zero when the platform is horizontal. The procedure for obtaining the roll angle error and the pitch errors θ1 and θ2 from the outputs β1 and β2 of the

accelerometers

24 and 25 will be described below.
[0043]
Since the period and spectrum of the first and second terms on the right side of Equation 5 are different, they can be easily separated by appropriate filter processing. Therefore, only the second term on the right side of Equation 5 can be extracted.
[0044]
On the other hand, on the right side of the equation (4), the platform tilt error H is ignored because it is sufficiently smaller than the gravity direction acceleration G ₀ . Further, since the gravity anomaly ΔG, the free air correction C, the vertical acceleration Z, and the Etobes effect E are sufficiently smaller than the reference ellipsoid gravity A at the gravity direction acceleration G ₀ , they are ignored here. Therefore, the following equation holds.
[0045]
[Formula 6]
Gt ≒ G ₀ ≒ A
[0046]
Assuming that the second term on the right side of Equation 5 is γ1 and γ2, respectively, it is expressed as follows.
[0047]
[Expression 7]
γ1 = Gt · sin θ1≈G ₀ · sin θ1≈A · sin θ1
γ2 = Gt · sin θ2≈G ₀ · sin θ2≈A · sin θ2
[0048]
Since the reference ellipsoid gravity A is obtained in advance, roll angle errors and pitch errors θ1 and θ2 are obtained. When the roll angle error and the pitch errors θ1 and θ2 are obtained, the platform tilt error H is calculated by the following equation.
[0049]
[Equation 8]
H = α1sin θ1 + α2sin θ2
[0050]
The reference ellipsoid gravity A, free air correction C, vertical acceleration Z, and Etobes effect E are obtained in advance. Therefore, when the platform tilt error H is obtained, the gravity anomaly ΔG is obtained from the equation (4). Note that the vertical acceleration Z on the right side of the equation (4) is obtained by integral calculation from the position information and time information of the navigation body detected by the mobile station 14 for carrier wave phase geodetic measurement.
[0051]
As shown in FIG. 3, the east-west acceleration αE and the north-south acceleration αN from the mobile station 14 for carrier phase measurement, the azimuth angle φ from the gyrocompass 21, the outputs β1 and β2 of the

accelerometers

24 and 25, and from the gravimeter 22 The output Gt is obtained on the navigation body 10 in real time together with the position information and time information of the measurement points.
[0052]
The roll direction acceleration α1 and pitch direction acceleration α2, the roll angle error θ1 and the pitch angle error θ2, the platform tilt error H, and the gravity anomaly ΔG of the navigation object are obtained on the ground offline. Various variables and constants A, C, Z, E, and the like included in the output Gt of the gravimeter 22 are obtained on the ground offline.
[0053]
Although the examples of the present invention have been described above, the present invention is not limited to the above-described examples, and various modifications can be easily made by those skilled in the art within the scope of the present invention described in the claims. It will be understood.
[0054]
【The invention's effect】
According to the present invention, when the gravity is measured by the gravity measuring device mounted on the navigation body, there is an advantage that the gravity abnormality can be accurately measured even if the platform supporting the gravity measuring device is inclined.
[0055]
According to the present invention, when gravity is measured by the gravity measuring device mounted on the navigation body, there is an advantage that gravity abnormality can be calculated by offline processing.
[Brief description of the drawings]
FIG. 1 is a diagram showing a main part of a gravity measuring apparatus according to the present invention.
FIG. 2 is a diagram showing a structure of a horizontal stabilizer of the gravity measuring apparatus according to the present invention.
FIG. 3 is a diagram showing a calculation procedure in the gravity measuring apparatus according to the present invention.
FIG. 4 is a diagram schematically showing a conventional gravity measuring apparatus and method.
FIG. 5 is a diagram showing a main part of a conventional gravity measuring device.
FIG. 6 is a diagram showing a structure of a horizontal stabilization base of a conventional gravity measuring device.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 10 ... Navigation body, 10A ... Base, 12, 13 ... GPS satellite, 14, 15 ... Mobile station, 16, 17 ... Reference station, 20 ... Gravity measuring device, 21. .. Gyrocompass, 22 ... Gravimeter, 23 ... Horizontal stabilizer, 24, 25 ... Accelerometer, 31, 32 ... Recording operation device

Claims

A gravimeter, a platform that supports the gravimeter, a horizontal stabilizer that has a gimbal mechanism that supports the platform horizontally, and a posture angle signal that is mounted in the same plane as the horizontal stabilizer and supplies an attitude angle signal to the horizontal stabilizer A gravity measuring device mounted on the navigation body,
The platform is equipped with an accelerometer to detect horizontal acceleration ,
A mobile station for carrier phase positioning by DGPS and a first recording operation device for recording position information and time information from the mobile station are provided,
A mobile station for code phase positioning by DGPS and a second recording operation device for recording position information and time information from the mobile station are provided,
The second recording arithmetic unit records the output signal of the gravimeter, the output signal of the accelerometer, and the output signal of the gyrocompass,
The east / west direction acceleration and the north / south direction acceleration are calculated from the position signal from the carrier phase positioning mobile station, and the roll direction acceleration and the pitch direction acceleration are calculated by the east / west direction acceleration and the north / south direction acceleration and the azimuth angle signal from the gyrocompass. Gravity measurement characterized by being configured to calculate and filter the roll angle error and pitch angle error of the platform by the horizontal acceleration signal from the accelerometer, thereby correcting the output of the gravimeter apparatus.

A gravity measurement method using the gravity measurement device according to claim 1 to measure gravity by a gravimeter attached to the platform of the navigation body ,
And computing the east-west direction acceleration and the north-south direction acceleration by the location information from the mobile station for the carrier phase positioning by DGPS,
And computing the roll direction acceleration and the pitch direction acceleration azimuthally signal from該東west direction acceleration and the north-south direction acceleration and the gyrocompass,
And computing the roll angle error and pitch angle errors of the platform by the horizontal acceleration signal from the accelerometer,
Correcting the output of the gravimeter by the roll angle error and the pitch angle error;
Gravity measurement method including.