捷聯(lián)式航空重力矢量測量低頻誤差機(jī)理和修正

蔡劭琨，曹聚亮，于瑞航

（國防科技大學(xué) 智能科學(xué)學(xué)院，長沙 410072）

重力信息對于國民經(jīng)濟(jì)建設(shè)和軍事應(yīng)用都具有重要的意義。地面單點(diǎn)重力測量是獲取地球重力場高頻信息的傳統(tǒng)方法，但是該方法在人跡罕至的地區(qū)是難以實(shí)施的，并且該方法效率很低，花費(fèi)較高；海洋重力測量同樣也不是一種效率很高的重力測量方法，因?yàn)槠浜剿佥^慢，并且只能在海面上進(jìn)行作業(yè)；衛(wèi)星重力測量可以方便地對全球重力信息進(jìn)行測量，但是只能得到地球重力場的長波部分，重力場的高頻信息不能通過衛(wèi)星數(shù)據(jù)獲得。相比這幾種重力測量方式，航空重力測量是快速獲取大面積地球重力場中高頻信息的一種有效手段。所謂的航空重力測量就是將重力傳感器、姿態(tài)傳感器、GNSS 等安裝在飛機(jī)上進(jìn)行重力測量的一種方式，這種測量方式用時少、效費(fèi)比高，是目前重力測量的一種主要手段[1]。

航空重力測量儀器的基本結(jié)構(gòu)可以分成三類：雙軸穩(wěn)定平臺結(jié)構(gòu)、三軸平臺式慣性導(dǎo)航系統(tǒng)結(jié)構(gòu)和捷聯(lián)式慣性導(dǎo)航系統(tǒng)結(jié)構(gòu)。

基于雙軸穩(wěn)定平臺的重力儀是將重力傳感器安裝在一個穩(wěn)定平臺上，通過穩(wěn)定平臺來保持重力傳感器的指向。典型的基于雙軸穩(wěn)定平臺的重力儀有LCR海/空重力儀、Chekan 重力儀。該種類型的重力儀可以達(dá)到1 mGal/5 km（1 mGal=10-5m/s2）的精度水平[2]。

基于三軸平臺慣導(dǎo)的重力儀其核心是一套三軸平臺慣導(dǎo)系統(tǒng)，通過三軸平臺慣導(dǎo)系統(tǒng)來將三軸加速度計(jì)組合保持在當(dāng)?shù)氐乩碜鴺?biāo)系中。典型的有AIRGrav重力儀和GT 系列重力儀。這種類型重力儀可以達(dá)到1 mGal/3 km 的精度[3,4]。

“業(yè)無高卑志當(dāng)堅(jiān)男兒有求安得閑，工程建設(shè)雖然苦累，但是我依舊選擇了這一行?！币磺粺嵫借F建參與祖國建設(shè)，這也許就是王維振的夢想。

將式(9)代入到式(8)中，并將慣性元件誤差簡化為只有加表零偏（bax、b ay、b az）和陀螺漂移（bgx、b gy、bgz）。假設(shè)載體坐標(biāo)系為前右下，x、y、z三個加表和陀螺分別沿載體坐標(biāo)系的前向、右向和垂向安裝，導(dǎo)航坐標(biāo)系為北東地坐標(biāo)系，測量載機(jī)由東向西飛行，由此可得重力儀原始測量誤差與慣性元件誤差之間的關(guān)系：

雖然現(xiàn)金流動性較強(qiáng)，但其盈利能力并不高。Jenson（1986）提出的自由現(xiàn)金流量假說認(rèn)為，如果公司持有過多不必要的現(xiàn)金，管理者會不知不覺地實(shí)施現(xiàn)金濫用，從而盲目地進(jìn)行投資項(xiàng)目，這在無形之間增加了企業(yè)經(jīng)營期間的財務(wù)風(fēng)險，也是對股東權(quán)益損害較大的一種自私自利的行為。由此可知，董事會需要抑制這種惡性投資行為來提高現(xiàn)金配置效率。但是當(dāng)前現(xiàn)狀是不論出于什么樣的原因，很多公司都愿意超額持有現(xiàn)金量以滿足自己的財務(wù)報表的需求。

一類誤差是由慣性導(dǎo)航系統(tǒng)帶來的，包括由姿態(tài)誤差ψ產(chǎn)生的誤差項(xiàng)[f n×]ψ以及由加速度計(jì)測量誤差產(chǎn)生的誤差項(xiàng)C bnδfb；

在可以進(jìn)行重力矢量測量的兩類重力儀中，基于捷聯(lián)慣導(dǎo)系統(tǒng)的重力儀相比于基于三軸平臺慣導(dǎo)系統(tǒng)的重力儀，其結(jié)構(gòu)較為簡單、體積和重量都比較小、安裝方便、操作簡易，是重力儀的一個重要發(fā)展方向。本文基于國防科技大學(xué)研制的SGA 系列捷聯(lián)式重力儀，介紹了捷聯(lián)式重力矢量測量的原理，并分析重力矢量測量低頻誤差的產(chǎn)生機(jī)理及相應(yīng)的修正方法。

本文在借鑒CDP項(xiàng)目分類標(biāo)準(zhǔn)的基礎(chǔ)上，參考環(huán)保部公布的 《上市公司環(huán)境信息披露指南》，并結(jié)合樣本公司碳信息披露的實(shí)際情況，建立了涵蓋戰(zhàn)略規(guī)劃、治理架構(gòu)、風(fēng)險或機(jī)遇識別、碳排放核算四大類的上市公司碳信息披露評價體系。在碳信息披露評價體系的框架下，運(yùn)用Python軟件實(shí)現(xiàn)對樣本公司所披露文件的內(nèi)容爬取與文本分析，最后通過算法匯總企業(yè)碳信息披露所獲得分。碳信息披露評價體系見表1。

1 重力矢量測量數(shù)學(xué)模型及誤差分析

根據(jù)牛頓第二定律，質(zhì)點(diǎn)的運(yùn)動方程在慣性坐標(biāo)系i中可以表示為：

通過FMA評分判斷，治療前兩組患者無明顯差異，治療后兩組患者相比，差異有統(tǒng)計(jì)學(xué)意義（P＜0.05）。見表1。

式中，Gi表示引力加速度；表示質(zhì)點(diǎn)在慣性坐標(biāo)系i中的運(yùn)動加速度；xi表示質(zhì)點(diǎn)在慣性坐標(biāo)系i中的位置矢量；fi表示慣性坐標(biāo)系i中的比力，由重力傳感器獲得。

式(1)表明，不能從重力傳感器測得的比力fi中將引力加速度和運(yùn)動加速度區(qū)分開來。對于靜態(tài)重力測量，此時質(zhì)點(diǎn)運(yùn)動加速度等于零，重力傳感器輸出值在大小上等于引力加速度。但是在動態(tài)的重力測量中質(zhì)點(diǎn)的運(yùn)動加速度非零，此時不能從重力傳感器的測量結(jié)果中將運(yùn)動加速度扣除，因此，動態(tài)條件下僅僅通過單一的重力傳感器并不能測量得到重力信息。

動態(tài)重力測量的基本原理可通過對式(1)進(jìn)行變換得到：

式(2)可以表示成在當(dāng)?shù)氐乩碜鴺?biāo)系n中的形式[9]：

式中，γn表示正常重力矢量在n系中的投影；v en表示載體的速度；表示載體的加速度；fb表示載體坐標(biāo)系b下的比力信息；Cbn表示導(dǎo)航坐標(biāo)系n（即當(dāng)?shù)氐乩碜鴺?biāo)系）和載體坐標(biāo)系b之間的轉(zhuǎn)換矩陣；ωenn表示n系相對于地球的轉(zhuǎn)動角速度在n系中的表示；ω ine表示地球相對于慣性系i的轉(zhuǎn)動角速度在n系中的投影；δgn表示重力矢量測量所獲得的擾動矢量信息。式(3)的右邊部分根據(jù)來源不同可以分為兩類：

步驟5：在步驟4 得到的重力擾動中存在大量的高頻噪聲，噪聲的幅值是重力信號的幾百甚至幾千倍，重力信號比噪聲大的頻帶只集中在低頻段的很窄范圍內(nèi)，需要通過低通濾波器來提取重力擾動信息（示意圖如圖5所示）。并且由于重力信號和噪聲之間沒有明顯的過渡，低通濾波器的截止頻率較大（截止頻率大，對應(yīng)的波長則較短）雖然有利于提高測量結(jié)果的空間分辨率，但會使得殘留噪聲增多，從而導(dǎo)致精度降低；相反，如果低通濾波器的截止頻率較?。ㄏ鄳?yīng)的波長較長），雖然可以提高濾波的效果，同時也會有更多的重力信號被一起濾除，從而影響測量結(jié)果的空間分辨率。因此，需要根據(jù)任務(wù)的需求綜合考慮測量精度和空間分辨率，選取合適的低通濾波截止頻率；

從圖2可看出，谷索的就位對吊索索力影響很??；脊索就位對不同索的索力影響不同；膜的安裝將提高吊索的索力。

式中，δgD表示重力擾動的垂向分量，即標(biāo)量測量的結(jié)果；δgN表示重力擾動的北向分量；δgE表示重力擾動的東向分量；vN、vE和vD分別表示載體運(yùn)動速度的北向、東向和垂向分量；和分別表示載體運(yùn)動加速度的北向、東向和垂向分量；ω ie表示地球自轉(zhuǎn)角速度的大??；RM和RN分別表示子午圈曲率半徑和卯酉圈曲率半徑；L和h分別表示緯度和高度；γ表示正常重力的數(shù)值大小。

改革開放前我們重視公平而忽視效率，把共同富裕當(dāng)作同步富裕、同時富裕和同等富裕去追求，導(dǎo)致社會的同步貧窮、同等貧窮，造成分配上的平均主義和“大鍋飯”。改革開放后我們追求公平與效率的統(tǒng)一，允許一部分人、一部分地區(qū)先富起來，打破了平均主義和“大鍋飯”，調(diào)動勞動者積極性，解放和發(fā)展生產(chǎn)力，引起經(jīng)濟(jì)、社會和精神生活的深刻變化。所以，先富起來的辦法比平均主義的老辦法好，有利于實(shí)現(xiàn)先富。

對式(3)進(jìn)行變分，可以得到捷聯(lián)式航空重力矢量測量的誤差模型：

同樣，式(5)右邊的誤差項(xiàng)可以分成兩類：

在這三種類型的重力儀中，基于雙軸穩(wěn)定平臺的重力儀通常只配置了一個重力傳感器，因此只能進(jìn)行重力標(biāo)量測量，基于捷聯(lián)式慣性導(dǎo)航系統(tǒng)結(jié)構(gòu)和三軸平臺式慣性導(dǎo)航系統(tǒng)結(jié)構(gòu)的重力儀不僅可以對重力的標(biāo)量信息進(jìn)行測量，還可以實(shí)現(xiàn)對重力矢量信息的測量（也就是對垂線偏差進(jìn)行測量）。目前，進(jìn)行重力標(biāo)量測量的重力儀已經(jīng)有成熟的商業(yè)產(chǎn)品（LCR、Chekan、AIRGrav、GT），重力矢量測量儀還沒有成熟的產(chǎn)品，大部分處于實(shí)驗(yàn)室階段。但是與標(biāo)量測量相比，矢量測量在大地測量、固體地球物理、海洋學(xué)和資源勘探等方面具有優(yōu)勢[7]。首先，雖然垂線偏差可以利用費(fèi)寧-梅內(nèi)斯（Vening-Meinesz）積分通過重力異常信息間接得到，但是這種方法理論上需要全球覆蓋的重力異常信息，這一條件是很難滿足的，尤其是在極點(diǎn)區(qū)域。如果可以直接對垂線偏差進(jìn)行測量的話，那么就不需要這么嚴(yán)格的條件，因此，重力矢量測量可以快遞、經(jīng)濟(jì)地獲得更精確的結(jié)果；其次，大地水準(zhǔn)面對于固體地球物理中的均衡補(bǔ)償具有很重要的意義，而大地水準(zhǔn)面剖面可以通過對測線上的垂線偏差進(jìn)行積分得到，并且沿測線對垂線偏差進(jìn)行積分的方法不存在邊界效應(yīng)的問題[7]；再次，借助于由垂線偏差得到的大地水準(zhǔn)面資料，冰川學(xué)家可以研究冰川移動中的重力不平衡因素；第四，在海洋學(xué)中，結(jié)合海洋大地水準(zhǔn)面信息和衛(wèi)星測高確定的海洋表面，可以計(jì)算出所關(guān)心的洋流信息，而其中的海洋大地水準(zhǔn)面可以通過垂線偏差來計(jì)算得到[8]。

內(nèi)分泌失調(diào)是臨床女性常見生理現(xiàn)象,由環(huán)境、生理及營養(yǎng)等方面因素引起。近年來,隨著生活節(jié)奏不斷加快,工作及家庭對女性造成的壓力急劇上升,導(dǎo)致內(nèi)分泌失調(diào)發(fā)生率增加。內(nèi)分泌出現(xiàn)失調(diào)后,患者極易出現(xiàn)痤瘡、肥胖、不育及婦科、乳房等表現(xiàn)和疾病,對其正常生活及工作造成影響。因此,有必要立即對其實(shí)施對應(yīng)治療,確保其健康。本研究選取300例內(nèi)分泌失調(diào)患者作為研究對象,對其中150例采用克羅米芬聯(lián)合絨毛膜促性腺激素、人絕經(jīng)促性腺激素治療取得滿意效果,現(xiàn)報道如下。

式(3)可以寫成分量的形式：

其中，在當(dāng)前GNSS 精度條件下，由GNSS 引起的誤差已經(jīng)很小了，在這里不做討論。因此，本文重點(diǎn)討論由慣導(dǎo)系統(tǒng)引起的誤差。在這種情況下，式(5)表示的航空重力矢量測量誤差模型可以簡化成：

展開成分量形式如下：

式中，fN表示加速度計(jì)比力測量值的北向分量；fE表示加速度計(jì)比力測量值的東向分量；fD表示加速度計(jì)比力測量值的垂向分量；δfN表示加速度計(jì)比力測量誤差的北向分量；δfE表示加速度計(jì)比力測量誤差的東向分量；δfD表示加速度計(jì)比力測量誤差的垂向分量；ψ N、ψ E和ψ D分別表示沿當(dāng)?shù)氐乩碜鴺?biāo)系三個軸的姿態(tài)測量誤差。

在航空重力測量中，要求飛機(jī)沿著測線水平勻速直線飛行。此時有fD≈g?f N,fE，其中g(shù)為重力加速度大小，式(7)可以寫成：

The level of A-CoA of the model group was increased obviously compared with that of the blank group (P< 0.05), and the level of A-CoA of the FZ and LWDF groups were decreased compared with that of the model group (P< 0.01; Fig. 3).

雖然采用組合算法能夠在很大程度上估計(jì)慣導(dǎo)的姿態(tài)誤差并予以抑制，但是并不能將姿態(tài)誤差完全消除，這也正是引起重力矢量測量的一個重要的誤差源。所以舒勒周期會一直存在于測量結(jié)果中。姿態(tài)誤差ψ主要是由陀螺漂移bg引起，可以近似表示為：

基于捷聯(lián)慣性導(dǎo)航系統(tǒng)的重力儀其核心部分是一套捷聯(lián)式慣性導(dǎo)航系統(tǒng)。通過捷聯(lián)式慣導(dǎo)系統(tǒng)計(jì)算出載體的姿態(tài)，然后將三軸加速度計(jì)的輸出投影到當(dāng)?shù)氐乩碜鴺?biāo)系中。典型的系統(tǒng)有SISG[5]和國防科技大學(xué)研制的SGA 系列重力儀[6]。此種類型重力儀目前可以達(dá)到1-3 mGal 的精度，相應(yīng)的分辨率為1-3 km。

從式(8)可以看出，對于重力標(biāo)量測量而言，誤差主要來自于加速度計(jì)測量誤差，這一誤差通常表現(xiàn)為常值誤差或者線性漂移誤差，可以通過機(jī)場停機(jī)坪處的重力真值來進(jìn)行校正。而對于重力矢量測量，重力擾動水平分量測量結(jié)果不僅受到加速度計(jì)測量誤差的影響，更重要的是還受到水平姿態(tài)測量誤差ψN和ψE的影響，1 ?的姿態(tài)測量誤差會引起4.75 mGal 的水平分量測量誤差，而且姿態(tài)誤差并不表現(xiàn)為常值或者線性的形式，很難通過機(jī)場停機(jī)坪處的垂線偏差真值來進(jìn)行校正。

實(shí)際上，如圖1和圖2所示，在慣性導(dǎo)航系統(tǒng)中，水平姿態(tài)誤差存在舒勒周期（藍(lán)色曲線），并且這一周期性振蕩的幅值受到傅科周期的調(diào)制（紅色曲線）。需要指出的是，由于不知道水平姿態(tài)誤差的初始相位，因此不能通過圖1所示的姿態(tài)誤差曲線對水平姿態(tài)誤差進(jìn)行補(bǔ)償。但是在水平姿態(tài)誤差曲線中存在舒勒周期和傅科周期兩個周期，舒勒周期的周期為84.4 min（5064 s），傅科周期的周期在兩極等于24 h，緯度越低，周期越大。那么重力水平分量的測量誤差中也將存在這兩個周期。對于典型的小型航測飛機(jī)，飛行速度大約為60 m/s，舒勒周期對應(yīng)的波長為304 km，也就是說，捷聯(lián)式航空重力矢量測量水平分量的測量誤差主要影響波長大于304 km 的重力場長波部分，在重力水平分量測量結(jié)果中，波長大于304 km 的部分除了有真實(shí)的垂線偏差信息外，還包含姿態(tài)測量誤差，是一種低頻誤差[5]?，F(xiàn)在需要做的就是在低頻段內(nèi)對重力水平分量測量結(jié)果中的誤差進(jìn)行修正。

圖1 慣導(dǎo)系統(tǒng)水平姿態(tài)誤差之滾動角誤差Fig.1 Roll error in Horizontal attitude errors of the inertial navigation system

圖2 慣導(dǎo)系統(tǒng)水平姿態(tài)誤差之俯仰角誤差Fig.2 Pitch error in horizontal attitude errors of the inertial navigation system

如果將水平分量測量結(jié)果中波長大于304 km 的低頻部分濾除，那么得到的重力水平分量信息只是相對重力信息，要想將估計(jì)出來的相對重力場與絕對重力場聯(lián)系起來，就必須知道重力場的一些先驗(yàn)信息。而地球重力場模型正好可以提供這樣一種重力場的先驗(yàn)信息。

2 地球重力場模型

EGM2008（Earth Gravitational Model 2008）是由美國發(fā)布的全球超高階地球重力場模型[10]，該模型是用一組球諧函數(shù)的系數(shù)來表示地球重力場，模型的階次完全至2159 階，相當(dāng)于模型的半波長空間分辨率為5 ′（約9 km）。

可以利用EGM2008 模型通過下面兩個式子來計(jì)算地球上指定點(diǎn)的重力擾動水平分量：

式中，GM= 3.986004415 ×1014m3/s2表示牛頓重力常數(shù)與地球質(zhì)量的乘積；a= 6378137 m 表示地球的長半軸長度；r、θ和λ分別表示計(jì)算點(diǎn)的地心矢徑、緯度和經(jīng)度；和表示正常化的地球重力場模型球諧系數(shù)；n和m分別表示模型的階和次；（cosθ）表示n階m次的勒讓德函數(shù)值；（cosθ） /dθ表示勒讓德函數(shù)的一階導(dǎo)數(shù)。

3 利用地球重力場模型修正水平分量中的低頻誤差

第二節(jié)的分析表明，在重力矢量測量結(jié)果的水平分量中存在低頻誤差，這一低頻誤差所處的頻段和慣導(dǎo)系統(tǒng)水平姿態(tài)誤差的頻譜特性有關(guān)?？梢栽谒椒至康臏y量結(jié)果中將具有誤差的低頻部分去除，然后用相應(yīng)階次的地球重力場模型計(jì)算值來代替具有誤差的低頻（長波）部分，這樣就能對水平分量測量結(jié)果中的低頻誤差進(jìn)行修正。數(shù)據(jù)處理過程如圖3所示。

圖3 數(shù)據(jù)處理流程圖Fig.3 Data processing flow chart

步驟1：利用捷聯(lián)式慣性導(dǎo)航系統(tǒng)中的加速度計(jì)和陀螺數(shù)據(jù)進(jìn)行慣性導(dǎo)航解算；

具體步驟如下：

通過訪談發(fā)現(xiàn)，當(dāng)問到“家鄉(xiāng)與昆明的異同點(diǎn)”時，留學(xué)生經(jīng)常會圍繞空氣、天氣、飲食、建筑、生活等來回答，例如：“昆明空氣很好，就是有點(diǎn)冷”，“昆明的天氣和我的國家一模一樣”，“我的家鄉(xiāng)很擁堵，交通也不是那么方便，在昆明坐公交可以去很多地方”，“昆明很多東西我都吃不慣，我喜歡吃泰國菜”，“昆明的很多街道都很像”等。在一定程度上可以歸納出，留學(xué)生在城市尺度上的地方感主要是在宜人的氣候條件、便捷的生活、不合胃口的飲食、無地方性的城市建筑中形成的。

步驟2：通過GNSS 獲得載體的位置、速度和加速度信息；

步驟3：將GNSS 的位置和速度信息作為外部觀測，結(jié)合慣性導(dǎo)航解算的結(jié)果進(jìn)行組合卡爾曼濾波，估計(jì)出導(dǎo)航誤差參數(shù)并進(jìn)行反饋校正以提高比力的測量精度；

步驟4：利用公式(3)計(jì)算出重力擾動三分量（示意圖如圖4所示，圖4～7 中所示橫坐標(biāo)皆為波長）；

一是由慣導(dǎo)系統(tǒng)測量或計(jì)算得到，包括比力fb和轉(zhuǎn)換矩陣Cbn；

1)施肥要求。畝施雨露水稻專用肥(15-15-15)40公斤作基肥，栽后一周畝施尿素5～6公斤，結(jié)合化學(xué)除草作分蘗肥，畝施尿素7～8公斤、氯化鉀7公斤作穗肥。

無論是西方文明發(fā)展演進(jìn)邏輯，還是中國文明自信演進(jìn)邏輯，都是人類建設(shè)美好家園的路徑選擇。研究兩大文明演進(jìn)的起點(diǎn)、邏輯展開和邏輯歸依，有助于厘清近代以來人類社會的發(fā)展邏輯，有助于世界各國探索和選擇發(fā)展路徑，也有助于構(gòu)建中國文明自信演進(jìn)的話語體系。

圖4 步驟4：根據(jù)公式(3)得到的重力信號Fig.4 Step 4: Gravity information calculated from Eq.(3)

圖5 步驟5：對步驟4 得到的重力信號進(jìn)行低通濾波以去除高頻噪聲Fig.5 Step 5: Low-pass filtering is performed for the gravity information obtained in step 4 to remove the high frequency noise

步驟6：利用高通濾波器將重力測量結(jié)果中受低頻姿態(tài)誤差影響的低頻部分去除（示意圖如圖6所示）。高通濾波器的截止頻率根據(jù)舒勒周期來確定，理想情況下為舒勒周期的倒數(shù)，由于濾波器過渡帶的存在，實(shí)際使用中可適當(dāng)增大高通濾波器的截止頻率（即減小相應(yīng)波長）以保證低頻部分能被充分濾除。

步驟7：利用EGM2008 模型通過公式(11)和(12)計(jì)算重力擾動的兩個水平分量，計(jì)算的階次由步驟6中高通濾波器截止頻率對應(yīng)的空間分辨率來決定，公式如式(13)所示，其中空間分辨率的單位是度。

步驟8：將經(jīng)過步驟5 和步驟6 之后得到的重力水平分量與步驟7 得到的重力水平分量低頻部分EGM2008 計(jì)算值相加即可得到最終的重力水平分量（示意圖如圖7所示），在此結(jié)果中，由姿態(tài)誤差引起的低頻誤差已經(jīng)得到了消除。

圖6 步驟6：對步驟5 得到的重力信號進(jìn)行高通濾波以去除包含姿態(tài)誤差的低頻部分（長波部分）Fig.6 Step 6: High-pass filtering is performed for the gravity information obtained in step 5 to remove the low frequency(long wavelength) part containing attitude error

圖7 步驟8：將經(jīng)過步驟5 和6 得到的重力水平分量加上步驟7 得到的EGM2008 低頻重力水平分量計(jì)算值Fig.7 Step 8: Adding the EGM2008 low-frequency gravity horizontal components calculated by step 7 to the gravity horizontal components by step 5 and 6

4 算 例

2010年5月，國防科技大學(xué)和航空物探遙感中心研制的捷聯(lián)式重力儀SGA01（如圖8所示）在山東的東北部海域進(jìn)行了飛行試驗(yàn)。

圖8 SGA01 示意圖Fig.8 Photograph of SGA01

重力儀所采用的慣性元件包括加表和陀螺，其中陀螺漂移為0.004 ° /h。飛行試驗(yàn)載機(jī)為Cessna208，駕駛員通過自動駕駛儀來控制飛機(jī)的飛行，為了避免強(qiáng)對流氣象條件對測量作業(yè)的影響，飛行試驗(yàn)選擇在天氣情況良好的時候進(jìn)行。測量過程中飛行平穩(wěn)，飛機(jī)保持勻速直線等高度飛行，飛行試驗(yàn)的平均飛行高度為400 m，平均飛行速度為60 m/s。

在這次試驗(yàn)中，SGA01 系統(tǒng)進(jìn)行了重復(fù)測線的飛行以檢驗(yàn)系統(tǒng)的內(nèi)符合精度（如圖9所示），包括6條東西方向的重復(fù)線，每條重復(fù)測線的有效長度為100 km。

圖9 山東試驗(yàn)重復(fù)飛行測線（使用相對坐標(biāo)）Fig.9 Repeat flight lines in the Shandong test(using relative coordinate)

這次試驗(yàn)中，在6 條東西向測線上得到的原始水平分量如圖10 所示，空間分辨率為4.8 km。從圖10中可以看出，由于低頻誤差的存在，各條重復(fù)測線得到的結(jié)果之間存在較大的偏差，內(nèi)符合精度統(tǒng)計(jì)見表1，標(biāo)準(zhǔn)差和均方差差別較大（尤其是北向分量）的原因正是由于低頻誤差的存在。本次試驗(yàn)中所采用的陀螺其漂移為0.004 ° /h，根據(jù)式(10)估計(jì)出的重力擾動水平分量的誤差為14.05 mGal，與表1所示的北向分量均方差基本一致。經(jīng)過步驟5 的低通濾波和步驟6的高通濾波后得到的重力擾動水平分量如圖11 所示，為了將水平分量中的低頻部分充分濾除，高通濾波器的截止頻率選為0.00079 Hz（1/1266 s）。那么，高通濾波所濾除部分的波長大于76 km（1266 s×60 m/s），對應(yīng)EGM2008 的階數(shù)用式(13)計(jì)算為512 階（1°對應(yīng)的空間分辨率大約是108 km）。EGM2008 計(jì)算得到的512 階重力擾動水平分量如圖12 所示。最后將EGM2008 計(jì)算得到的重力擾動低頻部分加上圖11 所示的部分得到如圖13 所示的最終結(jié)果，內(nèi)符合精度見表2。比較表1和表2可以發(fā)現(xiàn)，水平分量中的低頻誤差得到了有效的消除，北向分量的均方差由16.95 mGal 改進(jìn)為2.53 mGal，東向分量的均方差由5.43 mGal 改進(jìn)為2.99 mGal。

人人都說內(nèi)蒙古的羊肉好吃。為什么好吃？其實(shí)也沒什么奧秘，無非內(nèi)蒙古的羊是吃沙蔥的羊，沙蔥本身去膻氣，羊肉固然就少有膻味。

圖10 6 條東西向測線上得到的原始重力擾動水平分量Fig.10 The original horizontal components obtained from the six survey lines in the east-west direction

表1 山東海域試驗(yàn)原始重力水平分量內(nèi)符合精度（mGal）Tab.1 The repeatability of original horizontal components in the Shandong test (mGal)

圖11 6 條東西向測線上經(jīng)過步驟5 和步驟6 后得到的重力擾動水平分量Fig.11 The horizontal components obtained from the six survey lines in the east-west direction after step 5 and step 6

圖12 EGM2008 計(jì)算得到的512 階重力擾動水平分量Fig.12 512-order horizontal components calculated from EGM2008

圖13 6 條東西向測線上經(jīng)過步驟5 和步驟6 后得到的重力擾動水平分量Fig.13 The horizontal components obtained from the six survey lines in the east-west direction after step 5 and step 6

表2 山東海域試驗(yàn)最終重力水平分量內(nèi)符合精度（mGal）Tab.2 The repeatability of final horizontal components in the Shandong test (mGal)

5 結(jié) 論

在不同類型的航空重力儀中，基于捷聯(lián)式慣性導(dǎo)航系統(tǒng)和三軸平臺式慣性導(dǎo)航系統(tǒng)的重力儀不僅可以進(jìn)行重力標(biāo)量測量，還可以進(jìn)行矢量測量。其中基于捷聯(lián)慣導(dǎo)系統(tǒng)的捷聯(lián)式航空重力儀具有結(jié)構(gòu)簡單、體積重量小、操作簡便等優(yōu)點(diǎn)，是重力儀的一個發(fā)展方向。本文對捷聯(lián)式重力矢量測量的誤差機(jī)理進(jìn)行了理論分析，并通過實(shí)測數(shù)據(jù)對水平分量中的誤差進(jìn)行了處理，得出以下結(jié)論：（1）捷聯(lián)式航空重力儀測量得到的重力擾動水平分量主要受到低頻誤差的影響，這一低頻誤差主要是由加速度計(jì)測量誤差和姿態(tài)誤差引起的。

搜集相關(guān)文獻(xiàn)，總結(jié)以往研究成果，提出本研究的基本思路，并在此基礎(chǔ)上構(gòu)建學(xué)生數(shù)學(xué)基本活動經(jīng)驗(yàn)實(shí)踐研究的研究體系。

（2）用本文所述方法對捷聯(lián)式航空重力儀水平分量測量結(jié)果進(jìn)行處理，水平分量內(nèi)符合的均方差得到了明顯的提升，在測量空間分辨率4.8 km 的條件下，北向分量的均方差由16.95 mGal 改進(jìn)為2.53 mGal，東向分量的均方差由5.43 mGal 改進(jìn)為2.99 mGal。