張金煥,李曉斌

(1.中國地質(zhì)大學(xué) (北京)地球物理與信息技術(shù)學(xué)院,北京100083;2.煤炭科學(xué)研究總院,北京100013)

回顧100多年重力勘探的歷史走過的是一條“測量梯度—測量重力場—再測量梯度”的道路。目前測量重力場處于地面、海洋和航空并存的局面,傳統(tǒng)的靜態(tài)地面重力測量因地域、工作量等客觀因素制約,不僅效率低,而且局限于陸地上的部分地區(qū);海洋重力測量是一種低動態(tài)的重力測量方法,且只能在海面工作;衛(wèi)星重力測量開創(chuàng)了重力探測技術(shù)的新時代,但在中高分辨率重力測量中無能為力;而航空重力梯度測量代表了重力勘探領(lǐng)域的最新和最先進(jìn)的技術(shù)。

隨著CPS動態(tài)定位技術(shù)的逐漸成熟,航空重力測量技術(shù)得到迅速發(fā)展,成為在區(qū)域范圍內(nèi)獲取高精度、中高分辨率重力信息的有效手段,可以覆蓋衛(wèi)星重力測量和地面重力測量之間的頻帶。航空重力測量是以飛機(jī)作為運(yùn)載平臺,利用空重力儀 (或加速度計)可以快速測定近地空中重力加速度,對于獲取中高頻重力場信息具有不可替代的作用,具有快速、經(jīng)濟(jì)、靈活和高精度等特點(diǎn),幾乎可到達(dá)任何地區(qū)進(jìn)行作業(yè)。鑒于航空重力梯度測量的高分辨率特點(diǎn),它除了在精細(xì)地質(zhì)勘探領(lǐng)域有著廣泛的應(yīng)用前景以外,也能應(yīng)用于軍事和國防現(xiàn)代化建設(shè)。

1 航空重力測量的基本原理

在經(jīng)典力學(xué)中的慣性參考系中,根據(jù)牛頓第二定律,作用于單位質(zhì)點(diǎn)上的總加速度矢量 fi(稱為比力)、載體運(yùn)動加速度矢量¨ri和引力加速度矢量Gi之間的關(guān)系為:

按照引力加速度和慣性加速度的分離方法,航空重力測量分為航空重力梯度測量和航空重力測量。①重力梯度測量:將兩個共基線的加速度計的輸出求差,以共模方式消除載體運(yùn)動加速度的影響,如果共用基線是旋轉(zhuǎn)穩(wěn)定的,則由該差值觀測量可得到重力梯度分量。②航空重力測量:采用兩個不同的加速度測量系統(tǒng),其中一個系統(tǒng)的輸出含有引力加速度,也就是測量的是比力,而另一個系統(tǒng)輸出的是不含引力的加速度,于是在同一個坐標(biāo)系內(nèi)對兩組加速度輸出進(jìn)行求差,即可消除共有的載體運(yùn)動加速度。剩下的差值中包含引力加速度和傳感器的系統(tǒng)誤差。

按照測量對象的不同,又可以將航空重力測量分為航空標(biāo)量重力測量和航空矢量重力測量。航空標(biāo)量重力測量只需要測量重力擾動矢量垂直分量的大小 (重力異常);航空矢量重力測量需要測量重力擾動矢量所有的三個分量。通常航空重力測量指的就是標(biāo)量重力測量。

目前國際上航空重力測量、航空重力梯度測量均已進(jìn)入實(shí)用階段。由于航空重力梯度測量對重力異常的高頻分量非常敏感,因此在資源勘探等領(lǐng)域具有廣闊的應(yīng)用前景。但是,航空重力梯度測量技術(shù)屬西方國家限制出口的尖端技術(shù),技術(shù)難度大,短期內(nèi)難以掌握。因此,筆者著重討論航空重力標(biāo)量測量系統(tǒng)的研究進(jìn)展。

2 航空重力梯度測量儀器發(fā)展現(xiàn)狀

2.1 旋轉(zhuǎn)加速度計航空重力梯度儀

目前代表性的有兩種——Air-FTGTM 和FALCON。

全張量航空重力梯度測量Air-FTGTM系統(tǒng):Bell Aerospace公司 (現(xiàn)并入Lockheed Martin公司)成功研制的FTG是源于1970年美國空軍的發(fā)展計劃,出于導(dǎo)航和導(dǎo)彈發(fā)射的需要美國空軍投入十億多美元。研制計劃正式始于1972年,Bell是從當(dāng)時四個競爭對手中脫穎而出的,它的研制目的是進(jìn)行全球重力場的測繪。于1987年5月成功進(jìn)行了重力梯度的航空測量。隨著硬件和軟件的改進(jìn),精度在不斷提高。2002年研制成功全張量航空重力梯度儀,已經(jīng)獲得了許多商業(yè)合同,截止到2004年已在北美和非洲等地采集了超過60,000km的重力梯度數(shù)據(jù),還成功的應(yīng)用于潛艇導(dǎo)航等其他軍事領(lǐng)域。另外,英國ARKeX公司將全張量重力梯度儀 (FTG)做了特殊改進(jìn)形成FTGeX系統(tǒng),并于2005年在北美測試后進(jìn)行了首次商業(yè)運(yùn)行。FTGeX系統(tǒng)類似于Bell GeoSpace公司的FTG技術(shù),經(jīng)改進(jìn)適應(yīng)一些特殊要求。

部分梯度張量航空重力梯度測量FALCON系統(tǒng):目前已投入商業(yè)使用的還有澳大利亞BHP公司研制的部分梯度張量測量系統(tǒng)FALCON。該公司于20世紀(jì)80年代末引進(jìn)美國Bell公司 GGI技術(shù),經(jīng)過十年的研究開發(fā),研制成功FALCON航空重力梯度儀。于1999年正式投入使用,它測量時可同時附帶進(jìn)行航磁和航放測量。FALCON已在澳大利亞、北美、南美和非洲共飛行了28萬測線km,目的是普查油氣和各種其他礦產(chǎn)資源,取得了良好的效果。FALCON是美國出口管制產(chǎn)品。

2.2 航空超導(dǎo)重力梯度儀 (SGG)

M aryland大學(xué)超導(dǎo)重力梯度儀:超導(dǎo)重力梯度儀研制的權(quán)威機(jī)構(gòu)為美國的M aryland大學(xué)物理系,該大學(xué)長期以來一直從事重力梯度儀的開發(fā)和研究,并取得了許多有價值的研究成果。20世紀(jì)80年代初 Maryland大學(xué)研制出了靈敏度為0.3～0.7E/的單軸超導(dǎo)重力梯度儀實(shí)驗(yàn)室樣機(jī)。2002年M aryland大學(xué)已經(jīng)研制出單軸超導(dǎo)重力梯度儀實(shí)驗(yàn)室樣機(jī)并進(jìn)行了試驗(yàn),其靈敏度提高到了0.02E/,比常溫下重力梯度儀提高3個數(shù)量級,目前馬里蘭大學(xué)在研的Ⅲ型超導(dǎo)重力梯度儀預(yù)期的靈敏度為10～4E/。

英國 ARKeX公司超導(dǎo)重力梯度儀:英國ARKeX公司是由著名的英國牛津儀器公司和另一家地球物理公司的一些人員新組成的一家以生產(chǎn)重力梯度儀為主的新企業(yè)。牛津儀器公司為歐洲航天局設(shè)計制造超導(dǎo)重力梯度儀,但未按期完成任務(wù),因而未能用在衛(wèi)星上。2004年3月歐洲航天局繼續(xù)資助該項(xiàng)目。該公司改進(jìn)的新一代航空重力梯度儀具有卓越的穩(wěn)定性、低噪聲和高分辨率等性能。公司在航空重力梯度儀方面的研究取得了突破性的進(jìn)展,公司的最終目標(biāo)是生產(chǎn)經(jīng)營世界上首臺超導(dǎo)勘探重力梯度儀 EGG(Exp loration Gravity Gradiometer)。EGG的設(shè)計可放六個加速度計模塊,用來測量全部的梯度張量分量,但第一臺EGG只用了兩個加速度計模塊,只測量Γzz。該儀器計劃在2006年第3季度進(jìn)行商業(yè)勘探。用改進(jìn)后的 EGG系統(tǒng)進(jìn)行試驗(yàn),其靈敏度、分辨率和穩(wěn)定性空前卓越,適應(yīng)了動態(tài)航空重力勘探環(huán)境。其設(shè)計指標(biāo)為:靈敏度1E/,帶寬為200m～60km。EGG的組成包括:1296個超導(dǎo)接頭,94個超導(dǎo)變壓器,151個熱開關(guān),60個超導(dǎo)鍍錫卷板,14個SQU IDs,2個鈮檢測質(zhì)量,1個低溫恒溫器,總重量為500kg。

澳大利亞西澳大學(xué) (UWA)超導(dǎo)重力梯度儀:澳大利亞西澳大學(xué) (UWA)研制的低溫超導(dǎo)重力梯度儀是針對資源勘探而按照適宜帶寬設(shè)計的。該重力梯度儀主要依靠低溫和超導(dǎo)技術(shù),安裝在輕型飛機(jī)上,專用于低空(100m)飛行測量。它主要用于動基座低空范圍內(nèi)的地質(zhì)勘探,特別針對礦藏資源的勘探。這種重力場測量技術(shù)在很大程度上克服了傳統(tǒng)重力計測量重力的限制,可以直接測量重力場的梯度變化(而不僅僅是重力值)。這種信號盡管十分微弱但卻蘊(yùn)涵了豐富的地下信息,因此這種儀器無論是應(yīng)用前景還是技術(shù)難點(diǎn)都將是非常巨大的。

加拿大 GEDEX公司超導(dǎo)重力梯度儀:加拿大GEDEX公司開發(fā)的AGG-OQR重力梯度儀是綜合采用了M aryland大學(xué)、UWA大學(xué)和加拿大空間公司的技術(shù),也將于近期投入商業(yè)運(yùn)營。

2.3 其他航空重力梯度儀

美、俄等國的科學(xué)家還在積極研究多種重力梯度測量方案。80年代后期俄羅斯專家研制出了精度為0.1E的旋轉(zhuǎn)加速度計重力梯度儀實(shí)驗(yàn)室樣機(jī)。此外利用原子干涉現(xiàn)象測量重力的方案特別值得注意。諾貝爾物理學(xué)獎 (1997)獲得者、美籍華裔物理學(xué)家朱棣文領(lǐng)導(dǎo)的小組,利用原子噴泉構(gòu)成的原子干涉儀來測量g,測定g的分辨率己達(dá)到△g/g=10～10,即可分辨0.1μGal的重力變化。美國許多科學(xué)家正在努力使這種測量重力的設(shè)備小型化,并構(gòu)成梯度裝置使之適合在運(yùn)動平臺上進(jìn)行觀測。

表1為航空重力梯度測量系統(tǒng)一覽情況。

表1 航空重力梯度測量系統(tǒng)一覽

3 國內(nèi)航空重力測量發(fā)展現(xiàn)狀

由總參西安測繪研究所等單位承擔(dān)的我國首臺航空重力測量系統(tǒng)CHAGS(Chinese -airbo rne gravimeter system)的研究始于 1995年,引進(jìn)LCR(Lacoste&Rom berg) Ⅱ型航空重力儀(L&R126),自行研發(fā)了航空重力數(shù)據(jù)處理軟件,并于2002年11月通過鑒定,歷時7年。2002年3～4月CHAGS航空重力測量系統(tǒng)在山西大同地區(qū)進(jìn)行了3次飛行試驗(yàn),飛行高度2800～3400m、飛行速度400km/h,測線間距5～10km,測定山區(qū)平均重力異常的內(nèi)符合精度為3.32m Gal,異常半波長分辨率為10km。與地面重力測量數(shù)據(jù)的對比,外符合精度優(yōu)于5m Gal。目前該系統(tǒng)已經(jīng)基本可以滿足大地水準(zhǔn)面測量等應(yīng)用的要求,但是還不能達(dá)到地質(zhì)調(diào)查、資源勘探等高精度、高分辨應(yīng)用的要求。

到目前為止,國內(nèi)發(fā)展起來的航空重力測量主要依靠的是進(jìn)口的方式。而用于地球物理勘探的航空重力測量技術(shù)還處于起步階段,許多技術(shù)問題亟待解決。最主要的是航空重力儀制造相關(guān)技術(shù),難度較大,需國內(nèi)相關(guān)各方協(xié)作,爭取早日獲得突破。在航空重力梯度測量儀器方面,國內(nèi)相關(guān)研究基本尚處于空白領(lǐng)域。

4 航空重力梯度儀需解決的關(guān)鍵技術(shù)

由于重力信號很微弱,就對傳感器、信號轉(zhuǎn)換、信號放大,以及環(huán)境噪聲隔離、加速度計溫度控制等現(xiàn)有技術(shù)水平提出嚴(yán)峻的挑戰(zhàn),許多因素都應(yīng)予以考慮。

1)加速度計性能的影響。加速度計是系統(tǒng)中的核心關(guān)鍵部件之一,其性能的重要性是不言而喻的。被測信號的微弱性決定了所使用的加速度計必須具有更高的性能要求,如誤差模型需要更加精確化,其動態(tài)誤差系數(shù)要可測,并能夠快速補(bǔ)償;加速度計的熱噪聲需保持在一個相當(dāng)?shù)偷乃?否則會嚴(yán)重影響測量精度;多個加速度計標(biāo)度因數(shù)的匹配也是一個有待解決的問題。

2)平臺性能的影響。穩(wěn)定平臺是能否有效適應(yīng)航空環(huán)境的關(guān)鍵,平臺的加工、安裝定位精度對重力梯度儀系統(tǒng)是至關(guān)重要的。但是由于受到現(xiàn)有加工水平等因素的影響,要求從其他方面想方法來加強(qiáng)對平臺的運(yùn)動、振動和溫度等的控制,力圖減小平臺對系統(tǒng)整體性能的影響。同時三個平臺的穩(wěn)定同步旋轉(zhuǎn)也是一個比較大的難點(diǎn)。

3)安裝誤差的影響。重力梯度儀是非常精密的儀器,任何微小的安裝誤差都可能其對性能造成嚴(yán)重的影響,其中需要重點(diǎn)關(guān)注的安裝誤差主要包括:平臺的安裝定位誤差;加速度計角度安裝誤差;加速度計徑向安裝誤差。需要探索研究新型的安裝定位方法以保證系統(tǒng)精度。

4)電磁干擾和屏蔽技術(shù)。系統(tǒng)中有用的信號主要是二倍頻信號,可以采用數(shù)字濾波技術(shù)來去除噪聲干擾等,但同時由于重力梯度信號非常微弱,外界的電磁信號對系統(tǒng)的影響也就不可忽視,需要采用電磁屏蔽措施對系統(tǒng)加以防護(hù),使系統(tǒng)具有較高的電磁兼容性。

5)高速采集和處理。加速度計的輸出信號比較小,要抑制線加速度和角加速度的共模干擾,則必須實(shí)現(xiàn)多個加速度計的并行、高速和高精度數(shù)據(jù)采集。要將系統(tǒng)用于機(jī)載實(shí)時系統(tǒng)中進(jìn)行系統(tǒng)誤差修正,則系統(tǒng)必須實(shí)時解算重力梯度和重力向量值,同時換算出重力異常和垂線偏差,這些都要求旋轉(zhuǎn)加速度計重力梯度儀系統(tǒng)中必須采用高速的專用計算機(jī)進(jìn)行數(shù)據(jù)的采集和處理。有效的相關(guān)方法的研究也是很重要的。

航空重力梯度儀具有極高的研究和應(yīng)用價值,現(xiàn)實(shí)應(yīng)用還是有難度的,需要在加大基礎(chǔ)理論研究力度,加強(qiáng)先進(jìn)數(shù)據(jù)處理算法、解算方法和其他相關(guān)應(yīng)用程序等研究的基礎(chǔ)上,在各個環(huán)節(jié)上都求尖求精。

5 總結(jié)

正在使用和處于研制階段的重力梯度儀有20余種,針對航空應(yīng)用的還不多,如果關(guān)鍵技術(shù)問題得到解決,航空重力梯度儀就發(fā)揮出其無法比擬的優(yōu)勢。在這方面,我們又走在了部分西方國家的后面,那么買進(jìn)一臺相應(yīng)的設(shè)備就更是困難,而且價格肯定驚人。從各個方面講,目前緊迫和最有效的方法就是集中資金和力量搞自主研發(fā),盡快掌握這門我國急需的有用技術(shù)。

此外,開展地球重力梯度的研究,歷來是大地測量學(xué)的熱點(diǎn)和難點(diǎn)之一,因?yàn)樗坏从车厍虮韺蛹皟?nèi)部的密度分布和物質(zhì)運(yùn)動狀態(tài),而且在地下資源勘探、海洋、軍事應(yīng)用和環(huán)境科學(xué)等領(lǐng)域中有著重要的應(yīng)用。傳統(tǒng)的測量方法,需要消耗大量的時間和人力,而且不能滿足諸多實(shí)時系統(tǒng)的要求,性能也不佳。比較而言,航空重力梯度測量儀器省時省力,可以同時測量重力梯度的各個分量,進(jìn)而獲得準(zhǔn)確的重力值和垂線偏差,還可以應(yīng)用于實(shí)時場合進(jìn)行相關(guān)測量與誤差補(bǔ)償,具有廣闊的研究價值和應(yīng)用前景。

[1] 趙德軍.航空矢量重力測量的理論與方法 [D].鄭州:信息工程大學(xué)測繪學(xué)院,2005.

[2] 夏哲仁,石磐,孫中苗.航空重力測量系統(tǒng)CHAGS[J].測繪學(xué)報,2004,33(3):216-220.

[3] 孫中苗.航空重力測量理論、方法及應(yīng)用研究 [D].鄭州:信息工程大學(xué)測繪學(xué)院,2004.

[4] 周衍柏.理論力學(xué)教程 [M]2版.北京:高等教育出版社,1985.

[5] 董緒榮,張守信,華仲春.GPS/INS組合導(dǎo)航定位及其應(yīng)用 [M].長沙:國防科技大學(xué)出版社,1998.

[6] 楊元喜.抗差估計理論及其應(yīng)用 [M].北京:八一出版社,1993.

[7] 歐吉坤.誤差理論若干問題研究——兼論大氣對 GPS測量的影響 [R].中國科學(xué)院測量與地球物理研究所,1994.

[8] 胡明城,魯福.現(xiàn)代大地測量學(xué):下冊 [M].北京:測繪出版社,1994.

[9] 許其鳳.GPS衛(wèi)星導(dǎo)航與精密定位 [M].北京:解放軍出版社,1994.

[10] 孫廷萱,肖慶達(dá),萬榮頤.GPS譯文集 [C].北京:中國計量出版社,1989.

[11] 劉大杰,施一民,郭靜君.全球定位系統(tǒng) (GPS)的原理與數(shù)據(jù)處理 [M].上海:同濟(jì)大學(xué)出版社,1999.

[12] William R G.An Historical Review of Airborne Gravity[J].The Leading Edge,1998(1):113-116.

[13] 熊盛青.我國航空重磁勘探技術(shù)現(xiàn)狀與發(fā)展趨勢 [J].地球物理學(xué)進(jìn)展,2009,24(1):113-117.

[14] 孫中苗,石磐,夏哲仁,等.利用GPS和數(shù)字濾波技術(shù)確定航空重力中的垂直加速度 [J].測繪學(xué)報,2004,33(2):110-115.

[15] 孫中苗,夏哲仁.FIR低通差分器的設(shè)計及其在航空重力測量中的應(yīng)用 [J].地球物理學(xué)報,2000,43(6):850-855.

[16] Bruton A M.Imp roving the Accuracy and Resolution of SINS/DGPS Airborne Gravimetry[D].Calgary:Department of Geomatics Engineering at University of Calgary,2001.

[17] JUAN GILBERTO SERPAS.Local and regional geoid determi-nation from vector airbo rne gravimetry[D].Columbus:The O-hio StateUniversity,2003.

[18] DennisD,McCarthy.IERSTechnicalNote21[R].Centra Bureaau of IERS-Observatoire de Paris,1996.

[19] GrafarendEW,Krumm F,OkekeF.Curvi linear geodet datum transformations[J].Zfv,1995(7):334-350.

[20] Gravity Diamonds Limited,FALCON Summary[EB/OL].http://www.gravitydiamonds.com.au/BHPBFalcon-Brochure.pdf

[21] Principle of Gravity Gradiometry [EB/OL].http://www. physics. umd. edu/GRE/Gravity _Gradiometry.pdf

[22] Pullin J.,Meters of Gravity[EB/OL].http://www.phys.lsu.edu/mog/pdf/mog9.pdf

[23] Superconducting Gravity Gradiometer for Moving-Base Applications[EB/OL].http://www.physics.umd.edu/GRE/Moving_base_SGG.pdf

[24] Superconducting Gravity Gradiometers(SGGs)[EB/OL].http://www. physics. umd. edu/GRE/NASA _SGG.pdf

[26] William R G.An Historical Review of Airborne Gravity[J].The Leading Edge,1998(1):113-116.