張永旺 趙麗娟 趙婧 蘇鵬飛

摘要：重力固體潮改正是重力數(shù)據(jù)處理的基礎(chǔ)，在重力工作中重力儀的儀器試驗(yàn)、重力基點(diǎn)網(wǎng)聯(lián)測(cè)、重力測(cè)點(diǎn)絕對(duì)重力值計(jì)算等主要重力數(shù)據(jù)處理中，重力固體潮理論值的計(jì)算都必不可少。但重力固體潮的計(jì)算公式繁雜，為提高計(jì)算效率，簡(jiǎn)化計(jì)算過(guò)程，筆者應(yīng)用EXCEL中的Visual Basic編輯器，通過(guò)VBA代碼實(shí)現(xiàn)了重力固體潮理論值的一鍵計(jì)算。

關(guān)鍵詞：重力固體潮；VBA；重力計(jì)算

Abstract： the theory value of Gravity solid earth tides is the basis of gravity data processing， in the gravity test instrument， gravity basepointsnet conjunction， absolute gravity value of observation point and so on， the theory value of Gravity solid earth tides is necessary of gravity data-processing. But Gravity solid earth tides correction is complicated. In order to improve the computational efficiency， simplify the calculation process， the author take full advantage of VBA editor in the EXCEL， through VBA code calculating the theory value of Gravity solid earth tides by one key.

Key words： Gravity solid earth tides；VBA；Gravity Calculation

前言

隨著礦產(chǎn)資源的勘查與開(kāi)發(fā)日趨向深部發(fā)展，深入了解大地構(gòu)造，劃分?jǐn)嗔褬?gòu)造及其控礦特征，圈定隱伏巖體，對(duì)于開(kāi)展深部找礦具有決定性的作用，而重力工作對(duì)于解決以上地質(zhì)問(wèn)題具有明顯優(yōu)于其他物化探方法的優(yōu)勢(shì)，所以重力工作對(duì)于進(jìn)一步推動(dòng)地質(zhì)礦產(chǎn)勘查的轉(zhuǎn)型和發(fā)展都具有重要意義。

重力工作中數(shù)據(jù)處理的真實(shí)準(zhǔn)確是開(kāi)展解釋推斷的基礎(chǔ)，而重力固體潮理論值的計(jì)算又是重力數(shù)據(jù)處理的前提，所以準(zhǔn)確無(wú)誤的獲得重力固體潮的理論值顯得尤為重要。

目前，重力固體潮理論值的計(jì)算有幾款商業(yè)軟件可以實(shí)現(xiàn)，但不同軟件輸入格式多樣，其計(jì)算結(jié)果應(yīng)用于下級(jí)計(jì)算步驟繁雜，不利于資料的二次利用。通過(guò)VBA實(shí)現(xiàn)重力固體潮的計(jì)算可以簡(jiǎn)化重力資料處理，提高工作效率。希望文章對(duì)重力工作者有所幫助，同時(shí)也歡迎廣大重力工作者批評(píng)指正。

1. 概念

1.1 重力固體潮

由于地球和太陽(yáng)及月亮的相對(duì)位置不斷變化而引起的重力變化稱為重力固體潮。

1.2 重力固體潮理論值

月球和太陽(yáng)對(duì)地球上某一點(diǎn)作用值之和即為重力固體潮的理論值。

2. 格式要求

數(shù)據(jù)格式依次為點(diǎn)號(hào)、日期（年-月-日）、時(shí)間（60進(jìn)制，h：mm：ss）、經(jīng)度（60進(jìn)制，DDDMMSS.S）、緯度（60進(jìn)制，DDMMSS.S）

運(yùn)行VBA代碼后會(huì)在F列顯示點(diǎn)對(duì)應(yīng)的重力固體潮理論值。

3. VBA代碼及解釋說(shuō)明

4. VBA代碼的應(yīng)用

實(shí)際使用時(shí)，新建一個(gè)Excel表格，在Sheet1中輸入數(shù)據(jù)，然后在工具選項(xiàng)卡中打開(kāi)宏-Visual Basic編輯器，在代碼窗口中，粘貼上述代碼，點(diǎn)擊運(yùn)行（如圖1紅圈所示）按鈕即可在數(shù)據(jù)尾列獲得重力固體潮理論值（單位毫伽）。

5. 算例與分析

為確保VBA代碼的準(zhǔn)確性和可靠性，筆者針對(duì)性的選擇以下特殊情況的數(shù)據(jù)（見(jiàn)表2）進(jìn)行驗(yàn)算，并與“區(qū)域重力信息系統(tǒng)RGIS2006”、“跨平臺(tái)金維地學(xué)信息處理研究應(yīng)用系統(tǒng)GeoIPASV3.2”計(jì)算成果進(jìn)行對(duì)比。

從表2可知，由于計(jì)算過(guò)程中截取位數(shù)的差別，筆者編寫(xiě)的VBA代碼與RGIS2006的計(jì)算成果最大偏差0.002毫伽，與GeoIPAS V3.2最大偏差0毫伽，證明在不同特殊情況下，筆者編寫(xiě)的VBA代碼計(jì)算成果準(zhǔn)確可靠。經(jīng)過(guò)大量檢查計(jì)算（此處不一一列舉），VBA代碼的計(jì)算成果與RGIS 2006的最大偏差僅0.003毫伽，與GeoIPASV3.2基本沒(méi)有偏差。

6. 結(jié)束語(yǔ)

利用VBA代碼進(jìn)行重力固體潮的改正不僅方便快捷，便于進(jìn)行數(shù)據(jù)的下級(jí)計(jì)算，而且精度可靠，能夠滿足大比例尺重力工作，具有較高的應(yīng)用推廣價(jià)值，同時(shí)歡迎廣大重力工作者批評(píng)指正。

參考文獻(xiàn)：

[1] DZ/T0082-2006，區(qū)域重力調(diào)查規(guī)范[S].

[2] 馬玄龍，劉永生.重力固體潮在Excel工作表中的計(jì)算方法[A].資源環(huán)境與工程.2010，24（2）：203-206.

[3] 宋偉，吳建國(guó)等.中文Visual Basic6.0編程基礎(chǔ)[M].北京：清華大學(xué)出版社，1999.