光學(xué)干涉式重力儀的數(shù)據(jù)處理方法

黃時康，黃騰超

（浙江大學(xué) 光電科學(xué)與工程學(xué)院，浙江 杭州 310013）

光學(xué)干涉式重力儀的數(shù)據(jù)處理方法

黃時康，黃騰超

（浙江大學(xué) 光電科學(xué)與工程學(xué)院，浙江 杭州 310013）

本文研究了光學(xué)干涉式重力儀的關(guān)鍵技術(shù)中的數(shù)據(jù)處理部分，闡述了對重力儀測出條紋進行分析處理的兩種算法：過零探測法及二階差分法的原理.并在Matlab平臺下針對這兩種算法進行了編程仿真，對用兩種算法處理同一組數(shù)據(jù)的仿真結(jié)果進行了比較.

光學(xué)；干涉；重力儀；過零探測法；二階差分法

重力場是反映地球內(nèi)部物質(zhì)結(jié)構(gòu)及其變遷的地球物理基本場，其主要參數(shù)為重力加速度(g，常用值9.81m/s2).g值的精密測量是探測地球重力場的重要途徑.重力值測量是計量、測繪、地質(zhì)、地震與資源勘探等領(lǐng)域不可或缺的部分，其在大地測量、地球物理、精密計量和軍用導(dǎo)航等領(lǐng)域具有十分重要的意義和應(yīng)用價值.

本文的研究內(nèi)容是第五項關(guān)鍵技術(shù)數(shù)據(jù)處理與儀器控制系統(tǒng)中的數(shù)據(jù)處理以及誤差分析部分.我們根據(jù)兩種重力值的算法編寫程序，利用MATLAB軟件對兩種算法進行了仿真，對重力儀測量的數(shù)據(jù)進行處理，計算出了重力值g并進行了比較.

1 數(shù)據(jù)處理方法

1.1 過零探測法

落體下落距離S（t）的數(shù)學(xué)模型為：

其中，s0為物體的初始位置，v0為初速度.γ為垂直重力梯度.若重力梯度未知，由于其數(shù)值很?。?0-6數(shù)量級）則可以忽略不計.

過零探測法的工作原理如下：當物體下落時，重力儀的采點頻率為fs=100MHz，如圖1所示：

圖 1過零探測法原理圖

當?shù)趇個采樣點Ui和下一個采樣點Ui+1相乘小于0 時，在Ui與Ui+1之間必定存在一個零點.在Matlab平臺下，可以在數(shù)秒內(nèi)完成占存儲空間幾百Mbite數(shù)據(jù)量的過零探測.零點的準確時刻tk可由相似三角形關(guān)系算出：tk=ti+Ui·Ts/（Ui-Ui+1）.其中，Ts為相鄰兩點的間隔時間，其值為1/fs，同理，ti=（i-1）/fs.由于兩個相鄰的零點之間為1/2個條紋周期，其相位變化為π，故物體下落距離為λ/4.所以，在tk時刻物體的位置為Sk(tk)=(k-1)·(λ/4).將Sk(tk)與tk的值代入公式（1.1）中即可求出重力值g

2.2 二階差分法

二階差分法的核心思想為，當忽略重力梯度γ的影響時，取固定的Δt，下落距離S（t）的二階差分Δ2(S(t)=g(Δt)2為常數(shù).

地球的重力值在9.78ms-2至9.83ms-2之間，取 Δt= 0.2ms，則 Δ2在 3.912×10-7m 與 3.932×10-7m之間.測量光源為He-Ne激光器，波長λ=632.8nm.則條紋數(shù)差值為Δ2/（λ/2），其數(shù)值在1.236與1.242之間.也就是說，將條紋以Δt=0.2ms分為N段，每一段的條紋數(shù)Kn會比前一段的條紋數(shù)Kn-1多1.236至1.242個.

也就是說，Kn-Kn-1的整數(shù)部分是穩(wěn)定不變的C=1，我們要確定的是小數(shù)部分.

圖 2二階差分法原理圖

要確定小數(shù)部分，只需在每段Δt部分的同一位置取一個點作為特征點，用過零探測的方法求出與該點最近的兩個零點，根據(jù)該點與前后零點之間點的個數(shù)，來求出該點在一個周期條紋中的位置即Φ(n).舉個例子，若第n段的特征點是第100個點，距離最近的兩個零點是第90個點與第140個點，則Φ(n)=10/50=0.2.將每段的Φ(n)求出后，Φ(n)-Φ(n-1)即為條紋數(shù)差Kn-Kn-1的小數(shù)部分.

在求出Kn-Kn-1之后，下落距離，之后代入(1.1)中即可計算出重力值g.

2 實驗結(jié)果分析

下表為通過兩種算法計算的重力值（重復(fù)實驗20次，共20組數(shù)據(jù)）

表1 實驗數(shù)據(jù)記錄表

對兩種方法處理得到的重力值用Matlab畫成平滑曲線圖.

圖3 兩種算法計算出的重力值比較示意圖

其中，紅線代表過零探測法，藍線代表二階差分法.從圖中可看出，兩種算法得出的結(jié)果非常接近，幾乎沒有區(qū)別.

過零探測法：

方差：D（g）=8.9428×10-7ms-2

標準差：σ=9.4566×10-4ms-2

二階差分法：

方差：D（g）=9.0056×10-7ms-2

標準差：σ=9.4898×10-4ms-2

3 結(jié)論

本文提出了兩種提取重力儀測得的干涉條紋上的信息，計算重力值g的算法（過零探測法與二階差分法），對誤差來源進行了原理分析，并進行了數(shù)學(xué)建模，計算出了誤差的大小.并在Matlab平臺下，編寫程序，對兩種算法進行了仿真計算，并將計算結(jié)果進行了比較.結(jié)果表明，兩種算法的結(jié)果幾乎一致，但計算結(jié)果與真實重力值的差值還待進一步的研究.

〔1〕吳瓊.高精度絕對重力儀關(guān)鍵技術(shù)研究[D].中國地震局地球物理研究所,2011.

〔2〕Dr.Eberhard Bansch，Lijun Wang,Ph.D，Dr.Thomas Fauster.Ultra-high Precision,Absolute，Earth Gravity Measurements，26. Juni 2008.

〔3〕S Svitlov1,P Mas?yk,Ch Rothleitner,H Hu and L J Wang.Comparison of three digital fringe signal processing methods in a ballistic free-fall absolute gravimeter，11 November 2010.

〔4〕高景龍.NIM-3型新的輕小高精度可移式絕對重力儀[J].測繪學(xué)報，1993，22（3）.

〔5〕郭有光，黃大倫，方永源，徐進義，劉達倫.中國NIM型絕對重力儀及國際絕對重力儀比對[J].現(xiàn)代計量測試,1998（6）.

〔6〕W.TORGE,R.H.RDER,M.SCHNIJ LL,H.-G.WENZEL.FIRSTRESULTSWITH THE TRANSPORTABLE ABSOLUTE GRAVITY METER JILAG-3,1987.

〔7〕T M.Niebauel;G.S.Sasagawa,J.E.Fallel;R.Hilt and E Klopping，A new generation of absolute gravimeters，1995.

〔8〕J.Liard and C.Gagnon.The new A-10 absolute gravimeter at the 2001International Comparison of Absolute Gravimeters.2002.

〔9〕Giancarlo D’Agostino,Sergio Desogus,A-lessandro Germak, Claudio Origlia, Danilo Quagliotti,Giovanna Berrino,Gennaro Corrado,Vincenzo d’Errico and Giuseppe Ricciardi,The new IMGC-02 transportable absolute gravimeter:measurement apparatus and applications in geophysics and volcanology,2008.

〔10〕Ch Rothleitner1,S Svitlov,H M′erim`eche2,H Hu and L J Wang,Development of new free-fall absolute gravimeters,2009.

〔11〕HHu,S Svitlov,C Rothleitner,J Schafer,J Zhang and L J Wang,Improvements of the MPG -2 transportable absolute ballistic gravimeter,2010.

〔12〕Giancarlo D’Agostino,Sergio Desogus,A-lessandro Germak,Claudio Origlia,Danilo Quagliotti,Giovanna Berrino,Gennaro Corrado,Vincenzo d’Errico and Giuseppe Ricciardi,The new IMGC-02 transportable absolutegravimeter:measurementapparatusand applicationsin geophysicsand volcanology,2008.

〔13〕J.E.Faller,Thirty years of progress in absolute gravimetry:a scientific capability implemented bytechnological advances,2002.

TB96

A

1673-260X（2017）09-0036-03

2017-07-02